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It  is  necessary  to  know  how  much  protection  structures 
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radioactive  material.  This  information  is  well  known  for 
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fallout.  This  research  transfers  the  knowledge  gained  from 
weapon  fallout  work  to  protection  factors  for  any  gamma 
spectrum. 

Point  kernel  models  were  developed  for  both  fallout  and 
cloud  sources.  That  development  included  a  method  of 
accurately  combining  buildup  factors  in  multi-region 
problems  over  wide  ranges  of  energy  and  photon  mean  free 
path,  and  a  method  for  calculating  the  effect  of  ground 
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The  results  were  reported  for  six  spectra  as  well  as 
discrete  energies  from  15  KeV  to  15  MeV.  The  structures 
used  as  examples  include  small  wood  frame  and  large  brick 
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1-hr  weapon  fallout.  However  there  are  significant 
differences  for  other  spectra,  such  as  that  from  Three  Mile 
Island.  The  effects  of  varying  building  size  are  reported 
as  well  as  the  relative  importance  of  both  cloud  and  fallout 
sources  that  infiltrate  structures. 
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CHAPTER  I 


INTRODUCTION 


1 . 1  Qcigi n_sf _ih#_Prgb!f m 

In  the  event  of  an  unplanned  release  of  radioactive 
material,  it  is  necessary  to  know  the  doses  the  public  could 
receive  in  order  to  make  decisions  that  minimize  risk  to 
that  public.  The  research  reported  here  is  as  outgrowth  of 
preliminary  work  done  to  determine  what  doses  the  public 
might  receive  if  they  try  to  evacuate  or  seek  shelter. 

Figure  1  is  the  flow  chart  of  a  conceptual  program 
designed  to  determine  which  alternative  provides  the  lowest 
dose  to  Individuals.  The  alternatives  arei  taking  no 
action)  seeking  shelter  in  homes)  seeking  shelter  in  other 
buildings,  such  as  factories  or  schools)  or  evacuating.  The 
input  includes  weather  data  and  the  chemical  and  isotcpic 
nature  of  the  release.  A  dispersion  and  deposition  model 
determines  the  real  time  and  projected  concentrations  in  the 
radioactive  cloud  and  fallout.  The  infiltration  models 
determine  how  much  material  gets  inside  the  structures.  The 
structure  shielding  (and  vehicle  shielding)  models  determine 
how  much  the  unprotected  dose  is  reduced  by  available 
structures.  The  dose  models  compute  the  doses  from  the 
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■fallout  and  tha  cloud  sources.  An  avacuatlon  modal 
aatimataa  how  fast  and  in  what  diraction  tha  population  is 
expected  to  exit  in  tha  avent  evacuation  is  ordered.  A  time 
use  model  estimates  how  much  time  an  unwarned  papulation 
would  spand  indoors,  outdoors,  or  in  vehicles.  The  result 
of  tha  program  is  three  dose  estimates  for  the  populations 
in  each  sector  around  tha  accident  and  at  different  ranges 
from  the  accident.  These  doses  arai  a  sheltered  dose,  a  "do 
nothing  dose”  for  the  option  of  not  warning  the  population, 
and  a  dose  received  while  evacuating.  With  this 
information,  an  official  can  make  informed  decisions  as  to 
what  the  public  should  be  told.  Ideally  this  program  would 
be  contained  in  mini -computers  available  near  the  release 
site  and  would  provide  results  in  seconds  or  at  most  a  few 
minutes.  The  research  reported  here  was  done  to  provide 
data  for  the  structure  shielding  subroutines  of  su_h  a 
program.  However,  the  bibliography  included  in  this  paper 
is  a  survey  of  literature  applicable  to  the  total  problem 
concerning  which  action  alternative  provides  the  lowest 
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i  >  2  BsiMCS&_Qb Jisiiy* 

The  objective  of  this  research  it  to  develop  a  model 
for  tha  protaction  provldad  by  structures,  much  as  home*, 
against  ralaasas  of  airborna  radioactiva  material  as  a 
function  of  tha  gamma  energy  spectrum  of  tha  released 
material.  Doses  racaivad  through  tha  food  chain  or 
inhalation  are  net  considered. 

It  is  not  nacassary  to  do  a  detailed  analysis  of  each 
structure.  Tha  results  of  detailed  studies,  both 
mathematical  and  empirical  have  bean  raportad  for  a  few 
spectra  and  building  types.  Tharsfora  tha  purpose  of  this 
research  is  not  to  calculate  protection  factors,  but  rather 
to  show  how  they  vary  with  energy. 
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FIGURE  1  Elfi»_Gb«Ci-£fiE-*-BfttiiSi£UXf-B9ll«19.fi£&l0D 
BlElllSQ-BOfrlx 


s 


1  •  3  §2£&grsuQd_cf _&d9_Sbit IdiQg.Ecsklto) 

During  tha  1930 's  and  aarly  1960 'a  a  greet  daal  of  work 
waa  dona  to  predict  tha  protection  one  would  racaiva  from 
weapon  fallout  by  aaaking  ahaltar  in  a  home,  basement, 
factory  or  bomb  shelter.  This  work  was  the  result  of  the 
Federal  civil  defense  policy  of  the  time  Ell.  The  result  of 
this  effort  was  the  development  of  the  Office  of  Civil 
Defense  (OCD)  Standard  Method  in  several  versions  meant 
primarily  for  architects  or  for  surveying  existing  buildings 
C2,3,4,3,63.  Very  little  additional  work  was  done  until  the 
mid  1970s.  The  publishing  of  the  Reactor  Safety  Study  C7D 
prompted  a  new  interest  in  structure  shielding.  In  1973,  Z. 
6.  Burson  and  A.  E.  Profio  published  SfcELUEtUC9_Stli.ll.lli.Qg 
fn9<D_Clfiud_«Dd_EiiIgui_fi*n!M_3iy_SfiycEsi_fac_flass59iDfl_ibff 
CgQS»gytOE*S_gf_BMEigC_AESldiQi*  For  fallout 

protection  data,  they  relied  on  the  previous  bomb  fallout 
studies,  and  for  cloud  source  protection  data,  they 
developed  a  point  kernel  integration  technique.  Their  data 
were  used  by  G.  H.  Anno  and  M.  A.  Dors  in  "The  Effectiveness 
of  Sheltering  as  a  Protective  Action  Against  Nuclear 
Accidents  Involving  Gaseous  Releases",  which  was  published 
by  the  EPA  in  197B  C93.  That  report  raised  the  issue  of 
infiltration  of  radioactive  gases  into  the  structures. 

Their  report  was  further  amplified  by  work  done  at  Sandia 
Laboratories  in  1977  and  1978  C10,113.  The  Three  Mile 


6 


Island  accident  promptsd  a  new  public  interest  in  seeking 
protection  -from  radioactive  releasee.  In  1980,  the  National 
Bureau  o f  Standards  published  §£CU££U£t_Sbi9ldiQg_AggiQf£ 
ElIlfiyi-SiffiJM-BiYS-ECSffl-NyElsaC-BsiBQBiisnS  C1D.  This  work 
is  both  a  summary  and  a  complete  update  of  the  previous 
fallout  shelter  work,  but  it  deals  only  with  fallout  from  a 
nuclear  weapon.  Nearly  all  the  information  in  the  open 
literature  to  date,  including  the  works  just  mentioned,  rely 
on  gamma  energy  spectra  either  derived  from  bomb  fallout,  or 
from  scenarios  involving  the  release  of  a  great  amount  of 
nuclear  reactor  core  materials.  As  the  TM  and  8L-1 
accidents  pointed  out,  it  is  far  more  likely  that  only  noble 
gases  (which  have  a  much  lower  energy  spectra  than  weapon 
fallout)  and  a  limited  amount  of  volatile  materials,  such  as 
iodine,  will  escape  even  from  a  severe  accident  £12,133. 
Radioactive  releases  are  also  possible  from  other 
activities,  such  as  transportation  of  spent  reactor  fuel, 
and  manufacturing  accidents  involving  medical  sources.  I£ts 

c*lMi«*_wb*iiY»c_tb9_c*l«i»*sJ-*B»Eicym-.ffliflbt_b»A 


V  •  .  •  .  V  V  v  *„•  •  •  \-  v  V  v  V '  -  •.>  v  v  W  V  >V  V  *„•  V  V  V.v 

'-".v  v* * -"*  .*•  A  .*■  .**l*'^  -**  l'-  .**  .*  /’  /•  /•  <  v'.-  v 


7 


1 . 4  QcaiDilf iifiO 

This  paper  consists  of  a  discussion  of  ths  rational* 
behind  th*  s*l*ction  of  a  point  k*rn*l  model  for  both  th* 
fallout  and  cloud  sourc*  conditions  and  a  description  of  th* 
component  buildup  factor,  ground  roughness  and  geometry 
models  that  were  used.  The  results  of  the  models  are 
related  to  structures  such  as  homes  and  vehicles  and  to  the 
historical  cases  of  TM1-2  and  SL-1. 

The  three  appendixes  contain  descriptions  and  FORTRAN 
listings  of  the  programs  developed  to  provide  the  required 
shielding  data. 

The  protection  factors  reported  here  are  called  Dose 
reduction  factors  (DRFs).  This  term  is  adopted  to  prevent 
confusion  with  other  terms  such  as  shielding  factors  (SFs) 
and  Protection  Factors  <PFs>  which  ars  defined  differently 
by  different  researchers.  In  the  strict  sense,  a  DRF  should 
be  defined  as  the  ratio  of  the  actual  dose  received  to  the 
dose  that  would  be  received  if  there  was  no  protection  at 
all.  Unfortunately,  the  relationship  between  dose  and  gamma 
energy  for  a  standard  individual  is  not  easily  defined  for 
low  energy  gamma  rays.  The  unavailability  of  this  data 
required  the  use  of  exposure  instead  of  dose  in  calculating 
the  DRFs  reported  here.  Because  the  DRF  is  a  ratio  of  two 
exposures  (or  doses)  taken  at  the  the  same  gamma  energies, 
the  errors  in  using  exposure  instead  of  dose  cancel  each 
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other.  The  DRFt  ae  ueed  in  this  paper,  ie  defined  as  the 
ratio  of  the  protected  exposure  to  the  theoretical 
unprotected  exposure. 


CHAPTER  II 


STRUCTURE  SHIELDING 

2.  1  CBdBl_S*llEiiSD 

The  principle  objective  of  this  dissertation  is  to 
predict  how  tha  exposure  from  tourcaa  aaaociatad  with  a 
ralaaaa  of  radioactiva  matarial  variaa  am  a  function  of  tha 
anargy  of  tha  amittad  gamma  rays.  Tha  ability  to  pradict 
gamma  flux  as  a  function  of  anargy  then  bacama  a  prime 
factor  in  choosing  a  mathematical  modal.  Tha  ability  to 
easily  identify  tha  components  of  tha  model  and  ralata  them 
to  tha  physical  world  Mas  also  a  prime  salaction  criteria. 

Shielding  data  have  bean  calculated  or  experimentally 
found  for  a  few  shielding  geometries,  and  for  a  vary  limited 
number  of  energies  C143.  Tha  energies  used  in  tha  available 
literature  are  associated  Mith  weapon  fallout  or  the 
postulated  PWR-2  accident  identified  in  the  Reactor  Safety 
Study  (RSS)  CIS].  C  The  postulated  PWR-2  accident  Hincludes 
failure  of  the  cooling  systems,  and  core  meltdown  concurrent 
with  a  loss  of  containment  spray  and  heat  removal  systems. 
Failure  of  the  containment  barrier  occurs  through 
overpressure  causing  a  substantial  fraction  of  the 
containment  atmosphere  to  be  released  in  a  'puff'"  C16D.3 
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While  soma  good  data  ai  j  aval  labia  for  thaaa  high  anargy 
•pactra,  almost  nothing  la  aval labia  for  othar  anargias. 

Thus  a  mathematical  modal  only  naada  to  dat ar mi na  how  tha 
shielding  factors  vary  with  energy,  and  than  relate  them  to 
tha  known  values  for  tha  bomb  fallout  spectra. 

Three  methods  have  baan  developed  that  can  calculate 
tha  shi si ding  provided  by  a  structural  moments, 

Monte-Carlo,  and  point  kernel  integration.  Tha  method  of 
moments  is  tha  basis  for  tha  "standard"  method  C23.  Tha 
standard  mathod  is  vary  good  for  calculating  tha  protection 
provided  by  a  specific  structure  fraa  a  specific  anargy 
spectra,  usually  tha  1.12  hour  bomb  fallout  spectra. 

However,  each  new  anargy  spectra  requires  the  development  of 
a  set  of  tables,  nomographs,  and  curves  -  making  it  very 
difficult  to  use  the  technique  over  a  broad  energy  range. 

The  Monte-Carlo  method  Is,  in  theory,  an  exact  solution 
to  the  radiation  transport  equation  and  can  be  used  for  any 
geometry  or  energy  distribution.  It  has  the  advantage  of 
being  able  to  find  nearly  exact  answers  to  specific 
problems,  however  it  is  cumbersome  to  use,  requires  a  great 
deal  of  computer  time,  and  its  impl  seen  tat  ion  is  difficult 
to  relate  back  to  the  physical  problem.  The  Monte-Carlo 
method  is  unnecessarily  complex  for  the  problem  at  hand. 

The  point  kernel  method  has  the  advantage  of  being 
easily  related  to  simple  geometries.  While  it  cannot 
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provide  exact  solutions  to  real  problems,  it  can  mak® 
raasonabl®  estimates  and  is  aasily  adaptad  to  varying 
anargias.  Tha  point  karnal  mathod  can  ba  usad  to  provlda 
tha  trans-forms  nacassary  to  gat  from  known  shielding  values 
to  values  at  other  anargias.  For  these  reasons*  tha  point 
karnal  mathod  was  chosen  as  tha  math  modal  in  this 
dissertation. 

Evan  with  this  Beamingly  simple  choice,  tha  mathod 
required  a  great  deal  of  davalopmant  in  order  to  provide 
reasonably  accurate  estimates.  First,  buildup  factors  that 
are  accurate  over  tha  nacassary  rangas  of  energy  and  mean 
frae  path  had  to  ba  found.  Second,  an  accurate  mathod  of 
combining  the  buildup  factors  in  multi-region  problems  over 
a  wide  energy  range  had  to  ba  developed.  And  third,  a 
method  for  calculating  tha  affact  of  ground  roughness  on  tha 
attenuation  factor  for  fallout  sourcas  had  to  ba  developed. 
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2.2  PoiDi-^SCDfl-OltbSdl 

In  Shi  aiding,  tha  traatmant  of  uncollidad  photon  -flux 
ia  ganarally  a  simple  exponential  attanuation  calculation. 

A  point  aourca  of  atrangth  (So)  (*/aec) ,  in  an  infinite 
homogeneous  medium  characterized  by  a  linear  attanuation 
coefficient  /*.  (1/cm),  will  have  an  uncollidad  flux  0 
(#/cma/sec>  at  a  distance  R  given  by 

(j>  -  Bo  exp(-/c  R)/(4  ^  Ra).  (2.1) 

The  factor  wxpi-f*,-  R)  ia  dua  to  material  attenuation  and  tha 
geometry  factor  1/(4  •n'  Ra)  is  due  to  the  inverse  aquare 
la*.  The  combined  factor 

Q(R>  -  exp(-/t  R)/(4 -rr  R=»)  (2.2) 

ia  referred  to  aa  the  point  kernel.  The  uncollidad  flux  for 
any  geometry  can  be  obtained,  at  leaat  in  principle,  by 
integrating  the  point  kernel  over  tha  geometry. 

When  tha  effect  of  collided  flux  ia  taken  into  account, 
the  problem  becomea  a  difficult  tranaport  problem.  To  get 
around  thia  difficulty,  a  semi empirical  factor,  the  buildup 
factor  (BUF  or  B)  ia  uaed  to  correct  for  the  contribution 
from  acattered  flux.  There  are  many  typea  of  buildup 
factora  including  the  "number  buildup  factor" ,  the  "energy 
buildup  factor",  the  "energy  abaorption  buildup  factor"  and 
the  "doae  buildup  factor".  In  this  dlacuaaion  only  the  doae 
buildup  factor  la  of  intereat.  The  doae  buildup  factor  ia 
defined  aa  the  ratio  of  the  total  doae  rate  at  a  point  to 


the  dose  rate  due  to  the  uncollided  flux  at  the  same  point. 
The  point  kernel  is  a  function  of  the  photon  energy,  shield 
thickness,  and  the  atomic  number  of  the  shield  material  (Z) 
C173.  The  point  kernel  for  the  flux  at  a  distance  R  becomes 
0<R>  -  B  exp<-/«.  R)/<4  re  R=> .  (2.3) 

Here  the  value  of  the  calculated  flux  is  weighted  to  account 
for  the  reduced  dose  rate  due  to  the  lower  linear  energy 
transfer  rate  of  the  lower  energy  collided  photons.  Because 
the  total  dose  rate  is  always  greater  than  the  dose  rate  due 
to  uncollided  flux  only,  buildup  factors  are  always  greater 
than  1. 

As  explained  in  Section  1.4,  even  though  the  term 
"dose"  is  used  throughout  this  dissertation,  the  programs 
developed  here  calculate  exposure.  This  should  not 
represent  any  real  problem  because  this  dissertation  deals 
with  the  ratio  of  exposures  which  is  equivalent  to  the  ratio 


of  doses  (DRF) 
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2.3  A£SyC*tf_ByilduB_E*£tors 

Finding  accurate  equations  for  the  required  buildup 
factors  can  be  difficult.  The  models  require  buildup 

factors  over  a  range  of  energy  from  15  KeV  to  15  MeV  and  a 

range  of  mean  free  paths  (mfps)  from  0  to  greater  than  40 
mpfs.  Nearly  all  the  works  reporting  buildup  factor  data 
were  published  before  1970.  As  lata  as  1960,  even  the  most 
comprehensive  works  only  reported  data  for  G.5  to  10.0  MeV 
and  out  to  20  mfps  C1B3.  Further,  even  when  data  are 
available,  the  most  accurate  formulas  commonly  used  can 
deviate  from  the  true  buildup  factors  by  more  than  60X  C193. 

The  greatest  deviations  occur  between  0.1  and  2.0  MeV, 

which,  unf artunately,  is  where  the  most  interest  lies.  The 
three  most  common  formulas  for  estimating  buildup  factors 
arc  Taylor's  formula  (a  three  parameter,  two  term 
exponential  equation),  Berger's  formula  (a  two  parameter 
exponential  equation) ,  and  a  polynomial  formula  which 
usually  has  three  terms,  but  is  often  used  with  two  terms 
and  as  such  is  good  only  for  thin  shields  C20,213. 


IfiXlg£.ll_fS£!DUl«i 


(2.4) 


B  ■  Aexp<-  *ytcr)  +  <1  -  A)axp(-  b  r) 

ifCfllCll-fSCCUlii  <2. 3) 

B  ■  1  +  a  -fju  r •  exp (b  ytc  r) 

EcUnamUL-iacnuU*  (2.6) 

3 

B  ■  ^  aj  r)1 

t«o 

IwS_£sC!B_H2lyQS!!!iSl_fSEfl!lil«  (2. 7) 

B  *  1  +  k  <  ju.  r) 

A  further  complication,  and  a  littla  raportad  -fact,  ia  that 
tha  paramatara  -for  tha  Juat  mantionad  formulae  may  hava  baan 
calculated  uaing  criteria  daaignad  to  prevent 
underaati mation  of  doaaa  C22].  In  other  worda,  data  takan 
by  anginaara  aa  accurately  rapraaanting  phyaical  reality 
hava  baan  altarad  to  prevent  them  from  aerioualy 
underaati mating  doaaa.  Aa  an  example,  Taylor 'a  formula  at 
0.5  MeV  haa  a  total  deviation  of  approximately  44%  over  0  to 
40  mfpai  paramatara  reported  in  tha  literature  C233  cauaa  a 
deviation  of  from  +41. 3X  to  -  5.4%  inatead  of  +22.%  to  -22. % 
Z221 .  Bacauaa  ona  of  tha  goal  a  of  thia  paper  ia  to  compare 
the  expoaure  received  for  different  apectra  aa  accurately  aa 
poaaible  to  allow  an  intelligent  eel action  of  alternativea, 
data  with  thia  kind  of  inaccuracy  cannot  be  uaed. 

A  much  more  accurate  formula,  along  with  the  required 
paramatara  for  water  haa  recently  baan  publiahed  by  A. 
Fodararo  and  R.  J.  Hall  C243.  Thia  ia  a  three  term,  five 


16 


parameter  formula  similar  to  Taylor's  formula.  Ths  average 
deviation  of  this  formula  is  reported  to  be  less  than  l.X  at 
most  energies  and  never  exceeds  4.X  (the  maximum  deviation 
is  at  0.1  MeV  and  40  mfps). 

Iticsszsa.ESQSQi.if  L_i2C2JuLi  <  2  •  a  > 

3 

B  ■  y  Al  exp  <-Q(.i •/»  r ) 
isl 

where  A3  ■  1  -  A,  -  A*. 

The  buildup  factors  for  air  and  water  as  a  function  of 
energy  and  mean  free  path  are  very  nearly  identical  C253. 
This  is  because  the  average  "Z"  number  jf  air  and  water  are 
nearly  the  same.  A  comparison  of  the  two  buildup  factors  is 
shown  in  Figure  2  at  20  mfps.  The  error  for  fewer  mfps  is 
much  less  than  at  20  mfps.  Fortunately  the  greatest 
deviation  between  the  air  and  water  buildup  factors  occurs 
at  large  mean  free  paths  where  attenuation  makes  the  error 
in  calculating  exposures  insignificant.  The  program  used  in 
this  paper  uses  the  three-exponential  formula  (equation  2.8) 
with  the  parameters  given  for  water  to  find  the  buildup 
factors  for  air  C243. 


10  KeV 


100  KeV 


1  MeV 

ENERGY 


10  MeV 


100  MeV 


FIGURE  2  e_CB!DG*!:i*aD_Bf_i£]*_lyilriUB_E«EiBC»_C«l£Ul*i!d 
fCB!D_ibS_IbCB#rEaBBD«Dti*l_EBCBMll_fBC_!!lli*C_Hiib_i;t3«_Iru* 

BulldBB_E«£ifiE9_fBr_dlrx 
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Finding  accurate  buildup  factors  for  a  gcnsral 
structurs's  Malls  and  roof  is  also  difficult.  Fortunately, 
construction  materials  are  almost  completely  constituted  of 
low-Z  elements  and  thus,  after  adjusting  for  density 
variations,  wood  and  concrete  are  very  nearly  equivalent 
£263.  Therefore,  the  mass  attenuation  coefficients  and 
buildup  factors  of  what  is  referred  to  as  HNBS"  Concrete 
were  chosen  to  represent  building  materials.  The  formula 
contains  six  non-linear  terms  and  eleven  parameters  £273. 

It  was  developed  by  the  National  Bureau  of  Standards  and  is 
the  standard  for  concretei  as  such  there  is  no  error 
associated  with  it  as  there  would  be  if  an  approximation 
formula  had  been  used.  Using  concrete  as  an  approximation 
fo-  building  material  in  general  is  much  more  accurate  than 
using  the  buildup  factors  for  water  (or  compressed  air)  as 
is  commonly  done  £283.  The  practice  of  using  the  buildup 
factors  for  water  seems  to  be  a  continuation  of  an  earlier 
practice  developed  to  avoid  a  limitation  of  the  moments 
method  where  all  materials  are  considered  to  be  water  or  air 
of  varying  density  £29,3213.  However,  the  buildup  factors 
for  water  can  be  more  than  122)  times  that  of  concrete  for 
large  mean  free  paths  and  energies  near  0. 1-MeV  £313. 

All  of  the  buildup  factors  used  in  the  program  are  "Air 
Kerma  Response  Function"  buildup  factors.  These  buildup 
factors  are  meant  to  be  used  to  calculate  exposure  and  dose. 
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2.4  £ombiQing_giJildye_F2£igr! 

A  major  problem  in  uaing  the  point  kernel  method 
occurs  whan  tha  photons  pass  through  mora  than  ona  madia.  A 
singla  buildup  -factor  must  ba  -found  that  raprasants  tha 
buildup  through  all  tha  layars  of  material  batwaan  tha 
source  and  tha  detector.  In  order  to  understand  this 
problem  it  is  necessary  to  consider  the  definition  of  a 
buildup  factor. 

Tha  buildup  factor  is  tha  ratio  between  tha  detector 
response  to  tha  total  radiation  at  a  point  of  interest  over 
tha  detector  response  to  the  uncollided  radiation  at  the 
same  point  C323.  Buildup  factors  are  necessarily  a  function 
of  four  vari ablest  the  composition  of  the  material  through 
which  the  radiation  passes  (usually  Just  referred  to  as  1 
number  dependence) ,  the  geometry  of  the  source  and  tha 
detector,  the  number  of  mean  free  paths  (mfps)  that  the 
radiation  passes  through  and,  the  energy  of  the  uncollided 
radiation.  By  using  mfps  instead  of  physical  distance,  the 
formulas  for  calculating  buildup  factors  are  independent  of 
material  density.  The  number  of  mean  free  paths  remain  a 
function  of  material,  density  and  distance.  Buildup  factors 
are  calculated  using  the  Monte-Carlo  or  moment  methods,  or 
found  from  empirical  data.  Usually  the  medium  is  assumed  to 
be  infinite. 
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The  aquation  for  tha  doaa  from  a  simple  point  aourca  in 
a  continuous  homogeneous  infinita  madium  isi 

Dr  ■  Da/Ra  B(Ea,mfpTfZ)  •  exp(-mfpr)  (2.9) 
Where*  mfp-r  -  tha  number  of  mean  free  paths 

between  tha  source  and  tha  detector) 

R  ■  tha  gaomatrical  distance  between  tha 
source  and  the  datactor  in  units  of 
distance) 

Dr  =  the  dose  at  the  detector) 

Do  »  the  dose  at  a  unit  distance  from  the 

source  in  the  absence  of  any  medium  in 
the  units  of  dose  times  distance 
squared)  and 

B<Eo,mpf f Z)  or  B  *  the  Buildup  factor  as  a 
function  of  source  energy,  mean  free 
path,  and  material  Z  number. 


By  defining  Do  at  the  same  distance  from  the  source  as  the 
detector,  (this  can  be  done  because  R  is  in  arbitrary 
distance  units)  R  becomes  equal  to  1  and  the  equation  is 
simplified  toi 

D«  ■  Do  •  B(Eo, mfp-r, Z)  •  exp(-mfpx)  (2.10) 


or  i 


B  ■  B(Eo, mfp-r, Z)  9  Dr/Do  -  exp(mfpx)  (2.11) 
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If  we  taka  a  continuous  medium  and  brsak  it  into  two  rsgions 
(sae  Figurs  3) i 

B  *  Da/Di  •  Dx/Do  ■  exp(mfp»)  •  exp(mfp=). 

(2. 12) 

Buti  B(£o,mfpi  , Z)  ■  tha  buildup  -factor  -from  point 

"0"  to  point  "1" 

■  Di/Do  •  sxp(m-fpi).  (2.13) 

Thus  ws  can  think  o-f  ths  buildup  -factor  from  point  "1"  to 
point  "2"  as 

B(Eo,mf p3,Z)  -  Da/D,  ■  s?cp(,7.fp3)  (2.14) 

andt  B  ■  B(Eo,mfp,,Z)  -  B  (Eo,mf p3 , Z) 

5  B(Eo,mfpT,Z) .  (2. 13) 

Clearly  buildup  factors  muat  ba  multipliad  in  ordar  to  ba 
combined.  Tha  difficulty  lias  in  finding  B(Eo,mfp3,Z) .  For 
a  continuous  madiumi 

F(Eo,mf p3,Z)  S  B  (Eo,mf  p-r ,  Z)  /B  <Eo ,mf  p i  , Z) . 

(2. 16) 

B(Eo,mfp3,Z>  raprasants  tha  additional  buildup  from  point  1 
to  point  2.  It  takas  into  account  tha  buildup  and  shift  in 
anargy  spactrum  from  point  0  to  point  1  (saa  tha  insart  in 
Figura  3).  Tha  problam  bacomas  apparent  whan  ona  raalizas 
that  buildup  factors  arm  only  tabulated  for  continuous 
madia.  If  tha  madia  in  region  "Zi”  and  region  "Z3"  are 
different,  wa  can  find  B(E=>,mfp,  ,Z,)  from  tables,  but  not 
B(Eo,mfpv,Zi ,3)  or  B <Eo,mf p3 , Z3> .  Hare  Z,.3  represents  the 
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composite  region  made  up  raglona  Z>  and  Za.  In  fact  tha 
valua  of  B (Eo,mfpr, Z» .a)  is  depandant  on  tha  geometry  of  the 
problam  as  wall  as  tha  data! lad  composition  of  both  ragions. 
Thus  Ma  cannot  exactly  datarmina  B (Eo,mf pr,Z» .a)  by  tha 
point  karnal  mathod  uni ass  tha  sama  spacific  problam  is 
first  solvad  aithar  empirically  or  by  tha  Monta-Carlo 
mathod. 

Several  methods  have  baen  developed  over  tha  years  for 
estimating  B(EotmfpT,Zi.a> .  Tha  most  common  mathod  simply 
multiplies  tha  buildup  factor  B(Eo,mfp* , Z»>  by  B(Eo,mfp2,Za> 
[32,33,34,33,363.  This  method  ovarastimatas  tha  composite 
buildup  factor  B(Eo,mfpT,Zi .a)  [32,353,  and  therefore  gives 
an  upper  limit. 

Another  mathod  developed  by  D.  L.  Brodar  C373  for  a 
narrow  range  of  energies  is  to  combine  tha  buildup  factors 
as  followsi 

B  ( Eo , ftif  px ,  Z i ,  a)  ■  B  ( Eo , mf  p i  ,  Z  i ) 

+  B(Eo,mfp-r,Za)  "  B (Eo,mf p * , Za)  (2.17) 

Tha  mathod  can  be  extended  to  more  than  one  region.  Brodar 
did  not  discuss  the  development  of  Equation  2.17  except  to 
say  that  "it  can  be  recommended  for  calculating  buildup 
factors  for  heterogeneous  media  at  energies  near  1  Mev", 
and,  "...it  may  be  assumed  that  the  derived  equation  is  also 
applicable  for  larga  energies".  It  doas  fit  his  data,  but 
ovarastimatas  tha  buildup  factors  at  lower  energies  (see 
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Figure  3).  Broder's  method  and  a  few  slight  variations  of 
it  are  vary  popular  where  more  accurate  combined  buildup 
factors  are  sought  but  the  restrictions  on  whan  they  are  to 
be  used  are  seldom  mentioned  [32,33,34,353 . 

Another  commonly  used  method  takes  only  the  larger  of 
B(Eo,mfp, ,Zi)  or  B (Ea,mf p3 , Za) .  This  underestimates  the 
buildup  factor,  but  it  does  give  a  lower  bound. 

The  method  used  here  for  estimating  the  combined 
buildup  factor  is  derived  from  Equation  2.16.  The  buildup 
factor  for  region  "Za"  is  estimated  by  taking  the  ratio  of 
the  buildup  factor  found  by  using  the  material  of  region 
"Za"  and  the  total  number  of  mean  free  paths  from  the  source 
to  the  end  of  region  "Za"  divided  by  the  the  buildup  factor 
found  by  using  the  material  of  region  "Za"  and  the  total 
number  of  mean  free  paths  from  the  source  to  the  beginning 
of  region  "Za"* 

B ( Eo , mf pa , Za)  ■  B ( Eo , mf px , Za) /B ( Eo , mf p i , Za) . 

<2. IB) 

This  method  takes  into  account  the  development  of  the 
spectrum  as  it  penetrates  region  "Zi"  by  assuming  the 
spectrum  would  be  similar  to  that  developed  if  it  penetracad 
the  same  number  of  mean  free  paths  of  material  "Za".  If  the 
media  in  regions  "Zi"  and  "Za"  are  the  same  (or  even  if  they 
are  made  of  the  same  elements  but  with  different  densities) 
Equation  2.18  reduces  to  the  identity  Equation  2.16.  In 
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fact  this  assumption  should  bs  vary  good  for  raglons  of 
similar  Z  numbar,  bacausa  as  shown  by  plots  of  tha 
diffarantial  anargy  spactrum  for  watar  and  aluminum  dua  to 
mono-anargatic  gamma  rays,  thara  is  littla  diffaranca  in  tha 
shapa  of  tha  spactrum  produced  in  materials  of  Z  number  lass 
than  13  [303.  Tha  "Z"  number  of  concrete  is  usually  taken 
as  13  also. 

Determining  tha  exact  accuracy  of  this  method  would 
require  empirical  or  Monte-Carlo  solution.  However,  this 
method  passes  two  important  tests  that  are  required  of  an 
accurate  method.  First,  its  results  generally  lie  between 
the  known  upper  and  lower  bounds  over  the  entire  anargy 
range.  The  only  exception  is  at  approximately  150  KeV,  where 
the  shift  in  the  photon  energy  spectrum  that  occurred  in  tha 
first  medium  altered  tha  BUF  curve  of  the  second  medium  as 
expected.  Second,  for  the  available  data  (at  1.25  MaV) 
agreement  is  vary  good  as  shown  in  Figure  4  where  Border's 
data  and  his  formula  also  agree.  The  method  developed  here 
is  the  only  one  that  passes  both  tests. 


Iff  KeV 


100  KeV 


1  XeV 
ENERGY 


10  MeV 


100  MeV 


figure  4  EBfflfisriseo_ef_Bssyits_fif_B*l;bBffi_ef_GEfi!feiQiQa 
§yildyB_E*siflcs_fBC_i_iwB_caai9D_BCBbla(Dx 
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2.5  CgmB*ri3Qn_fif_§2yc£S_SB*EiC« 

While  it  is  not  a  purpose  of  this  dissertation  to 
speculate  on  the  protection  offered  by  a  shelter  for  any 
particular  release  event,  it  is  necessary  to  calculate  the 
protection  provided  for  a  few  well  documented  spectra  in 
order  to  compare  results  and  draw  conclusions  relevant  to 
the  current  literature.  How  the  spectra  for  comparison  are 
obtained  is  of  some  interest.  The  programs  used  here 
require  the  source  to  be  reported  by  each  gamma  energy  and 
the  fraction  of  disintegrations  resulting  in  that  energy  for 
each  curie  emitted  by  the  source. 

For  releases  reported  as  curies  of  individual  isotopes, 
the  Job  is  easy.  The  only  assumption  that  needs  to  be  made 
is  that  the  isotopes  are  uniformly  mixed.  Even  this 
assumption  may  be  relaxed  when  the  models  used  here  are 
coupled  with  a  dispersion  model.  Care  must  be  taken  that 
all  daughter  isotopes  are  found.  Table  1  gives  an  example 
of  this  kind  of  release  data.  It  contains  the  release 
reported  for  the  first  33  hours  of  the  Three  Mile  Island 
Accident  in  March  1979  C393.  In  Table  2  the  energy  peaks 
for  each  of  the  isotopes  and  their  daughters  have  been 
catalogued  C403.  Table  3  shows  the  data  reduced  to  a 
suitable  spectra  by  multiplying  the  fraction  of  the  isotope 
in  the  release  by  the  fraction  of  disintegrations  that 
result  in  that  particular  gamma  photon  energy.  Similar 
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energies  have  baan  summed.  Nota  that  tha  iodine-131  plays 
an  insignificant  rola  in  tha  cloud  sourcai  however ,  it 
comprises  the  antira  ground  sourcn  (Table  4).  In  a  similar 
mannar  tha  fallout  spectrum  for  tha  1961  SL-1  accident  can 
be  found  from  tha  10  Curias  of  1-131,  0.3  Curia  of  Cs-131 
and  0.1  Curia  of  Sr-90  ralaasad  C413.  Nota  that  Sr-90  plays 
no  rola  in  the  shielding  calculation  sinca  it  is  a  pure  beta 
emitter.  Bramsstrahlung  is  not  included  in  tha  spectrum 
calculations  dua  to  tha  low  average  energies  generated,  high 
shielding  factors  at  thasa  energies,  and  tha  low 
bramsstrahlung  yiald  in  tha  low  Z  number  materials  involved. 

When  ralaasaa  are  raportad  as  spactra  in  tha 
litaratura,  they  are  often  raportad  in  energy  groups  and  as 
tha  relative  energy  content  of  tha  gamma  rays  emitted  in 
earh  group.  Thasa  spectra  must  ba  converted  to  energy 
groups  and  tha  relative  number  of  gamma  photons  emitted  in 
each  group.  Tables  6  and  7  C423  give  examples  of  this  type 
of  conversion.  In  these  tables,  Mf"  is  the  fraction  of  the 
total  photon  energy  contained  in  each  energy  group.  To 
convert  to  the  relative  number  of  photons  in  each  group,  "f” 
is  divided  by  the  average  energy  of  each  group  and 
renormalized  to  1.0.  The  resulting  fraction  is  used  by  the 
programs  to  calculate  the  dose  reduction  factors,  but  cannot 
be  used  to  find  the  dose  per  Curie  of  release,  unless  the 
isotopic  composition  of  the  radioactive  material  is  known. 
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This  is  because  ths  number  of  curias  required  to  produce  a 
given  number  of  photons  is  unknown. 

A  much  easier  way  of  obtaining  spectra  is  to  find  them 
given  in  the  way  needed  -  by  energy  group  and  relative 
number  of  photons  in  each  group  as  was  found  for  the  weapon 
fallout  spectra  given  in  Table  8  C433.  Again  the  dose  per 
Curia  released  cannot  be  directly  found ,  but  the  dose 
reduction  factors  can  be  calculated.  When  these  spectra  are 
used,  the  exposure  for  each  average  group  energy  is  found 
and  multiplied  by  the  fraction  of  photons  in  that  group. 

The  results  are  summed  to  find  the  total  exposure  and 
divided  by  the  similarly  found  unprotected  exposure  to 
obtain  the  dose  reduction  factor. 

Figures  5  and  6  depict  the  relative  energy  of  each  of 
the  six  spectra  discussed.  In  order  to  provide  a  visual 
comparison  of  the  spectra,  all  six  spectra  were  converted  to 


energy  groups 
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TABLE  1 

ISOTOPIC  RELEASE  FOR  THE  FIRST  THIRTY-THREE  HOURS 
OF  THE  THREE  MILE  ISLAND  ACCIDENT  OF  MARCH  1979  C39D 


Isotope 

Curies 

Fraction 

Xe-133 

4.9  x  10**6 

0,729 

Xa-133m 

1.2  x  10**3 

0.018 

Xa-133 

1.3  x  10**6 

0.223 

Xa-135m 

1.4  x  13**5 

0.021 

K.—aa 

6. 1  x  10**4 

0.0091 

1-131 

1 . 908 

0.0000003 

total 

6i.Z2_S._i3**6 

1.9 

TABLE  2 


ENERGY  PEAKS  FOR  ISOTOPES  IN  TABLE  1 


I sot op* 

Enargy(MaV)  Fraction 

of  Diaintagrations 

Xa-133 

a.  081 

0.37 

Xa-133m 

a.  233 

0. 14 

0.081  (Xa-133) 

0.37 

X*— 133 

0.250 

0.91 

0.61 

0.03 

X«*-135m 

0.527 

0.80 

0.250  (Xa-133) 

0.91 

0.61  (Xa-135) 

0.03 

Kr — 88 

0.028 

0.07 

0.  166 

0.07 

0.191 

0.33 

0.36 

0.05 

0.85 

0.23 

1.55 

0.14 

2.19 

0.  18 

2.40 

0.35 

0.898  (Rb-88) 

0.13 

0.186  (Rb-BB) 

0.21 

2.68  (Rb-88) 

0.023 

1-131 

0.80 

0.026 

0.284 

0.054 

0.364 

0.82 

0.637 

0.06B 

0.723 

0.016 

0.164  (Xa-131m) 

0.0002 
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TABLE  3 

ENERGY  SPECTRA  DERIVED  FROM  TABLES  1  AND  2 


ENERGY  <K»V)  Fraction  of  gamma*  par  Curia  of  ralaaao 


2630. 

0.C00209 

2400. 

3. 20319 

2190. 

0.00164 

1863. 

0.00191 

13S0. 

0.00127 

898. 

0.00118 

850. 

0.00209 

610. 

0.00732 

527. 

0.0)68 

360. 

3.000455 

250. 

0.2220 

233. 

0.00252 

191. 

0.00319 

166. 

0.00637 

81. 

0.2764 

28. 

0. 000637 
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Tabla  4 

1-131  SPECTRUM  (TMI  PARTICULATES) 


Energy (KaV) 

Fraction  of  Gamma* 

723. 

0.816 

637. 

0.868 

364. 

0.82 

284. 

0.834 

164. 

0.0882 

88. 

8.826 

TABLE  3 

SL-1  FALLOUT  SPECTRUM  C41D 


Enargy(KaV)  Sourca  Fraction  of  Gamma*  par  Curia 


723. 

1-131 

0.0131 

662. 

C*-137 

0.8401 

637. 

1-131 

8.0642 

364. 

1-131 

0.774 

284. 

1-131 

0.8509 

164. 

1-131 

8.88019 

80. 

1-131 

0.0245 

1 
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TABLE  6 

RELATIVE  SOURCE  SPECTRA  FOR  A  CLOUD  SOURCE  FROM  A 
REACTOR  SAFETY  STUDY  PWR-2  EVENT  C42D 


E(MaV) 

Eava 

i 

i /Eava 

Fracti on 

0.01-0. 1 

0.05 

0.033 

0.  660 

0.3693 

0. 1-0.5 

0.30 

0.  164 

0.547 

0.3059 

0. 5-1.0 

0.75 

0.  197 

0.264 

0. 1476 

1. 0-2.0 

1.50 

0.279 

0.  186 

0. 1841 

2. 0-3.0 

2.50 

0.327 

0.,  131 

0.0732 

total 

laZiZ 

1.0 

i 

TABLE  7 

RELATIVE  SOURCE  SPECTRA  FOR  A  FALLDUT  SOURCE  FROM  A 
S  REACTOR  SAFETY  STUDY  PWR-2  EVENT  C423 


E(M«V) 

Eava 

i 

f /Eava 

Fraction 

0.01-0. 1 

0.05 

0.010 

0.200 

0. 1696 

0. 1-0.5 

0.30 

0.086 

0.287 

0.2432 

0. 5-1.0 

0.75 

0.244 

0.325 

0.2759 

1. 0-2.0 

1.50 

0.386 

0.257 

0.2183 

2. 0-3.0 

2.50 

0.274 

0.110 

0.0930 

total 

uaaa 

ii.iZ2 

li.9 

« 


1 


« 
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TABLE  B 

FINN-SIMMONS  1-HOUR  WEAPON  FALLOUT  SPECTRA  C43D 


E(M«V) 

Eava 

Photon  Fraction 

0.0-0.05 

0.025 

0.0271 

0.05-0. 10 

0.073 

0.0137 

0. 10-0.20 

0. 130 

0.0737 

0.20-0.30 

0.250 

0.0476 

0.30-0.40 

0.350 

0.0929 

0.40-0.60 

0.300 

0. 1373 

0.60-0.80 

0.700 

0. 1717 

0.80-1.00 

0.900 

0.1627 

1.00-1.33 

1.165 

0.0889 

1.33-1.66 

1.500 

0.0937 

1.66-2.00 

1.830 

0.0299 

2.00-2.30 

2.250 

0.0397 

2.30-3.00 

2.750 

0.0148 

3.00-4.00 

3.50L 

0.0042 

4.00-5.00 

4.300 

0.0001 

total 

ubbbb 

figure  s  CsoeftEiaBO-sf-Clcyd-Ssucca 
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2.6  grsund_s£f3S£3 

Thai  hypothetical  in-finite  smooth  plane  source  assumed 
for  point  kernel  integration  over  a  disk  source  does  not 
exist  in  nature.  The  fallout  particles  distribute 
themselves  over  the  real  terrain  surface.  Small  surface 
irregularities,  such  aa  grass,  gravel  and  concrete  are 
usually  called  ground  roughness  in  the  shielding  literature. 
For  an  infinite  plane  source,  most  of  the  gamma  rays 
originate  from  large  distances  and  travel  through  small 
grazing  angles  to  the  surface!  irregular! ties  in  the  surface 
cause  those  rays  to  be  attenuated  as  they  begin  their 
journey.  Larger  surface  irregularities,  such  as  hills, 
washes  and  buildings  are  called  terrain  effects.  These 
large  irregularl ties  nearly  always  reduce  dose  rates  (by  as 
much  as  30%  "for  a  person  standing  on  top  of  a  small  steep 
hill  that  falls  away  in  all  directions. ...  because  the  hill 
hides  much  of  the  fallout  beyond  the  immediate  area.")  E443. 
However  these  large  effects  are  not  general  in  nature,  being 
peculiar  to  specific  buildings  and  situations,  and  are 
therefore  not  considered  in  this  model.  This  dissertation 
is  concerned  with  modeling  effects  that  are  more  general  in 
nature  i.e.  the  shielding  provided  by  lawns,  streets,  etc. 
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2.6.1  Ecayiguajiggfis 

In  196S ,  the  Damans*  Civil  Preparedness  Agency 
[45 3  publiahad  reduction  -factor*  to  be  applied  to  varloua 
ground  roughness  conditions  for  a  fallout  source.  These 
data  are  repeated  in  Table  9.  These  data  Indicate  that 
ground  roughness  can  reduce  the  dose  rate  by  as  much  as  5 0X 
compared  with  1  meter  above  the  standard  hypothetical  smooth 
plane  source*  assuming  that  the  same  amount  of  fallout  is 
uniformly  deposited  in  both  cases. 

The  simplist  model*  and  the  one  most  commonly  found  in 
the  literature*  applies  the  above  factors  directly  to 
shielding  factors  found  for  buildings*  without  regard  to  the 
fallout  source  [44,463.  However*  because  the  reduction 
factor  is  caused  by  a  portion  of  the  gamma  photons  passing 
through  surface  irregularities  as  they  start  their  journeys* 
the  reduction  factors  must  be  source  energy  dependent. 

Another  possible  model  was  suggested  by  C.  M. 

Huddleston  in  1964  whan  he  observed  that  the  ground 
roughness  effect  of  various  Nevada  terrains  on  fallout  from 
an  atmospheric  bomb  test  was  equivalent  to  raising  his 
detector  20  to  40  ft  [473.  While  lifting  the  detector*  and 
thus  increasing  the  attenuation  caused  by  air  and  distance 
will  give  a  source  energy  dependent  model*  it  is  not 
satisfactory  for  the  following  reasons.  First.*  that  model 
exposes  all  the  gamma  rays  to  the  additional  material,  while 
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in  tha  raal  situation  only  that  -fraction  of  tha  photona 
behind  a  surfaca  irragularity  would  ba  affactad  by  it.  This 
would  offact  tha  rasultant  spectrum  at  tha  dotactor. 

Second,  that  modal  causa*  a  ralativaly  graatar  attanuation 
of  photons  originating  naar  tha  datactor  than  thosa 
originating  at  a  larga  distanca,  yat  tha  gamma  rays 
originating  from  long  diatancas  travel  through  small  grazing 
anglas  to  tha  surfaca  and  thus  should  ba  mo r?  graatly 
attanuatad  than  photons  originating  naar  tha  datactor  which 
have  a  naarly  unobstructad  line  of  sight. 


TABLE  9 

GROUND  EFFECT  DOSE  REDUCTION  FACTOR 
FOR  1.12  HR.  WEAPON  FALLOUT  C433 


Ground  Roughneaa  Condition  Reduction  Factor 


Smooth  plana  (hypothetical)  1.00 
Paved  areaa  1.00  to  0.83 

Lawna  0.83  to  0.73 

Gravelled  areaa  0.73  to  0.63 

Ordinary  plowed  field  0.63  to  0.55 

Deeply  plowed  field  0.55  to  0.47 
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2.6.2  £E22893d-tfSdal 

Thar*  is  vary  1  it'd  a  information  upon  which  to 
baae  a  ground  roughness  modal.  In  fact,  th»  reference  given 
in  Tabla  9  givas  only  ona  datum  point  -for  ona  particular 
■pactrum  (tha  1.12  Hr.  -fallout  spectrum)  for  aach  of  fiva 
surf  tea  typas.  In  ordar  to  ma^e  tha  bast  usa  of  aval  labia 
data  and  rasaarcht  a  modal  must  ba  vary  ganaral  in  nature  - 
but  must  hava  sufficient  data  available  to  define  all  its 
parameters.  Obviously  tha  ona  datum  point  available  will 
only  allow  tha  fitting  of  ona  parameter.  Any  additional 
parameters  must  ba  defined  by  other  mean:*. 

Tha  modal  used  in  this  paper  was  developed  from  ona 
used  to  comoare  tha  angular  distribution  of  dose.  rata  ovar  a 
plowed  dry  lake  bed  to  that  over  a  smooth  dry  lake  bed  after 
fallout  from  tha  explosion  of  a  nuclear  device  C4B1.  The 
model  reflects  the  geometry  of  a  plowed  field,  see  Figure  7, 
but  the  three  parameters  used  by  tha  model  make  it  quite 
general  in  nature.  The  parameters  aret  V'  «  the  trough 
argle;  w,  the  trough  width;  and  d,  the  width  of  the 
remaining  flat  surface.  The  major  difference  between  the 
model  used  in  this  paper  and  the  model  it  is  derived  from  is 
that  the  existing  model  assumed  the  concentration  of  fallout 
was  equal  on  all  surfaces  regardless  of  tilt,  while  the 
model  used  here  assumes  that  the  fallout  is  uniform  on  the 
horizontal  projection  of  the  surface.  The  use  of  the  two 
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modsls  is  also  diffsrant.  Whils  ths  prsvious  modal  mss  ussd 
to  prsdict  spscific  doss  valuas,  ths  modsl  ussd  in  this 
papsr  is  ussd  to  prsdict  doss  valuss  at  diffsrsnt  snsrgiss. 
Instsad  o-f  -finding  a  spscific  doss  rsduction  valus  from 
known  paramatars,  ths  modsl  must  usa  known  doss  rsduction 
valuss  to  find  ths  unknown  paramatars. 

Ths  problsm  is  to  dotsrmina  ths  valuss  of  thraa 
paramatars  whan  only  ons  datum  point  (Tabla  9>  is  givsn  for 
sach  surfacs  condition.  Bacauso  ths  modsl  is  ussd  as  it  is, 
this  problsm  can  ba  solvsd. 


v  ■ 


figure  7  J3acmficy_af.iba_ecauDil_Effacia_t!ori9l 
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The  concept  of  a  "half -shadow"  angle  wee  introduced  by 
Yu.  A.  Izrael  in  1963  C493.  The  ratio  d/w,  for  a  given 
value  of  y  ,  can  be  found,  at  least  in  many  cases,  from  the 
half -shadow  angle  •  The  half -shadow  angle  is  defined 

asi  the  angle  above  the  horizontal  from  which  the  ground  is 
illuminated  causing  half  the  surface  to  be  in  shadow  when 
viewed  from  above.  As  an  examplei  "when  meadow  land  is 
illuminated  from  an  angle  of  1-degree  above  the  horizontal 
and  photographed  from  above,  about  50  percent  of  the  surface 
is  in  shadow  C581.  The  half -shadow  angle  can  be  thought 
of  as  being  measured  from  the  same  reference  as  p  in  Figure 
7.  The  relationship  between  d/w,  ^  and  V*7  is  as  followsi 

d/w  «l.-2.  cosC?')  sin  (££)/»! nC^+^jt)  (2.19) 

for  y. 

As  can  be  seen  a  /3f  does  not  exist  if  d  is  greater  than  w 
— as  in  the  case  of  the  previously  referenced  plowed  dry 
lake  bed  where  d  was  18  inches  and  w  was  12  inches  C4B1. 
However  for  the  more  usual  cases  /?£  does  exist.  In  fact  it 
has  been  demonstrated  that  is  proportional  to  the  dose 
reduction  factor  C493.  Yu.  A.  Izrael 's  data  give  the 
following  equation  for  /3>  where  7  is  the  dose  reduction 
f actori 

pi  -  16.43  -  19.85  7  degrees.  (2.28) 

Izrael 's  data  for  /?£  and  7  for  various  ground  conditions  are 
given  in  Table  18. 


TABLE  10 


HALF-SHADOW  ANGLE  AND  DOSE  REDUCTION  FACTOR 
FOR  VARI0U3  GROUND  TYPE3  C493 


Ground  Typ® 

/£?>  ,dagra®a 

7  (DRF 

Vary  -flat,  virgin  grasay 
region  (meadow,  clearing) 

1.0 

0.31 

Vary  flat  region  of  arid  atapa 

1.7 

a.  77 

Arid  ataps 

3. a 

0.69 

Cultivated  fiald 

3.3 

0.60 
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Becauaa  tha  modal  ia  usad  as  it  is*  to  -find  tha 
paramatar  w  that  givaa  a  givan  doaa  raduction  -for  a  givan 
apactrum,  tha  valua  choaan  for  y  makws  vary  littla 
diffaranca  in  tha  raaulta  for  varying  anargiaa.  Tha  doaa 
raduction  factora  ovar  a  plana  with  no  atructura  praaant 
were  comparad  ovar  a  Mi  da  ranga  of  Y  'a  US  to  43  dagraaa) 
and  half-ahadow  anglaa  (1  to  5  dagraaa).  Tha  raaulta  ahow 
that  tha  diffarancaa  ara  insignificant  for  varying  y  'a  for 
anargiaa  abova  80  KaV  (oaa  Figuraa  8  and  9).  Whan  a  small 
cantral  void  was  inaartad  in  tha  plana*  rapraaanting  tha 
ramoval  of  tha  ground  sourca  under  a  atructura  (tha  fallout 
falling  on  a  atructura  ia  accounted  for  aa  a  roof  source)* 
tha  diffarancaa  ara  even  smaller  (Figure  10).  .Ind  whan  a 
wall  ia  added  tha  diffarancaa  become  totally  insignificant 
far  all  anargiaa  aa  shown  in  Table  11.  What  happans  ia  that 
a  specific  doaa  raduction  factor  is  specified  for  a  specific 
gamma  spectrum  (in  this  case  0.81  for  Co-60  C49D*  y  and  /?£ 
or  d/w  ara  chosen  (  Y  and  determine  d/w)  and  tha  "Ground 
Factor"  program  (see  Appendix  A)  iteratively  determines  a 
valua  of  w  such  that  tha  DRF  equals  tha  givan  valua.  Tha 
values  of  w*  d/w  and  Y  ara  usad  in  tha  "Fallout"  program 
(Appendix  B)  to  determine  tha  DRF  for  other  conditions  where 
structures  ara  prasent.  For  different  values  of  Y  and  a 
givan  /?£ *  tha  program  simply  finds  a  diffarant  valua  of  w  — 
which  results  in  vary  nearly  tha  sama  DRF  vs  Energy  curves. 
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y  The  only  significant  diffarancas  (for  flat  planaat  saa 

Ju 

Figuraa  8  and  9)  ara  at  low  anargias  whara  tha  attenuation 
of  building  walla  is  ao  graat  that  tha  variance  in  the  DRF 
~j  causad  by  amall  changaa  in  ground  af  facta  ia  no  1  on  gar 

significant,  as  ;hcwn  in  Tabla  11. 

] 
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J 
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10  KeV 


100  KeV 


1  MeV 

ENERGY 


10  MeV 


100  MeV 


figure  8  Bai*_BiductiaQ_EiciBC*_fcc_*_biiifrSb*riD»_flDalf_Bf 
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1  MeV 
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figure  9  Bea*_B*duEtiDD_EaEiQcs-fBc_*_B*IfrSbadsH_0Qfll3_Gf 

S-Baacaai-Svsr  *Q_lQ£iQi£i  Plana 


10  KeV 


100  KeV 


1  MeV 

ENERGY 


10  MeV 


100  MeV 


FIGURE  IB  CMB_B!dUCUDO-E*£&BCl-fS£-i-tiAlf::SbldS*-&Qglf-Bf 
l-tai£M..oyrc.AQ_lQUal£i-eiioi-sUb-*-£iD£cil_¥Bid-uUb-i-3 
M&rc-Bftdiui 
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TABLE  11 

EFFECT  OF  VARYING  THE  TROUGH  ANGLE  ON  THE 
DCSE  REDUCTION  FACTCR  FOR  A  SMALL  HOUSE 


5 

V 

V. 

I  ,■ 

v  ' 

] 


i 


K 

r 

J 


i 

£■ 


f 

I 


\  i 


> 

f. 


v--. 

V-; 


□N  A  LAWN 


Enargy 

(KaV) 

DRF  inr 
y  ■  45  dag. 

DRF  for 

y  -  15  dag. 

p-jrcant 
di  t-f  aranca 

12300. 

0.329 

3.326 

3.234 

13000. 

0.533 

8.230 

8.574 

8000. 

0.533 

3.330 

0.361 

6080. 

0,329 

0.326 

0.535 

2300. 

8.524 

3.522 

8.522 

4000. 

0.521 

0.519 

0.319 

3000. 

0.508 

3.536 

0.506 

2000. 

0.487 

0.486 

0.486 

1230. 

0.471 

0.471 

8.471 

1300. 

0.458 

0.438 

0.037 

800. 

0.452 

0.452 

-0.044 

600. 

0.444 

0.445 

-0. 109 

200. 

0.442 

0.443 

-0.281 

430. 

0.439 

0.440 

-0.267 

300. 

0.  433 

0.434 

-8.347 

230. 

0.379 

3.400 

-0.643 

153. 

0.359 

8.362 

-0.366 

130. 

0.244 

0.248 

-1.783 

88. 

8.  137 

0.  162 

-2.738 

60. 

8.0592 

0.0617 

-4.247 

20. 

0.0238 

8.0220 

-3.390 

40. 

0.00215 

0.00230 

-7,810 

30. 

0.342E-2 

0.384E-5 

-12.569 

23. 

0.  U7E-16 

0. 117E-16 

0*0 

12. 

3.3 

0.0 

0.0 

The  Fallout  program  needs  a  representati ve  value  for 
tha  numbar  of  mean  fraa  paths  that  tha  ground-attenuated 
gamma  rays  pass  through.  It  also  raqulras  tha  fraction  of 
gammas  that  pass  through  tha  ground.  Tha  fraction  of  gammas 
that  pass  through  tha  ground  1st 

(a  +  w)/(d  +  w) .  (2.21) 

where: 

a  -  w  cos  if)  ■  sin  (/3)  /sin  .  (2.22) 

The  problem  is  in  finding  a  raprasantati va  value  of  t  (sea 
Figure  7).  Choosing  the  average  value  of  t  works  very  well. 
The  maximum  value  of  t  isi 

tnAx  «  w  sin(^>/(2  cos(^)  sin(yS)).  (2.23) 

The  average  value  of  t  1st 

t«v,c  »  ti-,«x /2.  (2.24) 

In  order  to  find  a  raprasantativa  value  of  t  (tnc») ,  tha 
ground  buildup  and  attenuation  must  be  considered. 

Therafora  the  mean  value  of  tha  kernel  was  found  over  tha 
integral  from  x  to  Umax  (Figure  7)  and  a  value  of  t  was 
iteratively  found  such  that  its  kernel  and  the  mean  Just 
found  were  equal. 
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(2.25) 


Where 


X  MAX  **  Nsin^-/3)/(2’COS(n) 


(2.26) 


and 


t(x)  -  cot  ()»'-£>  +  cot  (^ ) 

-  x  »in(y')/(*in(/9)  »in(/-^))  (2.27) 

Tha  results  ara  given  in  Table  12  and  are  compared  to 
the  results  obtained  by  using  tAVE  as  the  raprasentati vu 
value  o-f  t.  Clearly  there  is  no  significant  error  in  using 
tAWE<  However  there  is  a  significant  savings  in  computer 
time.  Again  the  program  simply  finds  a  value  of  w  such  that 
the  givan  DRF  is  obtained.  Using  t«Ef*  instead  of  tAVE 
simply  changes  the  value  of  w  —  but  not  the  resulting  DRFs. 


* 
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TABLE  12 

EFFECT  GF  USING  T^«  VS  T„„-  IN  DETERMINING 
THE  DOSE  REDUCTION  FACTOR  FOR  VARYING  ENERGIES 
OVER  A  MEADOW  LAND 


Enargy 

(KaV) 

DRF  for 

tncp 

DRF  for 

t*VE 

Parcant 
dif f aranca 

13000. 

0.S32 

0.828 

0.41 

10000. 

0.833 

0.831 

0.33 

8000. 

0.834 

0.832 

0.31 

6000. 

0.834 

0.831 

0.26 

3000. 

0.833 

0.831 

0.24 

4000. 

0.832 

0.830 

0.19 

3000. 

0.826 

0.823 

0.13 

2000. 

0.817 

0.816 

0.087 

1300. 

0.809 

0.808 

0.032 

1000. 

0.811 

0.811 

-0.005 

800. 

0.810 

0.811 

-0.049 

600. 

0.807 

0.808 

-0.13 

500. 

0.811 

0.812 

-0.17 

400. 

0.813 

0.815 

-0.23 

300. 

0.827 

0.820 

-0.40 

200. 

0.797 

0.802 

-0.62 

130. 

0.793 

0.799 

-0.78 

100. 

0.745 

0.748 

-0.48 

80. 

0.703 

0.706 

-0.091 

60. 

0.665 

0.662 

0.44 

30. 

0.649 

0.642 

1.03 

40. 

0.644 

0.633 

1.75 

30. 

0.647 

0.631 

2.48 

20. 

0.675 

0.662 

1.91 
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2.  6. 3  32aulia_2f_3CSyDd_C2UafcD333_3Sd2l 

The  buildup  f actcri  and  mams  attenuation 
coaficients  for  concrata  are  uaad  in  tha  modal  to  represent 
tha  ground  surface.  Tha  "ground  -factor"  is  determined  by 
iterative  maans  using  tha  Finn-Simmons  13-group,  1-hour  fast 
-fission  dalayad  gamma  ray  spactrum  which  is  raprasantati v« 
of  a  typical  bomb  fallout  spactrum  C311. 

Tabla  13  givas  tha  nacassary  valuas  of  d/w  and  w  to 
obtain  tha  DRFa  givan  in  Tabla  9.  Tha  trough  angla  is 
assumad  to  ba  43  dagraas. 

Of  coursa,  tha  ground  roughnass  shielding  factor  cannot 
ba  usad  diractly  on  tha  disk  sourca  karnal .  First  tha  total 
buildup  factors  for  both  tha  diract  and  indiroct  photons 
must  ba  datarminad  for  tha  ground,  outaida  air,  building 
walls,  and  insida  air-  as  demonstrated  in  tha  saction  of  this 
papar  on  combining  buildup  factors.  In  thia  section,  t h  i 
building  itsalf  was  laft  out  (except  whara  tha 
insignificance  of  varying  trough  angles  waa  discussed);  only 
attenuation  by  tha  ground  and  air  was  considered.  The  total 
affect  of  the  ground  factor  on  tha  attenuation  provided  by  a 
building  will  ba  discussed  in  Chapter  3.  A  listing  and 
description  of  tha  program  that  finds  the  ground  factor  w  is 
givan  in  Appendix  A. 
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SUGGESTED  GROUND  EFFECT  DOSE  REDUCTION  FACTORS 


Surface 

w 

d/w 

dagraas 

Smooth  plan* 

a. 

1. 

43. 

Pavad  Araaa 

0.394 

1. 

45. 

Lawns 

2.370 

0.959 

43. 

Gravsllsd  Araas 

7.  468 

0. 897 

43. 

Plowad  field 

23.821 

0.839 

45. 
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2.7  It!3_Fallout_Modal 

Tha  -fallout  sourca  program  models  a  structure  as  having 
cylindrical  walls  and  a  flat  roof.  The  error  in  assuming 
tha  building  is  round  is  small  and  is  not  anargy  dapandsnt 
C523.  Tha  purpose  is  not  to  accurately  predict  the  dose 
reduction  factor  for  any  particular  building,  rather  to 
predict  how  dose  raduction  factor  will  change  with 

sources  of  different  energies.  The  reference  position  is 
taken  as  1  meter  off  tha  floor  in  the  cantar  of  the 
structure.  Infinity  is  takan  as  12  mfps.  Tha  program  was 
tested  to  prove  that  the  error  in  only  integrating  to  12 
mfps  is  less  than  0.001%.  The  program  does  three 
integrations!  One  for  the  reference  unprotected  exposure; 
one  for  the  exposure  from  material  deposited  on  the  roof; 
and  one  for  the  exposure  from  material  deposited  on  the 
ground. 


FIGURE  11  0MBttry_Bf_tbf_ElllBUi_ECSflCi!l!_Bfldfl 
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The  point  kernel  for  exposure  from  a  disk  sourca,  as  in 
th«  cas a  of  either  tha  roof,  floor,  or  ground  (aaa  Figure 
11) ,  ias 

G  <r)  =  0.5  k  3  r/(r3  +  a3)  ■  exp <-mf p ) ■ dr  <2.28) 
wharei  "r"  is  tha  incremental  radius  of  tha  disk; 

"a"  is  tha  length  of  the  normal  from  tha  canter 
of  tha  disk  to  the  datactori 
"B"  is  tha  appropriate  buildup  factor  for  shield 
<roof  or  wall)  and  air  combined; 

"k"  io  the  energy  dependent  conversion  constant 
that  convarts  "flux"  to  exposure; 
and,  "mfp"  is  the  total  number  of  mfps  from  the 
ring  of  integration  to  the  detector. 

The  kernels  are  integrated  using  a  form  of  Simpson's  rule 
using  rings  of  varying  width  to  reduce  computer  time. 

The  program  requires  the  input  of  the  "equivalent 
height",  "equivalent  radius”,  wall  mass  thickness,  roof  mass 
thickness  and  the  ground  factors  Y  (degrees)  ,  d/w,  and  w  as 
discussed  in  section  2.6.  Values  for  the  mass  thickness  of 
the  walls  and  roof  are  given  in  the  literature  C533.  The 
mass  thickness  of  lightly  constructed,  wood  frame  houses 
varies  from  about  5  to  17  gr/cm3  with  a  median  value  of 
about  10  gr/cm3.  For  brick  and  block  houses  the  mass 
thickness  varies  from  about  10  to  32  gr/cm3  with  a  median  of 
about  22  gr/cm3.  The  mass  thickness  of  a  1  foot  thick 
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concrete  wall  is  about  61  gr/cm3.  The  equivalent  height  o f 
a  structure  may  ba  -found  L»y  intagrating  tha  uncoil idad  gamma 
photon  path  langth  ovar  tha  atructura's  di mansions  or 
astimated  for  raasonably  square  structures  by  finding  tha 
radius  of  a  circle  with  an  equivalent  area  as  tha  structures 
floor  plan. 

The  program  calculates  tha  unprotected  exposure,  tha 
exposure  from  sources  on  the  roof  and  tha  axposura  from 
sources  on  the  ground.  It  than  calculates  the  DRF  for 
source  energies  from  15  KeV.  to  15  MeV.  and  for  the  RSS 
PWR-2  Fallout,  tha  One  Hour  Weapon  Fallout,  the  TMI  Fallout, 
and  the  5L-1  Fallout  Spectra.  Included  in  the  c-  tput  is  the 
exposure  that  would  be  received  from  fallout  internal  to  the 
structure  if  the  fallout  was  in  the  same  concentration  as 
outside  the  structure.  The  internal  exposure  assumes  no 
shielding.  However  internal  fallout  is  not  considered  in 
the  calculation  of  tha  DRFs. 

A  listing  and  description  of  the  Fallout  program  is 
given  in  Appendix  B.  Table  14  gives  a  typical  output  from 
the  Fallout  program.  The  results  of  the  Fallout  program 
will  be  discussed  in  the  next  chapter. 
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TABLE  14 

TYPICAL  OUTPUT  OF  FALLOUT  PRCBRAM 


Saill  Hood  Hcuit  on  Laan.  21  August  1934 


Integration  to  12  flFPi  in  ilr. 

Ditector  ii  1  altar  abovi  saooth  plant 
Equivalent  radius  of  structun  *  5.231  attars 

Equlvalant  hiight  of  structure  ■  5.231  eaters 

Hall  mi  thickness  •  11. 330  gr/ca2 

Rod  am  thicknm  ■  11.321  gr/c*2 

The  Ground  Fictori  ant 

The  trough  angle  PSI  «  45. Mil  degrees 

The  charictoriitie  trough  width  M  *  2. 3693 

The  ratio  of  flat  to  trough  D/F  *  .9591 


ElKiV)  UNPROTECTED 

ROOF 

GROUND 

TOTAL  E1P. 

ROOF  GOUT. 

DRF  RATIO  TO  1.12  Hr. 

INSITE 

1.12  Hr. 

14.83 

.9441 

5.287 

6.751 

.1393 

.4553 

1.338 

4,443 

RSS-Fallout 

14.64 

.9198 

5.752 

6.672 

.1379 

.4557 

1.211 

4.347 

Tfll-Fal  lout 

7.319 

.4716 

2.563 

3.332 

.1556 

.4319 

.94R5 

2.141 

SLl-Fallout  7.134 

.4779 

2.619 

3.087 

.1548 

.4327 

.9513 

2.175 

15331. 1 

133.4 

8.772 

63.32 

72.59 

.1ZM 

.5244 

1.152 

39.23 

12Ud. 3 

133.9 

6.482 

<6.77 

53.25 

.1217 

.5278 

1.159 

28.73 

0121.1 

84.19 

5.443 

38. 9C 

44.43 

.1226 

.5277 

1.159 

24.18 

6*11. 1 

67.71 

4.481 

31.18 

35.48 

.124* 

.5241 

1.151 

19.51 

5111.3 

59.49 

3.872 

27.32 

33.91 

.1253 

.5194 

1.141 

17.24 

4313.8 

51.59 

3.322 

22.81 

26.13 

.1271 

.3165 

1.134 

14.75 

3313.1 

41.61 

2.698 

13.23 

21.92 

.1291 

.5331 

1.135 

12.28 

2183.3 

31.63 

2.111 

13.25 

15.26 

.1318 

.4826 

1.361 

9.384 

1598.3 

25.91 

1.614 

11.48 

12.19 

.1334 

.4668 

1.825 

7.594 

1332.3 

23.42 

1.472 

9.398 

18.87 

.1354 

.4642 

1.119 

6.922 

<173.1 

21.11 

1.331 

3.341 

9.671 

.1375 

.4612 

1.111 

6.267 

1918.1 

18.34 

1.166 

7.152 

8.318 

.1412 

.4534 

.9958 

5.522 

391.1 

15.15 

.9711 

5.882 

6.773 

.1434 

.4478 

.9817 

4.633 

662.3 

12.75 

.8214 

4.825 

5.646 

.1455 

.4428 

.9726 

3.361 

631.1 

11.66 

.7511 

4.375 

5.126 

.1465 

.4395 

.9651 

3.524 

533.3 

9.644 

.6345 

3.584 

4.219 

.1514 

.4375 

.9617 

2.939 

491.1 

7.644 

.5135 

2.833 

3.316 

.1548 

.4338 

.9528 

2.335 

333.1 

5.533 

.3328 

1.982 

2.365 

,)M? 

.4274 

.9386 

1.696 

231.1 

3.688 

.2378 

1.289 

1.446 

.1631 

.3922 

.3614 

1.356 

153.3 

2.587 

.1573 

.7596 

.9169 

.1715 

.3545 

.7735 

.7277 

131.1 

1.819 

.72J4E-11 

.3661 

.4381 

.1644 

.2418 

.5238 

.4645 

81.3 

1.661 

•3876E-31 

.2199 

.2536 

.1499 

.1557 

..7423 

.4846 

61.1 

1.619 

. 1216E-I1 

.8242E-81 

•9447E-I1 

.1276 

.58712-11 

.1289 

.4861 

53.0 

1.629 

.3791E-82 

.2992E-I1 

.3371E-11 

.1125 

.2I69E-I1 

.4545E-31 

.4545 

41.1 

1.739 

.3695E-I3 

.3348E-I2 

.3718E-I2 

.9941E-I1 

.2137E-I2 

.4694E-02 

.5854 

33.8 

2.875 

•8638E-36 

. 6282E-05 

.7866E-I5 

.1222 

.34355-15 

.7478E-B5 

.9195 

21.1 

2.754 

.  3283E- 16 

.8188 

.321BE-I6 

1.131 

, 1165E-I6 

•2559E-16 

1.769 

15.3 

3.145 

.£631 

.8118 

.1881 

.88B9E+33 

.3811 

.8388 

2.541 
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2.8  It3e_Clgud_§gyr£f_Mgg§l 

The  cloud  source  program  used  here  models  the  structure 
as  a  hemisphere  (see  Figure  12).  The  exact  height  of  the 
reference  position  for  this  volume  source  is  of  little  or  no 
consequence  C543.  Therefore  the  reference  position  is  taken 
as  ground  level.  The  basic  source  configuration  is  taken  as 
semi -inf ini te  because  the  ground  excludes  its  lower  half. 
Infinity  is  taken  as  12  mfps.  The  program  was  tested  to 
prove  that  the  error  <n  only  integrating  to  12  mfps  is  less 
than  0.01%.  The  error  would  be  less  than  0.0007%  except 
that  the  buildup  factors  are  large  at  large  mfps.  There  are 
no  ground  roughness  or  terrain  effects  for  a  cloud  source. 
The  "Cloud"  program  (Appendix  C)  assumes  the  radioactive 
material  is  uniformly  distributed  in  the  cloud.  This 
assumption  may  be  relaxed  when  the  program  is  combined  with 
a  dispersion  model. 

The  point  kernel  for  a  hemispherical  structure  in  a 
semi-inf ini te  uniformly  distributed  cloud  source  is 
particularly  simplei 
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Which  reduces  toi 

E  *  k/2  axp(-m-fpw) 


°o 


Bcomb  •  exp  (-  /*. ■  r )  •  dr 


(2.3 a> 


Tha  buildup  -factor  has  throe  regions:  (1)  tha  air  between 
tha  hemispherical  cloud  shall  (sea  F’gu'-e  12)  and  tha 
building;  (2)  tha  structure  doma  and;  (3)  tha  interior  air. 
Thus  as  discussed  in  section  2.4,  tha  combined  buildup 
■factor  becomes! 


3ccmb  =  Baib (r )  ■  3wall !r  +  t) /3wall (r ) 

■  3  „-r«(r  +  t  +  a'/3Ai«(r  +  t).  (2.31) 

Hin  "r",  "t" ,  and  "a“  represent  tha  air  between  tha 
.structure  and  the  shall  of  integration,  the  roof  and  walls, 
and  the  intarior  air  respectively.  3« i «  is  tha  buildup 
factor  for  air  and  Buaul  is  the  buildup  factor  for  the 
walls.  Tha  thickness  of  the  walls  in  mean  free  paths  is 
represented  by  mfpw. 

Equation  2.30  is  integrated  using  Simpson's  rule.  The 
integration  is  done  in  regions  varying  from  3.2  mfp  to  3.0 
mfps  in  width  to  reduce  computer  time.  The  total 
integration  error  is  lass  than  0.01%.  Table  13  gives 
typical  results  for  a  structure  with  an  equivalent  radius  of 
5  meters  and  walls  10  gr/cm3  thick. 


Bunsen  and  Profio  presented  data  showing  that  far 
lightly  constructed  wood  frame  houses  the  wall  mas’s 
thickness  varies  from  about  5  to  17  gr/cm2  with  a  median 
value  of  about  10  gr/cm2j  and  for  brick  and  block  houses 
from  about  10  to  32  gr/cm2  with  a  median  valua  of  about  22 
gr/cm2  C333.  Thay  also  assumed  that  the  roof  and  wall  ara 
of  equal  thicknass  -  which  seems  raaaonabla  for  naw,  energy 
efficient  construction.  Thair  data  ware  used  for  their 
ground  source  model,  and  should  be  valid  for  a  cloud  source 
also.  However,  astimates  as  low  as  3.4  gr/cm2  appear  in  the 
literature  C5&3.  Tha  10  gr/cm2  figura  was  chosen  as 
raprasantati va  of  small  wood  frame  houses  and  22  gr/cm2  as 
reprosentati ve  of  brick  faced  houses.  This  discrepancy  is 
mentioned  not  to  pass  judgment,  but  to  explain  why  the 
numbers  arrived  at  hare  may  be  diffarant  than  those  baaed  on 
other  works  C37,3B3.  It  is  not  the  purpose  of  this  paper  to 
calculate  protection  factors,  but  to  show  how  they  vary  with 
energy. 

This  paper  is  not  the  first  to  use  a  hemispherical 
shell  modal  to  calculate  protection  factors  for  cloud 
sourcas  using  point  karnels  [343,  There  ara  several 
important  differences  however.  First  of  all,  the  modal  used 
here  integrates  to  12  mfpsf  the  referenced  modal  only  goes 
to  3  mfps.  For  a  disk  source,  3  mfps  might  be  far  enough, 
however  with  a  volume  source  there  is  no  1/R  factor,  thus  3 


67 


mf  ps  excludes  more  than  3%  of  the  source  Cexp<-3)  *  .04903. 
When  the  high  buildup  factors  from  these  distances  are 
included,  the  error  becomes  significant.  Second,  the 
previous  model  uses  the  buildup  factors  for  Mater  for  the 
Malls,  this  model  uses  concrete.  Third,  the  previous  model 
uses  the  linear  formula  for  calculating  the  buildup  factors 
for  the  air  outside  the  structure.  As  already  discussed, 
that  formula  is  only  good  for  thin  shields.  Fourth,  that 
model  multiplies  the  tMO  buildup  factors  together,  rather 
than  combining  them  -  thus  introducing  a  conservatism. 

Fifth,  the  previous  model  ignores  buildup  and  attenuation  by 
air  interior  to  the  structure. 

The  hemispherical  cloud  model  is  easily  extended  to 
include  cases  Mhere  the  cloud  is  limited  in  size,  either  in 
height  by  an  elevated  inversion  layer  or  similar  atmospheric 
phenomenon)  or  in  radius  as  it  might  be  if  the  structure 
Mere  near  the  source.  In  the  first  case,  the  integration  is 
carried  out  as  Mith  an  infinite  cloud,  however  those  shells 
Csee  Figure  123  that  are  partially  outside  the  height  limit 
are  modified  by  a  factor  of  h/r,  where  h  is  the  height  of 
the  top  of  the  cloud,  and  r  is  the  radius  of  the  shell. 

This  factor  gives  the  fraction  of  the  volume  of  the 
hemispherical  shell  beloM  the  height  of  the  inversion  layer 
(h).  In  the  second  case  the  integration  is  terminated  at 
the  cloud's  radius.  Both  methods  may  be  combined  for  clouds 


with  haight3  that  ar»  mora  limited  than  thair  radii. 

A  description  and  program  Hating  for  thia  modal  is 
givan  in  Appondix  C.  Table  13  contains  an  axampla  o-f  a 
typical  program  listing.  Ths  rssults  o-f  this  program  are 
givan  and  discussed  in  tha  naxt  chapter. 
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TABLE  13 

BSSUUia-QE.IbE_CLQUB-BQUBCE-HQBEL 


Typical  Bull  Hoed  Frui  Houn  in  infinits  cloud.  I  August  84 

tntsgritioA  to  12  RfPs  in  sir, 

Equivilsnt  radius  of  structurs  *  3.988  Mtsrs 

data  thlcknsss  of  sills  *  II. IM  gr/cs**2 
Rnliufl  cloud  height  ■  hmsshhm  niters 
Huiiui  cloud  rsdius  ■  mhhhhm  uteri 


E(KsV) 

8/hr/Ci*a3 

R/hf/Ci«i3 

l/hr/Ci«t3 

8Rf 

R/R1R88) 

Exclusion 

Unprotsctsd 

Protsctid 

fro*  int. 

Factor 

888-CLOUD 

348.1 

372.7 

1.922 

.4982 

l.lll 

.9744 

i.12  HR. 

814.2 

373.4 

13.44 

.7143 

1.121 

.9748 

THI-CL0U8 

184.3 

37.13 

1.832 

.3343 

.7771 

.9484 

13N4.I 

•1277E+I3 

•1I42E+I3 

128.4 

.1312 

1.214 

.9888 

1NN.I 

8188. 

7393. 

88.19 

.8299 

1.212 

.9882 

8811. 1 

7144. 

3984. 

73.94 

.8247 

1.199 

.9874 

6HI.I 

3438. 

4433. 

59.84 

.B192 

1.187 

.9847 

5811. 1 

4371. 

3712. 

32.92 

.8121 

1.177 

.9839 

48N.I 

3487. 

2939. 

43.27 

.8824 

1.143 

.9849 

3188.1 

2829. 

2225. 

37.45 

.7844 

1.139 

.9834 

21  HI 

1939. 

1483. 

28.34 

.7388 

1.898 

.9811 

13N.I 

1512. 

1112. 

23.31 

.7332 

1.143 

.9793 

1332.1 

1324. 

943.3 

21.24 

.7244 

1.833 

.9784 

1173.1 

1133. 

124.2 

19.24 

.7143 

1.138 

.9772 

1IM.I 

949.4 

488.3 

14.98 

.7117 

1.817 

.9754 

8N.I 

771.1 

324.1 

14.12 

.4131 

.9897 

.9739 

442.1 

439.4 

424.4 

11.83 

.4471 

.944! 

.9738 

4N.I 

381.2 

382.2 

11.81 

.4575 

.9327 

.9723 

SM.I 

471.8 

312.1 

9.821 

.4414 

.9278 

.9711 

4N.I 

371.9 

229.3 

7.147 

.4144 

.8931 

.9497 

388.1 

271.8 

134.9 

3.213 

.3794 

.8395 

.9479 

281.1 

187.8 

94.41 

3.241 

.5134 

.7439 

.9475 

158.1 

133.1 

41.23 

2.235 

.4441 

.4444 

.9442 

188.8 

91.79 

29.34 

1.424 

.3228 

.4443 

.9348 

88.8 

75.21 

17.34 

1.242 

.2313 

.3348 

.9332 

48.8 

37.77 

3.993 

1.243 

.1137 

.1513 

.8288 

58.8 

47.78 

1.947 

1.393 

•4123E-81 

.5974E-81 

.3853 

48.8 

34.79 

.1831 

1.794 

.5831E-82 

.7289E-I2 

.9353E-I1 

38.8 

24.43 

.4121E-83 

2.817 

•1339E-84 

. 2239E-84 

.1442E-83 

28.8 

17.28 

.1384E-13 

3.434 

. B822E-13 

.1142E-14 

•2349E-14 

15.8 

13.24 

.81H 

7.899 

.8818 

.8888 

.8811 
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2. 9  5°yCC33_Ioi2Ci2C-.i3_J_§icy£iy!I3 

Whan  radioactivs  matarial*  infiltrate  into  a  atructura, 
they  contributa  to  tha  axpoaura  received  by  aomaona  inaida 
tha  atructura.  If  a  passing  radioactiva  cloud  dapositad 
particulata  matarial,  than  thera  could  ba  two  aourcaa 
intarior  to  a  structural  tha  particulata  matarial  that  is 
dapoaitad  on  tha  floor  of  the  atructura  and  an  intarior 
cloud  of  gacaous  matarial. 

Infiltration  is  tha  uncontrolled  flow  of  air  through  a 
building  as  opposed  to  ventilation  which  is  controlled  flow 
(that  is,  it  can  be  shut  off).  For  a  typical  building, 
infiltration  causes  from  0.5  to  2.  air  changes  par  hour 
C593.  Infiltration  is  proportional  to  the  indoor-outdoor 
temperature  difference  and  windspeed  as  wall  as  being  a 
function  of  building  construction  and  condition  159,603. 

Therefore,  the  rats  of  increase  of  radioactive  gases 
inside  a  building  is  a  function  of  several  parameters  in 
addition  of  the  outside  concentration.  However,  unless 
radioactive  decay  is  rapid,  the  inside  concentration  will  be 
essentially  equal  to  the  outside  concentration  in  a  manner 
of  minutes  to  at  most  a  few  hours. 

The  infiltration  of  particulate  matter  must  also  be 
considered.  It  has  been  shown  that  the  amount  of 
particulate  matter  in  a  building  is  roughly  the  same  as 
outside  the  building  C613.  However  that  work  did  not 
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determine  how  much  o-f  the  particulate  matter  cam*  -from  the 
outdoor  anvironmant.  Tha  amount  coming  from  tha  outmida  may 
ba  vary  amall.  It  was  shown  by  maaauramants  takan  aftar  tha 
1957  Windscala  incident  that  "tha  deposition  within  a 
building  of  iodine  .  .  .  from  a  cloud  passing  around  it  is 
generally  only  1  or  2  percent  of  tha  deposition  on  the 
surrounding  ground"  C623. 

Tha  remainder  of  this  section  deals  with  tha  exposure 
recaived  from  sources  inside  a  building  that  are  of  tha  same 
concentration  as  outside.  True  exposure  rates  can  be  found 
from  these  data  by  calculating  the  actual  interior 
concentrations  from  infiltration  data. 
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2.9.1  Int^r  l.cr_F  3l.l_gut 

Tha  Fallout  modal  also  calculataa  tha  axposura 
racaivad  from  a  fallout  aourca  interior  to  a  structure.  It 
assumes  that  tha  concantrati on  of  dapasited  matarial  ia  the 
sama  as  that  extarior  to  tha  structure.  Tha  program  allows 
no  shielding  for  the  intarior  source.  The  exposure  is 
calculated  using  the  same  routine  as  for  the  infinite 
unshielded  smooth  plane  source  that  ia  used  as  the  standard 
for  calculating  DRFs.  The  only  difference  is  that  tha 
source  is  limited  to  the  intarior  radius  of  tha  structure. 
Figure  13  gives  the  ratio  of  the  exposure  received  from 
intarior  fallout  to  the  sum  of  the  exposure  from  both 
interior  and  exterior  fallout,  assuming  that  the 
concentrations  are  the  same  for  both.  Whan  this  ratio  is 
equal  to  1.8,  all  the  exposure  is  from  the  intarior  fallout. 
When  it  is  equal  to  0.,  all  tha  exposure  is  from  external 
sources. 

As  shown  in  Figure  13,  if  intarior  fallout  is  of  tha 
same  order  of  magnitude  as  extarior  fallout,  it  is  always 
important.  Below  approximately  130  KeV,  intarior  fallout 
can  account  for  nearly  all  the  exposura. 
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2.9.2  Intariar_Cl oud 

Tha  Cloud  modal  also  calculatas  tha  axposura 
r^caivad  -from  a  cloud  aourC3  interior  to  a  structure.  It 
also  assumes  that  tha  concantration  in  tn»  interior  cloud  is 
tha  same  as  that  of  tha  axta^ior  cloud.  Tha  lntarior  cloud 
is  also  used  to  calculate  tha  "exclusion  -factor"  which  is 
useful  in  datarmining  tha  importanca  of  infiltration  to  tha 
total  dose  racaivad.  Tha  axclusion  factor  is  defined  as  tha 
ratio  of  tha  axposura  from  tha  axtarior  cloud  to  tha  sum  of 
the  exposurss  from  tha  intarior  and  axtarior  clouds.  An 
axclusion  factor  of  1  maans  tha  intarior  cloud  contributas 
nothing  to  tha  total  axposura.  A  vary  small  axclusion 
factor  means  nearly  all  the  axposura  comas  from  tha  internal 
source.  Figure  14  shows  tha  value  of  one  minus  the 
exclusion  factor  as  tha  source  energy  varies.  One  minus  tha 
exclusion  factor  is  the  ratio  of  the  exposure  from  tha 
intarior  cloud  to  the  total  exposure. 

As  the  energy  of  the  source  photons  decrease  so  does 
thair  mfp.  Thus,  at  low  energies  the  interior  source  is 
vary  important,  but  at  higher  energies  it  is  always 
insignificant  as  shown  in  Figure  14. 
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3.  1  Introduction 

Tha  rtjsuiu*  that  follow  have  baan  obtainad  from  tha 
programs  that  ara  documented  in  tha  Append! ass.  They  ara 
not  tha  only  results  that  can  ba  obtainad.  Tha  programs  ara 
designed  to  allow  a  vary  wida  ranga  of  input  data.  Tha 
Fallout  and  Cloud  programs  ara  intandad  to  provide  DRFs  that 
can  ba  usad  as  ratios  to  translate  known  DRFs  to  naw 
spectra.  If  a  particular  structure  had  a  known  DRF  for  tha 
1-hr.  bomb  spectra  but,  for  example  tha  DRF  was  needed  for 
tha  TMI  spectra;  similar  structura  dimensions  should  be  used 
as  input  to  the  programs.  The  results  would  includa  a  1-hr. 
bomb  spectra  and  a  TMI  spectra  value.  Tha  ratio  of  these 
two  values  should  be  multiplied  by  tha  known  value  to  obtain 
the  required  DRF  for  the  structure.  However,  tha  results 
found  by  the  programs  are  very  close,  often  within  1 
percent,  to  tha  valuas  found  for  similar  structures  in  the 
literature.  Therefore  the  results  reported  here  do  not 
includa  any  such  translations.  DRFs  for  other  spectra  can 
be  obtained  from  the  valuas  for  specific  energies  that  are 
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given  as  part  o f  tha  program  outputs. 

In  ordar  to  compara  rasults  for  specific  structura 
typas,  tha  'following  structuras  ara  defined  from  valuas 
found  in  tha  litaratura  C551.  A  "small  wood  house"  has  a 
roof  and  walls  of  1C!  gr/cma  mass  thickness  and  a  S  m. 
equivalent  radius  and  height  (about  1000  sq.  ft.  of  floor 
area).  A  "large  brick  house"  has  22  gr/cma  walls  and  a  22 
gr/cma  roof  with  a  10  m.  equivalent  radius  (about  3000  sq. 
ft.).  A  "Small  Block  House"  is  1  meter  in  radius,  2  maters 
high  and  has  1  foot  thick  concmte  walls.  Tha  "thin  walled 
house"  has  walls  of  3.4  gr/cma  mass  thickness.  The  "thin 
walled  house"  is  included  to  allow  comparison  with  valuas 
found  in  tha  litaratura  £56,57,583  as  discussed  in  Section 
2.8. 

Tha  ground  roughness  surface  types  use  the  parameter 
valuas  reported  in  Table  13. 
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3.2  gffS£t_Sf„eC2yD^_RSiiShDf«l 

Ground  roughness  factors  arm  usad  in  the  currant 
shielding  litaratura  as  multipliars  far  DRFs  found  far 
structures  alone.  Hrwuvsr,  because  the  ground  roughnase 
changss  tns  spactra  of  tha  fallout  radiation  source,  and 
further  because  the  change  in  the  spectra  is  energy 
dependent,  the  ground  roughness  factor  cannot  be  separated 
from  the  total  DRF,  except  for  specific  spectra.  In  order 
to  3oa  the  snergy  dapandant  nature  of  tha  ground  roughness 
factors,  the  DRFa  obtained  for  a  structure  surrounded  by  a 
particular  ground  type  can  be  divided  by  the  DRFs  obtained 
for  the  same  structure  surrounded  by  the  hypothetical 
infinite  smooth  plane.  Figure  15  shows  the  effect  of 
Ground  Roughness  on  the  exposure  inside  a  small  wood  framed 
house.  Tha  figure  represents  the  ratio  of  the  DRF  found  for 
tha  house,  with  the  ground  taken  into  account,  to  the  DRF 
found  far  the  same  house  over  a  smooth  plana.  Thus  tha 
plotted  result  is  tha  DRF  associated  witn  the  ground  effect 
in  tha  way  it  is  usually  ussd,  that  iu  to  multiply  the  DRF 
found  far  the  structure  itsslf.  Figure  16  does  the  same  for 
a  large  brick  house. 

Tha  Lawn,  Gravel?  and  Plowad  Field  would  give  DRFs  of 
3. 80,  0.73  and  0.60  respectively  if  found  ■‘or  an  infinite 
plane  with  no  structure  present.  With  a  structure  present, 
the  DRFs  are  smaller  (the  affect  greater)  for  energies  above 
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3D  ;<.iV  -for  j  small  weed  houoa  and  33  KoV  -for  a  large  brick 
house.  Th  i  DRFa  are  smaller  bacausa  tha  structure  walls 
provide  greater  attenuation  -for  tha  photons  that  have  lost 
energy  in  interactions  with  ths  ground.  Bolow  the  30  and  50 
KaV  energies,  the  m-fps  o-f  tha  gammas  ara  so  short  that  the 
ground  plays  almost  no  role  in  tho  shielding. 
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3.  3  Ef  f3£t_of_3;_d  1  ding_3i23_cn_Fal  lout_E}<2oa|ura 

Figure  17  shows  the  energy  dependent  DRFs  for  "small 
wood  houses",  "larga  brick  houses",  and  "small  block 
houses."  For  a  small  wood  house  on  a  lawn,  the  DRF  is 
nearly  constant  for  energies  greater  than  300  KeV.  Below  40 
KeV,  the  DRFs  (and  therefore  the  doses)  ar®  near  zero.  With 
the  larga  brick  house  and  the  "block"  house,  the  DRFs 
increase  nearly  linearly  with  the  logarithm  of  energy  over  a 
wide  energy  range. 

Figures  18  and  19  demonstrate  the  effect  of  increasing 
structure  dimensions.  Figure  IS  shows  the  effect  of 
increasing  the  building  size  (effective  radius)  for  a  wood 
structure.  rt  doubling  of  the  building  radius  reduces  the 
DRFs  by  about  22‘i  uniformly  for  all  energies.  This  is 
because  increasing  the  building  size  effectively  removes 
part  of  the  source.  Figure  19  shows  the  effect  of 
increasing  the  wall  thickness  for  a  small  house.  Note  the 
sharp  drop  in  the  DRFs  at  lower  energies  and  with  increasing 
wall  thickness.  The  effect  of  increasing  the  wall  thickness 
is  much  more  energy  dependent  than  increasing  building  size. 
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FiGURE  17  BBE_y3i_g03CaY_f2C_5uil^iDS3_£iiS335^_tg_F«llgut 
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FIGURE  19  20E_V3i_lD3CaY_fSC_IO£!:3*fiog_y3ll_BJ3§_IbiEKD93» 
fQE_E3ll3yi_S3ynE3S  (S»0  matar  af-factiva  radiua) 
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3.4  §if3£i3_sf_Ql.syd_3sycs23 

The  effective  DRF  vs.  Enargy  curva  -for  four  diffarcsnt 
structures  is  given  in  Figure  23.  Tha  results  of  increasing 
tha  building  siza  are  pronouncad.  Figure  21  shows  tha 
•ffact  of  increasing  wall  thickness  on  tha  DRF  for  an 
infinite  cloud  source.  Tha  affsct  is  nearly  tha  same  as 
i«an  in  Figura  23.,  Figura  22  shows  that  affect  of 
increasing  tha  auilding  affective  radius  for  a  structure  is 
vary  small.  This  is  due  to  tha  nature  of  a  spherical 
source,  which  has  no  “l/R"  dependence. 

Figures  23  and  24  demonstrate  the  effect  of  limiting 
the  cloud  size.  Limited  cloud  size  has  very  little  affect 
an  the  DRFs.  However,  while  the  DRFs  are  not  effected,  tha 
doses  are  greatly  effected  as  shown  in  Figure  23.  The 
kernel  plotted  in  Figura  23  includes  tha  exposure  rate  term, 
which  is  energy  dependent.  The  effect  of  limiting  tha  cloud 
size  is  to  reduce  the  number  of  higher  energy  gamma  photons 
that  can  reach  the  structure.  Both  the  protected  and 
unprotected  exposures  are  limited  in  Figures  23  and  24. 
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FIGURE  21  P85_y3x>lDSC3S3iDfl_l!l«U_e?«aa>Ibi!i'*S0faa_fBC-f 
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CCNCLUSICNS 

Th®  contri buti ana  of  this  di 3  ;  irtati cn  t n  tha  currant 
literature  coma  in  throe  wayai  first,  in  tha  questions  it 
anawersj  second,  in  tha  modal*  davalopwu;  and  third,  in  tha 
results  of  tha  modal*. 

Tho  primary  raault  of  this  work  i*  it*  answering  of  the 
question,  "Can  tha  DRFa  found  for  ana  apoctra  b®  uamd  for 
castas  involving  other  spectra  without  serious  a-ror?" 

Tha  Reactor  Safety  Study  assumed  a  DRF  of  0.33  for  all 
photon*  from  all  iaotopam  either  from  cloud  or  fallout 
aourcas  C633.  Tha  results  roportad  hare  claarly  show  that 
this  assumption  is  not  strictly  in  accordance  with  raality. 
Other  works  have  used  tha  DRFa  found  for  weapon  fallout  to 
predict  the  outcome  of  a  reactor  accident  CS, 9, 10,113.  In 
all  these  works,  the  underlying  assumption  was  that  their 
numbers  wore  "conservative".  This  work  can  be  used  to  prove 
or  disprove  those  assumptions  and  make  the  results  closer  to 
real i ty. 

Table  16  gives  DRFs  found  for  fallout  using  the  model 
just  presented.  Table  17  does  the  same  for  cloud  sources. 
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TABLE  16 

DRFa  FOR  FALLOUT  SPECTRA 


Spactra  Structura  DRF 


l.  Hr 

Fal lout 

Small 

Houaa  on  Lawn 

13.455 

PWR-2 

<RSS) 

0.  456 

TMI  ~ 

Fallout 

3.  432 

SL-1  - 

■  Fallout 

0  a  433 

1.  Hr 

Larga 

Houaa  on  Lawn 

3.  134 

PWR-2 

3.  192 

TMI 

3.  123 

SL-1 

0.  125 

1.  Hr 

Smal  1 

Block  Housa/graval 

3.374 

PWR-2 

3.  385 

TMI 

3.321 

SL-1 

3.322 

1.  Hr 

Small 

houaa/  amooth  pi  Ana 

3.611 

PWR-2 

3.613 

TMI 

3.538 

SL-1 

0.589 

1.  Hr 

Larga 

Houaa/  plowad  Fiald 

0.  144 

PWR-2 

0.  149 

TMI 

0.  132 

SL-1 

0.  103 

« 
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TABLE  17 

DRFs  FDR  CLOUD  SPECTRA 


Spectra 

Structuro 

DRF 

PWR-2  (cloud) 

Smal  1 

Wood  HoUS» 

0.690 

TMI  (cloud) 

2.536 

PWR-2 

Largos 

Brick  Houaa 

0.  443 

TMI 

3.255 

PWR-2 

Smal  1 

“Block  Houaa" 

0.  138 

TMI 

3.336 

PWR-2 

Smal  1 

"Thin  Wall  ad  Houoa" 

0.371 

TMI 

0.  333 
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Tha  daca dm  old  question,  "la  tha  DRF  for  a  PWR-2  avant 
approx i mat al y  aqual  to  that  of  fallout  from  an  'atomic 
bomb'?",  is  answarad  -  yaa !  That  was  a  good  assumption. 

This  answar  may  seam  a  bit  surprising  until  tha  dafinition 
of  a  PWR-2  avant  la  racallad.  A  PWR-2  avant  (sea  tha 
introduction)  lnvolvas  tha  auddan  release  of  a  aignificant 
amount  of  cars  material  without  decay  tima  C161. 

Tha  DRFs  for  accidants  that  hava  happanad  (SL-1  and 
TMI)  ara  small ar  than  thosa  of  a  PWR-2  avant.  This  ia 
aspacially  tha  casa  for  largo  buildinga  with  thick  walls. 
This  is  bacauaa  moat  of  tha  radioactiva  matarial  ralaasad 
consisted  of  Xa-133  which  haa  a  81  KaV  photon  that  is 
comparati valy  aaay  to  attanuata. 

It  is  not  a  unlqua  obsarvation  that  largo  buildings 
attanuata  photons  mora  than  small  buildings,  howavar  tha 
magnituda  and  anargy  dapandanca  of  tha  differences  is 
uniqua.  Tha  valuaa  glvan  hara  can  influanca  decisions 
concarning  whara  to  aaak  ahaltar  and  whan  to  avacuata.  Tha 
values  found  hare  can  ba  axtendad  to  vehlclas  also.  A 
vahicla  haa  vary  llttla  matarial  for  shielding  (about  2 
gr/cma>  CUD  and  a  very  high  infiltration  rata  (about  IB  air 
changes  ovary  hour  at  low  speeds)  C643.  Therefore,  if  the 
public  ia  evacuated  while  a  radioactive  cloud  is  passing, 
they  may  receive  a  much  higher  dose  than  if  they  sought 
shelter.  This  work  should  influanca  that  decision.  In  a 
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particular  situation,  -for  example,  this  work  would  halp  show 
that  it  may  bo  oa-far  to  soak  shaltar  in  a  largo  offics 
building  or  factory  that  is  naar  by  than  it  is  to  attempt  to 
ovacuato  on  a  crowded  highway. 

Thara  were  othar  important  contributions  made  during 
the  development  of  this  work.  First,  there  is  now  a  simple, 
aaoily  understood  mathod  for  accurately  combining  buildup 
factors  in  mul t. i -regi on  problems  that  works  for  all 
reasonable  gamma  energies.  Pravious  methods  wars  either 
based  upon  conservatism  rather  than  accuracy,  or  had  a 
narrow  range  of  energies  for  which  they  were  applicable. 

Second,  there  is  now  a  modal  for  determining  the 
shielding  affect  of  ground  roughness  that  is  applicable  to 
all  spectra.  Pravious  models  ware  only  applicable  to  a 
particular  spectrum. 

Third,  there  is  now  a  model  for  determining  tha  change 
in  shielding  factors  with  changing  source  energy  for  a 
particular  shield.  Previously,  each  new  source  required  tha 
entire  shielding  problem  to  ba  re-solved,  now  the  problem 
needs  to  be  solved  only  once  —  and  the  results  for  all 
energies  can  be  determined. 

The  results  reported  in  this  work,  and  the  methods 
presented  for  extending  them  to  other  situations  are  also 
vary  important  contributions. 
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It  has  not  known  that  the  ground  ■’•oughnass  -factor  was 
essentially  anargy  independent  abova  230  KaV  (aao  Figures  15 
and  16) . 

It  was  not  known  how  th«  DRFs  providad  by  structures 
■for  -fallout  varied  with  energy  —  it  was  not  known,  for 
instance,  that  abova  200  KaV  tha  DRF  for  a  small  housa  is 
essentially  constant,  but  below  53  KaV  it  is  assantially 
zero  (s am  Figure  17).  It  was  not  krnswn  that  tha  DRF  for  a 
large  housa  varies  somawhat  linaarly  with  tha  logarithm  of 
ansrgy  over  a  large  anargy  range. 

It  was  not  reported  that  tha  DRF  of  a  structure  in  a 
cloud  saurca  was  assantially  indapandant  of  tha  building's 
size  (affective  radius).  It  was  not  reported  that  the  DRFs 
for  a  cloud  source  were  assantially  independent  of  «ny 
limitations  on  tha  clouds  height  or  radius.  That  is,  the 
same  DRFs  can  be  used  during  inversions  or  very  close  to  tha 
point  of  release,  and  during  all  stability  conditions 
(provided  the  clouds  radioactive  composition  remains 
uniform)  that  are  used  for  infinite  clouds. 

The  relative  importance  of  interior  fallout  and  cloud 
sources  as  a  function  of  energy  was  not  known.  It  is  now 
known  that  for  a  large  house,  at  energies  below  about  100 
KeV  assantially  all  the  dose  from  fallout  could  coma  from 
dirt  that  infiltrates  into  tha  housa,  rather  than  from 
exterior  fallout  (sea  Figure  13).  It  was  not  reported  that 
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interior  clouds  ara  unimportant  in  small  houses  for  nearly 
all  anergi as  axcapt  below  about  73  KaV  (339  Figura  14). 

In  summary,  a  modal  -for  ths  protaction  provided  by 
structural,  such  as  homes,  against  releases  of  airborne 
radioactive  materials,  as  a  -function  of  tha  gamma  anorgy 
spectrum  of  tha  released  material  has  bean  developed  and 
evaluated  --  that  was  the  objective  of  this  research. 


CHAPTER  V 


RECOMMENDATIONS 

This  modal  is  but  a  part  of  a  larger  work  which  should 
jj  complatad  and  should  include  a  dispersion  modal,  an 
evacuation  modal  and  an  infiltration  modal  (ia«  Figure  1). 
Tha  data  darivad  from  this  modal  should  than  ba  used  to 
pradict  tha  exposure  tha  public  could  expect  to  racaiva  from 
a  release  of  radioactive  material  of  any  kind.  This  should 
ba  done  on  a  real  time  basis  as  outlined  in  the  introduction 
so  that  informed  decisions  can  be  made  regarding  evacuating 
or  sheltering  the  public  in  order  to  raduce  exposure  and 
injury. 

A  simple  method  cf  measuring  the  effective  mass 
thickness  of  the  walls  of  existing  structures  should  bo 
developed  to  allow  rapid  surveying  of  populated  areas  where 
radioactive  releases  are  most  likely  to  happen. 

A  Monte-Carlo  study  should  be  conducted  to  determine 
the  exact  accuracy  and  limits  of  applicability  with  regard 
to  varying  Z  number  and  energy  of  the  method  presented  here 
for  combining  buildup  factors  in  multi-region  problems. 
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APPENDIX  A 


GROUND  FACTOR  PROGRAM  DESCRIPTION  AND  LISTING 

Tha  Ground  Factor  program  calculates  tha  valua  of  tha 
trough  width  w  (aaa  Saction  2.6.2)  that  gives  a  known  CRF 
for  a  specified  spectra.  Tha  ratio  d/w  and  tha  angle  must 
also  be  specified  (see  Equations  2.19  and  2.20).  Tha  value 
of  w  found  is  required  as  input  to  the  Fallout  program. 

When  the  Ground  Factor  program  is  run,  it  will  a3k  for 
input  data  as  shown  in  Figure  26.  After  asking  for  a  title 
and  date  for  the  output  listing,  the  program  asks  for  a 
"Radius  to  the  inside  of  an  infinite  disk  in  meters. "  This 
value  is  normally  zero,  indicating  integration  over  an 
infinite  plane.  A  value  other  than  zero  can  be  used  to 
study  the  effect  of  a  central  void  in  tha  source,  ouch  as 
that  caused  by  a  house.  The  next  question  asked  is  "Do  you 
wi ah  to  calculate  a  TERRAIN  FACTOR  from  a  known  DRF?" . 
Normally  the  answer  should  be  yes  (y) .  The  no  <n)  option 
allows  a  quick  listing  of  the  output  for  a  known  value  of  w. 
Only  a  "yes"  will  calculate  w.  The  program  can  presently 
find  w  for  values  of  the  DRF  found  for  the  Cs-137,  Co-60, 
PWR-2,  and  1-Hr  bomb  fallout  spectra.  After  choosing  the 
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prcpor  spactra,  you  must  input  tha  givan  DRF,  d/w  and  Y. 

The  pregram  then  itarativaly  finds  w  and  gives  tha  resultant 
DRF 's  for  energies  from  15  KaV  to  15  MaV  and  for  tha  PWR-2 
and  1-Hr  bemb  fallout  spactra  as  shown  in  Tabla  10.  Tha 
program  rafars  to  w  as  tha  "tarrain  factor"  and  Y as 
"thata" . 

Tha  program  is  writtan  in  £ji£!I232£i  FORTRAN  3.13.  This 
is  a  version  of  FORTRAN  77.  Thera  ara  no  non-standard 
FORTRAN  commands  in  tha  program.  Tha  program  was  run  on  a 

J2D»iO0i3i_k£O2yiDC  <^C)  which  uses  an 
2353  procassor  and  an  535Z  co-processor.  Tha  program  was 
run  under  tha  fjS-DOS  2.11  operating  aystam.  However ,  it 
should  run  on  any  computer  capabla  of  using  FORTRAN  77.  Tha 
program  takas  approx i matal y  33  minutas  running  time. 

Caution,  tha  program  is  writtan  is  FORTRAN,  tha  input 
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The  MAIN  program  directs  tha  input  and  output  of  the 
program  data.  Subrouting  SIMPS  carrias  cut  th«  integrations 
by  using  Simpson's  rula.  Tha  nu/nbar  of  incr amenta  its 
automatically  increased  until  tha  dasirad  accuracy  is 
obtained.  Subroutine  PROTECT  speeds  the  integration  of  tha 
disk  source  by  dividing  the  disk  into  rings  which  converge 
to  integrals  -faster.  Function  DISK  -finds  tha  kornel  -for  tha 
hypothetical  infinita  smooth  disk  sourca.  Function  GROUND 
calculates  tha  karnel  for  the  fallout  using  tha  methods 
dr.vol cpad  in  Section  2.  □.  2. 

Function  AMU  finds  the  mass  attenuation  coaficisnt  for 
air.  Function  BUF  finds  the  buildup  factor  for  air. 

Function  MUCONC  finds  the  mass  attenuat.or.  coeficiant  for 
tha  ground.  Function  BCCNC  finds  the  buildup  factor  for  tha 
ground.  Function  ROENT0  finds  the  factor  tar  converting 
flux  to  exposure. 

The  listing  of  tha  Ground  Factor  program  follows. 


PROGRAM  GRDFCT 
C 

CHARACTER  CT I TLE*  1 16 , YORN* 1 
INTEGER  IN,  OUT 
REAL  E  (29) ,  MU,  MFP1,  LININT 
REAL  EGRCUP  (23,4) ,  FRACT (23,4) 

C  EXTERNAL  DISK,  GROUND 

COMMON  /ENERGY/  EKEV,  MU,  BU,  W , THETA , DOVERW 

COMMON  /INO'JT/  IN,  OUT 

COMMON  /OLD J /  EOLD,  JOLD 

COMMON  / 0 L D J 2 /  E0LD2 ,  J0LD2 

COMMON  / F  A G /  RMAX 

COMMON/KERN AL/TOTKRN, TQTDOS , TD I RT,  MF'P  1  ,D2,R1 
CAAAAAA  NOTE  RADII  ARE  IN  METERS,  ENERGIES  IN  KeV,  DENSITIES  IN  gr/crc*2 
C 

DATA  E/ 15. ,23. ,30. ,43. ,50. ,60. ,80. , 130. , 150. ,200. ,300, , 

1  430. ,503. ,603. , 662. , B30 . , 1 030 . , 1  1 73 . , 1 332 . ,  1 500 .  , 

2  2300. ,3333. ,4000.  ,5333. ,6030. ,8303. ,10033, , 15030.  / 

DATA  EGROUP  /661 .643, 19*0. , 1  173.  2, 1332. 5, 18*0.  , 

1  53. ,333.  ,753.  ,  1533. ,2533. ,15*3. , 

2  25. ,75.  ,153. .250. ,350. ,530. ,730. ,933. , 1  165.  , 

3  1500. ,1830. ,2253. ,2753. ,3530. ,4503. ,5*0./ 

DATA  FRACT  /I .  ,19*0. , 

1  1 . , 1 .  ,  19*3.  , 

2  . 1 696,. 2431,. 2763,. 21 83,. 0933, 15*3. , 

3  .0271  ,. 01 37,. 0737,. 0476,. 0929,. 1373,. 1717,. 1  627, 

4  .3839, .0957,. 3299,. 0397, .3143, .0042,. 0301 ,5*0. / 

C*******  LINEAR  INTERPOLATION  FUNCTION 

LININT  <X,X1,X2,Y1,Y2)  *  ( Y2-Y1 ) * ( X-X 1 ) / ( X2-X 1 )  +  Y1 
CAA'AAAAAATTEMPT  TO  CONTROL  PRINTER  FROM  FORTRAN  TO  COMPRESS  PRINT 
OPEN  (1  ,FIL£='LPT1 '.STATUS*' NEW' ,ACCESS=' DIRECT' 

1  ,F0Rf1=  'UNFORMATTED'  ,RECL*1 ) 

WRITE  (1)27,15 
CLOSE  (1) 

OPEN  1 1 ,  F ILE=  ' LPT  1 ' ) 

C 

C **************** '**********#************************************** 


C*f*************INPUT  TITLE  FOR  EACH  LISTING*********************** 
337  FORMAT  '  1H3 , 32X , 'SEMI-INFINITE  DI5K  SOURCE  IN  AIR', 

1  '  WITH  NEW  GROUND  FACTOR  USING  Tave.  13  JULY  1984') 

C******#**#***********************#***********#****#****#*****#**»* 

OPEN  < 1 , F I LE=  '  LPT  1 ' ) 

C 


OUT  *  1 
EOLD  «  0.0 
E0LD2  =  0.0 
IN  «  0 
NSPEC  *  0 

cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

WRITE!*, 290) 

290  FORNATflX, 'INPUT  TITLE/DATE  OF  OUTPUT') 
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READ  ( *  ,  29 1 ) CT I TLE 

291  FORMAT ( A 1 16) 

292  FORMAT  ( 1 H0 , 1 5X , A 1 1 6 ) 

C . 

WR I  TE  <  OUT  ,  S37 ) 

WRITE (OUT ,292) CTITLE 
WRITE!  OUT , 808 ) 

WRITE! OUT , 804 ) 

300  WRITE (*, 301 ) 

301  FORMATdX, 'INPUT  RADIUS  TO  INSIDE  OF  INF  DISK  IN  METERS  ' ,  \ ) 

302  READ  (*,303)  RR1 

303  FORMAT (F13. 6) 

C*  .***********#****#****************#********#*********#******* 

C*****ROUTINE  FOR  FINDING  TERRAIN  FROM  KNOWN  SPECTRA 

729  WRITE!*, 730) 

730  FORMATdX, 'Do  you  wish  to  calculate  a  TERRAIN  FACTOR  ', 

1  '■from  a  known  DRF?  (y  or  n)',\) 

READ ( * , 731 ) YORN 

731  FORMAT  ( A 1 ) 

IF (YORN. EQ. 'Y'.OR. YORN.EQ. 'y')  THEN 
GOTO  750 

ELSE  IF (YORN.EQ.  '  N  ' . OR . YORN. EQ . 'n')  THEN 
GOTO  722 
ELSE 

GOTO  729 
END  I F 

r33  W;  TTI  -*.732) 

73?  FORMAT d X ,' Choose  the  spectrum  from  which  your  DRF  was  optained:' 

1  Cs-137  [input  13',/,'  Co-60  [input  23',/, 

2  '  Ruactor  Safety  Study  [input  33',/, 

3  '  1.12  Hr.  Bomb  deposition  [input  43',/, 

4  '  "’ther  -limit  20  groups  [input  53  ?  ' ,  \ ) 

E A !  (  • ,  733 )  NSPEC 

733  FORMA  i(  ID 

IF (NSPEC. GE. 1 . AND. NSPEC. LE. 5)  GOTO  734 
GOTO  750 

734  IF (NSPEC. LE. 4)  GOTO  760 
WR I TE  <  * , 799 ) 

799  FORMATdX, 'THIS  FEATURE  IS  NOT  IMPLEMENTED  YET  -  TRY  AGAIN') 

GOTO  750 
760  CONTINUE 

735  WRITE!*, 736) 

736  FORMATdX, 'Input  given  DRF  ’,\) 

READ  <* , 303) GIVEN 

722  WRITE!*, 766) 

766  FORMAT  ( 1 X ,'  Input  given  D  over  W  ' , \ ) 

READ  (*,303) DOVERW 

WRI TE ( * , 767) 

767  FORMAT (IX ,' Input  given  theta  ' , \ ) 

READ ( * , 303) THETA 


*~rr^r*a*?** r* 
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IF ( YORN. EQ. ' N ' . OR. YORN. EQ.  ' n  ' )  GOTO  721 

r i  *  a. 

N3  -  0 
GUESS  1  =  1. 

W1  =  0. 

W2  =  2.0 
143  =  W2 

700  CONTINUE 
SUNDEN  *  0. 

SUMNOM  =  0. 

W  =  W3 

WRITE (*,777)W,GUE3S2 
777  FORMAT  ( 1H  ,/,'  GUESS=  '  , G 1 4 . 7  , 

2  '  PROTECTION  FACTOR  =  ' , G 1 4 . 7) 

DO  701  IEKEV  =  1,20 
EKEV  =  EGROUP ( I  EKE  V , NSPEC ) 

FRCT  =  FRACT ( I EKEV , NSPEC) 

RMA X  *  1.E3 

C . RESET  RMAX  IN  FLAG  EACH  TIME  ENERGY  IS  CHANGED 

I F  ( FRCT . EQ . 0 )  GOTO  701 
WRITE ( * , 806 )  EKEV 

C . LET  YOU  KNOW  PROGRAM  IS  RUNNING 

CALL  PROTECT  (DRF) 

WRITE!*, 702 )  DRF 
SUMDEN  *  SUMDEN  +  TOTKRN  *  FRCT 
SUMNOM  «  SUMNOM  ♦  TDIRT  *  FRCT 
702  FORMAT  (IX,  '  PF  =  ',614, 7, X) 

701  CONTINUE 

IF (N3. EQ. 0)  THEN 

GUESS2  »  SUMNQM/SUMDEN 
N3  =  1 

W3  =  LININT (GIVEN, GUESS], GUES32,W1,W2) 

GOTO  700 
ELSE 

GUE553  ■  SUMNOM/S'JMDEN 
END  IF 

IF (ABS (GUESS3-GIVEN) . LT. i.E-4)  GOTO  720 
GUESS  1  «  GUESS2 
GUESS2  «  GUESS3 
W1  =  W2 
W2  =  W3 

W3  «  LININT (GIVEN, GUESS1,GUESS2,W1,W2) 

GOTO  700 

720  CONTINUE 

WRITE  (*,777) W,GUESS3 

721  CONTINUE 

C***************************** ***************************************** 

I F ( YORN . EQ .  'Y'.OR.YORN.EQ. ' y ' >  GOTO  740 
WRITE (*,304) 

READ  ( * , 303 )  W 

304  FORMATdX,  'INPUT  TERRAIN  FACTOR  *,\) 
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740  CONTINUE 
R 1  =  RR1 

WRITE ( OUT , 809 ) R 1 

WRITE (OUT , 820 )  W 

IF (NSPEC. EQ. 0)  WR I TE ( OUT , 79 1 ) 

I F ( NSPEC . EQ .  1 )  WRITE  (OUT ,  792)  GIVEN 
I F (NSPEC . EQ . 2 )  WR I TE  ( OUT ,  793 )  GIVEN 
IF (NSPEC. EQ. 3)  WR I TE  (OUT ,  794 )  GIVEN 
IF(NSP:(\EQ.4)  WRI  TE  ( OUT ,  795)  GIVEN 
IF (NSPEC. EQ. 5)  WR I TE  (OUT  ,  796 ) 

WRITE (OUT , 797)  DOVERW 
WRITE (OUT ,798)  THETA 

791  FORMAT ( 16X ,' Terrai n  Factor  given  as  input') 

792  FORMAT ( 16X Terr ai n  Factor  calculated  from  Cs-137  spectrum.', 
1  '  Gi ven  ■  ' , F 1 0. 4) 

793  FORMAT ( 16X Terrai n  Factor  calculated  from  Co-60  spectrum.', 

1  '  Given  *  ',F10.4> 

794  FORMAT ( 16X ,' Terr ai n  Factor  calculated  from  RSS  spectrum.', 


1  '  Given  =  '  , FI 0. 4) 

795  FORMAT! 16X, 'Terrain  Factor  calculated  from  1.12  Hr.', 

1  '  fallout  spectrum.  Given  ■  ' ,  F 1 0 . 4 ) 

796  FORMAT ( 16X ,' Terr ai n  Factor  calculated  from  given  spectrum') 

798  FORMAT ( 16X ,' Theta  (ground  sawtooth  angle)  ■  ' , F 1 0 . 4 , '  degrees') 

797  FORMAT ( 16X ,' D  over  W  (ratio  of  flat  surface  to  groove)  * 


1  F10.4) 

411  CONTINUE 

WR I TE ( OUT , 803 ) 

cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

CXXXXX  ROUTINE  FOR  GENERATING  ENERGY/DRF  TABLE 

c xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


DO  80 1  IEKEV  »  28, 1,-1 
EKEV  *  E(IEKEV) 

WR I TE ( * , 806 )  EKEV 

806  FORMAT (IX,  '  E  ( KeV)  =  ' , G 1 3 . 5 , \ ) 

RMAX  »  1.E8 

C . RESET  RMAX  IN  FLAG  EACH  TIME  ENERGY  IS  CHANGED 

C . LET  YOU  KNOW  PROGRAM  IS  RUNNING 

CALL  PROTECT  (DRF) 

WRITE (OUT, 802)  EK£V,T0TKRN,T0TD0S,MFP1,D2, 

1  TD I RT , ORF 
801  CONTINUE 

800  CONTINUE 

C8888889898a9a48888a88888888«88a88988998888898088888f980899898@9@9@@0Q9 
089989  ROUTINE  TO  FIND  DRFS  FOR  SPECTRA 
DO  300  NSPEC  *  3,4 
SUMDEN  *  0. 

SUMNOM  ■  0. 

SUMDOS  »  0. 

DO  501  IEKEV  «  1,20 
EKEV  »  EGROUPdEKEV, NSPEC) 

RMAX  «  1.E8 


4 
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C - RESET  RMAX  IN  FLAG  EACH  TIME  ENERGY  IS  CHANGED 

FRCT  =  FRACT  ( I EKEV , NSPEC ) 

IF(FRCT.EQ.3)  GOTO  331 
WRITE!*, 806)  EKEV 

C . LET  YOU  KNOW  PROGRAM  IS  RUNNING 

CALL  PROTECT  ( DRF ) 

WRITE!*, 732!  DRF 
SUMDEN  *  SUMDEN  +  TQTKRN  *  FRCT 
SUMNOM  *  SUMNOM  +  TDIRT  •  FRCT 
SUMDOS  =  SUMDOS  *  TOTDOS  *  FRCT 
501  CONTINUE 

DRF  =  SUMNOM/SUMDEN 

IF (NSPEC.EQ. 4) WRITE (OUT, 535) SUMDEN, SUMDOS,  SUMNOM, DRF 
IF (NSPEC. EO. 3) WRITE (OUT, 510) SUMDEN, SUMDOS, SUMNOM, DRF 
535  FORMAT ( IX , 15X, ' 1.  12  Hr  .',  5X , 2 <3 1 0. 5 , 3X >,'*****»****', 3X , 

1  '*****♦*♦**', 3X,G10. 5, 6X,G10. 5) 

513  FORMAT ( 1 X , 1 5X , ' RSS  ',  5X , 2 ( G 1 0 . 5 , 3X ) ,'*♦********', 3X , 

1  '***♦**♦***', 3X,G13. 5, 6X,G10. 5) 

500  CONTINUE 

322  FORMAT!  IX,  15X ,6 (G13. 5,3X) ,3X,G10.5) 

933  FORMAT ( 1 H0 , 15X, 'E(KeV) ' ,6X. 

1  '  KERNAL  '  ,6X, ' R/hr /Ci *m2 ' , 4X , 

2  '  MFP l ' , 4X , ' D2  (maters) ',2X, 'TERRAIN  KERNAL', 9X , 'DRF') 

829  FORMAT ( 1 H0 , 15X ,' Intsgr at i on  to  12  MFPs  in  air.') 

809  FORMAT  ( 16X , 'Radi  us  to  inside  of  disk  in  infinite  plane  *  ', 

1  F9. 3, '  meters  ' ) 

804  FORMAT  ( 1 6X ,' Detector  is  1  meter  above  smooth  plane') 

620  FORMAT! 16X, 'Terrain  factor  ■  ' ,  F 1 2 . 4 ) 

STOP 

END 

C 

SUBROUTINE  PROTECT (DRF ) 

REAL  TOLER (10) , DMFP (13) , KERN, MU 

EXTERNAL  DISK, GROUND 

COMMON/ ENERGY/ EKE V , MU , BU , W , THETA , DOVER W 

COMMON/KERNAt /TOTKRN , TOTDOS , TDIRT , XX JM1 , D2 , R 1 

DATA  TOLER/ 1 . E-4 , i . E-4 , 1 , E-4 , 1 , E-4 , 1 . E-4 , 1 . E-4 , 

1  l.E-4, l.E-4, l.E-3, l.E-2/ 

DATA  DMFP/. 003 ,.01 ,.05, . 1 , . 2 , . 5 , 1 . , 2. , 5 . , 1 2. / 

C . RHO  *  DENSITY  OF  AIR  AT  STP  0.031293  gm/cm**3 

RHO  *  0.001293 
AMUEKV  ■  AMU(EKEV) 

MU  =  AMUEKV  *  RHO 

C . DISTANCES  ARE  IN  METERS,  MU  IN  INVERSE  CENTIMETERS,  XJ  IS  IN  MFPS 

TOTDOS  «  2. 

TOTKRN  *  0. 

TDIRT  «  0. 

N 1  *0 
N2  =  0 

XDOSE  =  R0ENT3 (EKEV) 

DO  805  J  «  1,10 
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XJ  «  DMFP(J) 

TOL  ■  TOLER ( J ) 

C  *  *  *  *  #  IF  R 1  IS  GREATER  THAN  XJ  DO  NOT  INTEGRATE  BACKWARDS 
XJM1  *  SI  *  MU  *  120. 

IF ( X J . LE. XJM1 )  GOTO  335 
N1  »  N1  +  1 
D2  =  X J /MU/ 103. 

IF(Nl.EQ.l)  THEN 
XXJM1  *  XJM1 
ELSE 

X JM1  ■  DMFPiJ-1) 

END  IF 

D  1  *  X JM 1  /MU/ 1 00. 

CALL  SIMPS (01 .02, TOL, KERN, DISK) 

DS  «  X DOSE  #  KERN 
TOTKRN  =  TOTKRN  +  KERN 
TOTDOS  =  TOTDOS  +  DS 

CALL  SIMPS (D 1 ,02, TOL, DIRT .GROUND) 

TDIRT  *  TDIRT  ♦  DIRT 
805  CONTINUE 

DRF  »  TDIRT/TOTKRN 
END 
C 

FUNCTION  GROUND ( R ) 

C********  POINT  KERNAL  FOR  A  DISK  SOURCE 

REAL  MU,  HYP,  MUCQNC,  MUDIRT,  MFPG ,  MFPA IR , MFP , I ND IR 
COMMON  /ENERGY/  EKEV,  MU,  BU,  W,  THETA, DOVERW 
COMMON  /FLAG/  RMAX 
C  NSKIP  -  0 

PI  »  3.  141592654 
HEIGHT  ■  1.0 

C . STANDARD  REFERENCE  POINT  FOR  A  DOSE  IS  1  METER  •  'OM  GROUND 

C .  FIND  THETA  FROM  R  AND  HEIGHT 

IF  <R. GT. 0, )  THEN 

BETA  »  AT  AN (HE  I GHT/R) 

ELSE 

BETA  *  PI/2. 

END  IF 

THETAR  »  THETA/180. «PI 
TAVE  *  0. 

A  »  W 

C . IF  BETA  GREATER  THAN  THETA  THERE  IS  NO  ROUGH  GROUND  EFFECT 

IF ( BETA. GT. THETAR)  GOTO  100 

C . FIND  A  FROM  W,  THETA,  AND  BETA 

A  »  W  #  SIN (BETA)  *  COS  (THETAR) /SIN ( THETAR  ♦  BETA) 

C .  A  IS  THAT  PART  OF  GROOVE  THAT  IS  NOT  IN  SHADDOW 

C .  FIND  AVERAGE  THICKNESS  OF  EARTH  INDIRECT  PART  PASSES  THROU 

TAVE  »  W/4,  #  SIN  (THETAR)  /  (COS  (THETAR) *SIN  <  BfETA)  ) 

C*#***  DO  NO i  WASTE  TIME  CALCULATING  UNNECESSARY  VALUES  OF  TAVE 

C .  FIND  D,  THAT  PART  OF  GROUND  THAT  IS  STILL  FLAT 

100  D  »  DOVERW  *  W 
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DIRECT  =  '  A  +  D)  /  (W  +  D) 

1  INDIR  =»  (W  -  A)  /  CM  +  D) 

C  WRITE (*,933) THETA, THETAR , DOVERW , W,TAV£, SETA, R, HEIGHT 
HYP  =  SQ.RT  (HEIGHT  *  H  £  1 3HT  +  R*R ) 

MFPAIR  =  MU  *  HYP  *  133. 

BUAIR  =  3UFIEKEV,  MFPAIR) 

RHQDRT  *1.3 

SPEC.  GRAY.  OF  DIRT  13  APPROX  1.3  SOURCE:  MECH.  ENG.  HANDBOOK 
MUD  I RT  =  RHQDRT  *  MUCONC ( EKEV > 

MFPG  =  MUD  I  RT  *  TAVE 
IF (MFPG. GT. 12)  THEN 
TERAIN  «  3. 

GOTO  20 

E  N  D I F 

3UF 1  *  3UF  < EKEV , MFPC  +  MFPAIR) 

BUF2  =  BUF  (EKEV, MFPG) 

DUF3  «  BCONC (EKEV , MFPG ) 

TERAIN  «  INDIR*  EXP(-MFPG)  *  BUF3  *  BUF1/DUF2 
23  D 1 5 K 1 =  .5  *  BUAIR*EXP (-MFPAIR)  *  R/HYP**2. 

GROUND  *  DISKI  *  (DIRECT  +  TERAIN/BUAIR) 

RETURN 

END 

C 

SUBROUTINE  SIMPS (LOWER, UPPER ,  TOL,  SUM ,  F) 

C  NUMERICAL  INTEGRATION  BY  SIMPSON'S  RULE. 

C 

INTEGER  IN,  OUT,  PIECES,  I,  P2 
REAL  X,  DELTA,  LOWER,  UPPER,  SUM,  TOL 
REAL  ENDSUM,  ODDSUM,  SUM1,  EV3UM 
COMMON  / INOUT/  IN,  OUT 
C 

PIECES  =  2 

DELTA  =  (UPPER  -  LOWER)  /  PIECES 
ODDSUM  »  F (LOWER  +  DELTA) 

EV3UM  *  3.0 

ENDSUM  =  F  (LOWER)  +  F(UPPER) 

SUM  =  (ENDSUM  +  4  *  ODDSUM)  *  DELTA  /  3.3 
3  PIECES  «  PIECES  *  2 

P2  =  PIECES  /  2 
SUM  1  =  SUM 

DELTA  »  (UPPER  -  LOWER)  /  PIECES 
EVSUM  »  EVSUM  +  ODDSUM 
ODDSUM  *  0.3 
DO  10  I  =  1,  P2 

X  =  LOWER  +  DELTA  *  (2  *  I  -  1) 

ODDSUM  =  ODDSUM  +  FIX) 

10  CGNTINUE 

SUM  =  (ENDSUM  +  4.3  *QDDSUM  +  2.0  *  EVSUM)  *  DELTA  /  3.0 

IF ( ABS ( SUM  -  SUM1 )  .GT.  ABSfTOL  *  SUM))  GOTO  3 

RETURN 


o  n 
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END 

C 

FUNCTION  DISK(R) 

C»*****«*  POINT  KERNAL  FOR  A  DISK  SOURCE 
REAL  MU,  HYP 

COMMON  /ENERGY/  EKEV ,  MU,  9U,  W,  THETA,  DQ''ERN 
A  =  1.0 

C . STANDARD  REFERENCE  POINT  FOR  A  DOSE  IS  1  METER  FROM  GROUND 

HYP  =  SORT ( A*A  +  R*R ) 

BU  =  BUF ( EKEV ,  MU  *  HYP  *  100. ) 

DISK  *  .5  *  BU*EXP(-MU*100.*HYP)*R/HYP**2. 

RETURN 

END 

C 

FUNCTION  AMU (EKEV) 

C  MASS  ATTENUATION  IN  cm*2/gr 

REAL  £  ( 26 ) ,  MU0RH0<26) ,  LOGINT 

C . AIR  MASS  ATTENUATION  DATA  FROM  RADIOLOGICAL  HEALTH  HANDBOOK 

-,  JAN  1970,  PG139 

DATA  E/10.  , IS. ,20. ,30. ,40. ,50. ,60. ,80.  ,  100.  ,  150. ,200. ,300. , 

1  400.  ,500. ,600. ,800. ,  1000.  ,  1500. ,2000.  ,3000.  ,4000.  ,5000.  , 

2  6000. ,8000. ,10000. ,15000./ 

DATA  MUORHO/  4 . 99 , 1 . 55 , 0. 752 , 0. 349, 0. 248 , 0.  208 , 0. 1 88 , 0.  1 67 , 0.  1 54 , 

1  0.136,0.123,0. 107,0.0954,0.0870,0.0805,0.0707,0.0636, 

2  0.0519,0.  0445,0.0358,0.  0308,0.0275,0.0252,0.0223, 

3  0.0204,0.0181/ 

C 

C - USE  LOG-LOG  INTERPOLATION 

LOG  I  NT  <  X , X 1 , X2 , Y l , Y2>  ■  E  X  P (  (AL0GIY2)  -ALOG ( Y 1 )  )  * 

1  (ALOG(X)-f.LOG(Xl)  )  /  (AL0G(X2)-AL0G(X1)  )  +  ALOG(YD) 

C 

C 

J  =  2 

C - use  lowest  two  data  points  for  energies  below  table 

IF (EKEV. LT. E  ( 1 ) )  GO  TO  212 
200  IF (EKEV  -  E  ( J ) )  212,211,210 

210  J  »  J  ♦  1 

I F ( J . LT. 26)  GOTO  200 

C - use  highest  two  data  points  for  energies  above  table 

GOTO  212 

211  AMU  ■  MUORHO  ( J ) 

RETURN 

212  CONTINUE 

AMU  »  LOGINT (EKEV,E(J-1)  ,E(J)  , MUORHO ( J-l ), MUORHO (J ) ) 

RETURN 

END 

C 

FUNCTION  ROENTG (EKEV) 

cmmmmnnmmnmxmmnnnmmmnmmnnmmnn 

C  This  function  converts  BUF  adjusted  flux  to  R*m**2/hr / C i 
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7.  7. 7. 7.  X  X  XXX  X  X  X  X  X  X  X  X  X  X  X  XXXX  X  X  X  7.  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  ■/.  X  X  X  X  X  X  X  X  X  X  7. 7.  X  X  X  X  X  X 
REAL  E  (26) ,  MU0RH0(26) ,  LOS  I MT 

—  A i r  mass  absorption  data  from  NSRDS-NBS  29,  1V69  pg  20  k  21 
DATA  E/1J. ,15. ,23. ,33. ,43. ,50. ,60. ,83. ,133. , 153. ,233. ,333. , 

1  433. ,533. ,603. ,383. , 1333. , 1533. ,2303. ,3330. ,4030. ,5330. , 

2  6333. ,8330. ,13333.  ,15333./ 

DATA  MUORHO/  4.61 , 1.23,3.511 ,3. 143,3.3669, .9436, .3335, .3243, 

1  .3234,. 3253,. 0263,. 3239,. 0295,. 0296,. 3295,. 0239,. 3278, 

2  . 0254, .3234, .0235, .0136, .8174, ,8164, ,3152, .8145, .0132/ 


--USE  LOG-LOG  INTERPOLATION 

LOG  I  NT ( X , X 1 , X2 , Y 1 , Y2)  =  EXP  (  ( ALOG  ( Y2 ) -ALOG  ( Y 1 ) ) * 

1  (ALOG  (X) -ALOG  (XI ) )/ (ALOG 1X2) -ALOG (X l ) )  +  ALOG  ( Y1 ) ) 


J  =  2 

C - u s s  lowest  two  data  points  tor  energies  below  table 

IF  < G K E V . LT. E ( 1 ) )  GO  TO  212 
233  I F ( E K E V  -  E  ( J ) )  212,211,218 

210  J  =  J  +  l 
!F(J.LT.26)  GOTO  233 

C - use  highest  two  data  points  tor  energies  above  table 

GOTO  21.7 

211  ROE  NT  (i  =  MUORHC(J)  *  EKE  V  *  244.3365  / 1033. 

RETURN 

2)2  CONTINUE 

ROENTS  =  LOGIN T(EKEV,E(J-1) , E ( J ) , MUORHO(J-l) , MUORHQ  ( J ) )  * 

1  EKEV  f  244.3365  /  1000. 

RETURN 

END 

C 

FUNCTION  BUF (EKEV, MUR) 

REAL  E  ( 25 )  ,  Al(25),  A2 ( 25 ) ,  ALPHA  1 ( 25 ) ,  ALPHA2  ( 25 )  ,  ALPHA3<25> 
REAL  MUR,  LOGINT,  SEMINT 
COMMON  /OLD J /  EQLD,  JOLD 

DATA  E/ 15. ,20. ,30. ,40. ,50. ,60. ,80. ,130. ,150. ,200. ,330.  , 

1  4  33.  ,530. ,630. ,300. ,  1033. ,  1530. ,  2300. ,3030. ,4003. ,5000.  , 

2  6303. ,  3003. ,  10000.  ,15300. / 

DATA  A1/1.258E1,4.963,1.039E1,1. 133E2.5. 106E2, 

1  1.641E3, 1. 477E3, 1.583E3, 1.242E3, 1.206E3, 

2  1.251E3,!. 182E3, 1 . 232E3, 4 . 3 1 6E3, 1. 102E3, 

3  1. 123E3,2. 941S2.4. 159E2, 1. 162E2, 1.923E1 , 

4  1. 251  El ,  1  .  0471.1  ,  1.01  1E1 ,8.839,6.661/ 

DATA  A2/-3. 984E-1 , -6, 395E-1 , -6. 924E-1 , -1 , 462E2 , -6. 189E2, 

1  -2. 712E3,-2.543E3,-2.736E3,-2.317E3,-2. 149E3, 

2  -1 . 756E3 , -2. 019E3,-1 . 664E3 , -4. 588E3 , -1 , 303E3 , 

3  - 1 . 1 74E3 , -2.  179E2,-2.687E2,-6.740E1 ,-1.417El , 

4  -6.071,-3.919,-3.046,-2.360,-1.496/ 

DATA  ALPHA1 /-2. 509E-2, -1 . 058E-3, -3. 1 74E-2 , -2. B52E-2 , -4 . 23 IE-2 , 

1  -4. 888E-2,-7.303E-2,-8. 1 90E-2 , -8. 536E-2 , -7. 780E-2 , 

2  -5. 541E-2.-3. 850E-2,-2. d43E- 2,-1 . 75lE-2,-l . 1 41E-2, 

3  - 1 . 03 1 E-2 , -4. 964E-2, -3.  784E-2, -2. 395E-2 , -2. 575E-2 , 
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4  -4.  129E-2.-5. 1 76E-2  , -4 . 734E-2 , -4. 339E-2  , -5 . 032E-2/ 

DATA  ALPHA2/3. 955E-I ,5. 535E-1 ,9. 233c- 1 ,6.613E-3,-1.852E-2, 

1  -2.589E-2.-4. 1 09E-2 , -4 . 753E-2 , -5. 023E-2 , -4 .  773E-2 , 

2  -3. 475E-2,-1.673E-2,-1.351E-2,-l. 4069-2,-2.  14 IE-3, 

3  -5.925E-3,-3.988E-2,-2.974E-2,-7.309E-3,2.034E-2, 

4  5. 330E-2, 7. 394E-2t  9. 32 1 E-2 , 9. 1 46E-2 , 1 . 01 3E- 1  / 

DAJA  ALPHA3/-2. 6S9E-2.3. 524E-2,6. 484E-2.6. 739E-2f6. 537E-2, 

1  6.546E-3,1.458E-3,-6.323E-3,-9.9?7E-3,-9.998E-3, 

2  1 . 52SE-2, 1 . 292E-2.2.  686E-2,3.682E-2, 4. 24E5E-2, 

3  7. 813E-2, -6. 392E-2, -4. 662E-2»-3. 262£-2,  -  4.  182E-2, 

4  -5.231E-2,-6.205E-2,-5.452E-2,-5.524E-2,-5.774E-2/ 

BUILD  UP  FACTOR  DATA  FROM  NUC  SCI  S<  END  78  PG74  1981 

BF ( A , B , C , D , E , UR ) =  A*EXP  ( -C*UR ) +B*EXP ( -D*UR )  +  ! 1 -A-B ) *EXP ( -E*UR ) 

C - USE  LOG-LOG  INTERPOLATION 

LOGINT ( X ,  X 1 , X2 , Y 1 , Y2 )  =  E XP ( ( ALOG ( Y2) -ALOG ( V 1 ) ) * 

1  (AL0G(X)-AL0G(X1))/(ALDG(X2)-AL0S(X1))  +  ALOG  C  Y 1 ) ) 

SEMINT ( X , X 1 , X2 , Y1 , Y2)  »  EXP ( (ALOG (Y2) -ALOG ( V 1 >) * 

1  (X-Xl )  /  (X2-X1 )  +  ALOG  (VP) 

DO  NOT  LOOK  UP  CONSTANTS  IF  YOU  ARE  NOT  CHANGING  ENERGY 
I F ( EKEV. EQ. EOLD) THEN 
J  a  JOLD 
GOTO  212 
END  I F 

EOLD  «  EKEV 
C 

J  =  2 

C - use  lowest  two  data  points  For  energies  below  table 

IFIEKEV.LT. E(l) )  GO  TO  212 
200  IF (EKEV  -  E  ( J )  )  212,211,210 
210  J  »  J  +  1 

IF ( J. LT, 23)  GOTO  200 


C - use  highest  two  data  points  For  energies  above  table 

GOTO  212 

211  IF (MUR. GT. 48. )  GOTO  213 

3UF  a  9F ( A 1 (J) ,A2(J> , ALPHA  1  (J) ,ALPHA2(J) , ALPHA3  <  J ) , MUR ) 

C  WRITE ( 1 ,900) EKEV, BUF , J , A1  (J)  , A2 (J ), ALPHA  1 ( J ) , ALPHA2 ( J ) , ALPHA3 ( J ) 


JOLD  »  j 
RETURN 

212  CONTINUE 

IF (MUR. GT. 40. )  GOTO  213 

X 2 3  BF ( A 1 ( J ) , A2 ( J ) , ALPHA  1 ( J ) , ALPHA2 ( J )  ,ALPHA3(J) , MUR ) 

X 1 =BF ( A l ( J- 1 ) , A2  !  J-l ) , ALPHA  1 C  J - 1 ) ,AL?HA2(J-1> , ALPHAS  <  J- 1 >  , MUR ) 
BUF  -  LOGINT (EKEV,E (J— I )  ,E(J),X1,X2) 

JOLD  ■  J 
RETURN 

213  CONTINUE 

C . IF  MFPS  ARE  GT  40  INTERPRET  FROM  END  OF  RANGE 

X 2 *  BF ( A 1 ( J ) , A2  <  J ) t  ALPHA  1 ( J ) , ALPHA2 ( J ) , ALPHA3 ( J  > ,40. ) 

X 1 a8F ( A 1 < J- 1 ) , A2  ( J- 1 ) , ALPHA  1 ( J- 1 ) , ALPHA2 ( J- 1 ) , ALPHA3 ( J- 1 ) ,40. ) 
BUF2  »  LOGINT ( EKEV , E  <  J - 1 )  ,E(J) ,X1,X2) 

X2»  BF ( A1 (J) ,A2(J) , ALPHA  1  (J) , ALPHA2 ( J )  ,  ALPHAS  C J ) ,35. ) 
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X  1  =  BF ( A 1 (J-l ) ,A2(J-11 , ALPHA  1 (J-l) , ALPHA2 ( J - 1 ) , ALPHAS (J-l) ,33. ) 

BUF 1  =  LOGINT  (EKEV,E(J-1) ,ECJ) , X 1 , X  2 ) 

5UF  =  SEMI  NT ( MUR , 35. ,43. , SUFI , BUF2) 

J  OLD  =  J 
RETURN 
END 
C 

FUNCTION  BCONC  (EKEV , RMFP) 

C . THIS  FUNCTION  CALCULATES  THE  BUILD  UP  FACTOR  IN  CONCRETE  OR 

C . CONCRETE  LIKE  MATERIALS  FROM  THE  NBS  SOLUTION  USING  12  PARAMETERS 

C - BY  E I SENHAUER  AND  SIMMONS  REF  NUC  SCI  AND  ENS  56,263-270,1973 

C 

C . R  IS  IN  MFPS 

r 

REAL  £  (25)  ,  BO  ( 25) , MUC (257  ,L05 JNT , MUCONC 
REAL  A *3  (25)  ,  A 1  (25  /  ,92(235  ,  A 3  (25)  ,  A 4  (25) 

REAL  32(25) ,81  (25) ,32  (25) ,33(25) ,94(25) 

REALMS  BFC 

COMMON  / 0 L D J 2 /  ECLD2 ,  J0LD2 

DATA  E/  1  3. ,  23. : 30. ,  40. ,50. ,63. ,30. , 103. , 150. ,  223. ,300.  , 

1  4  33.  ,  533. ,603. ,823. ,  1332. ,1503.  ,2330. ,3330. ,  4003. ,3330.  , 

2  6333. ,3300. ,10023. ,15223./ 

DATA  MUC/3. 3123,3. 4453,1. 1133, 0.5539,3. 3633, 7.2734, 3. 2304,0. 1734, 

1  0.1399,0.1250,0. 1073,3.0933,0.0873,0.3807.2.3739, 

2  3.3637,3.3519,0  0443,3.0365,0.3319,3.3290,0.0270,0.3245, 

3  3.3231,3.3215/ 

DATA  DO/ 3. 1539,3. 1437,3. 1327,0. 1233,3. 1123, 

1  0.1118,3.1215,0.1376,0. 1791,0.2149,0.2693, 

2  3.3333,0.3409,3.3674,0.4385,0.4404,0.4961 , 

3  3. 5233, 3. 3t 79, 3. 5933, 0.6 153, 3.6241  ,0.6425,0.6560,0.6676/ 

DATA  A3/i.  11 35E-2, 2. 4283E-2, 0,36645, 3.  1  1734,3. 1591  3, 

1  3. 18644,0. 22372,0.24651 ,0.27760,0.29820,3.31937, 

2  3,32929,3. 33559 , 0.  33965 , 3. 344 1 3 , 3.  34673 , 0. 35 1  1 3 , 

3  3.35445,0.35691,3.35425,3.34739,0.33878,3.31686, 

4  3.29369,3.24338/ 

DATA  A 1  /-0.  03 1  1 583, -0.332461 3, -0.3060622, -0.  3092903, 

1  -0. 31 5356, +3. 025044,-0. 010144,3.043674,0.021783, 

2  3.029393,0.312225,0.0395543,0.0044976,0.0373127, 

3  0.0391422,3.0098983,0.010016,0.0099221,0.0078632, 

4  3.310371,0.343629,0.330309,-0.013397,-0.016613,-0.0029303/ 

DATA  A2/5.493E-04,1.2409E-03,3.9389E-03,8.6343E-03, 

1  1. 5 169E-02,1.1833E-02 ,3.79328-02, 3. 2059E-02. 2. 9878E-02, 

2  4 . 8446E-02  , 

3  5.3031E-02,5.0474E-02,2.2018E-02,2.7012t-02,2.2857E-02, 

4  2.0376E-02,  1.9785E-02,1.9236E-02,2.  1 80 1 E-02 , 3. 9089E-02 , 

5  -9. 3940E-03 , -8. 8533E-03 , 2. 883BE-02, 2.  7261E-02 , -1 , 91 01E-02/ 

DATA  A3/2. 1 S03E-04 , 5. 0984E-04 , 1 . 7224E-03 , 4 . 0328E-03 , 

1  7.5754E-03, 1 . 3 1 53E-03 , 3. 0337E-02 , 2. 6207E-02 , 3. 19I5E-02, 

2  8.2357E-02,  1 . 0486E-0 1  ,  1 . 1029E-01  , 3. 4039E-02 , 3 . 1744E-02, 

3  4.4123E-02,5.0435E-02,5.2275E-02,4.9103E-02,3.900b'E-02, 

4  2. 9853E-02 , 2. 88B8E-02 , 2. 5095E-02 1 1 . 8330E-02 , 1 , 4079E-02 , 
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3  2.281 7E-02/ 

DATA  A4/4.2382E-05, 1 . 0829E-04 , 3. 9905E-04 , 1 . 0308E-03 , 

1  2. 2S27E-33, 0.02003000, 1 . 2457E-02 , 3. 1 5 1 4E-36 , 5 . 804 1 E-02  , 

2  — 1  .  4633E  —  04  ,—8. 641  IE  —  34  ,  — 6. 3551E-38, 1 . 3539E-31  ,  9.  4332E-02, 

3  7.3339E-02,5.9493E-02,4.071  4E-02,3.2172E-02)2.  2137E-02, 

4  2.3373E-33,6.7343£-35,3.5731E-03,1.5436E-03,3,6397E-07, 

5  9 . 3532E-03/ 

DATA  90/0.57471,0.56334,0. 54853,0.53350,0.53737, 

1  0.54733,0.37172,0.60362,0.63913,0.69464,0.  73934 , 

2  3.  76  455,0.78166,0. 79363,0.82851,0.31671,0.32277, 

3  0.31861,0.80467,0.79201,0.78176,0.77269,0.73988, 

4  0.75116,0.73974/ 

DATA  32/0.72706,0.71518,0.69814,0.68724,0.68395, 

1  0.90637,0.69136,0.86298,0.77287,0.84634,0.82393, 

2  0.32999,3.76467,3.79645,3.32369,3.30432,3.79932, 

3  3.79546,3.81733,0.85749,0.72380,0.63341,3.85214, 

4  0.86119,3.55119/ 

DATA  B 1  /0. 42393. 3.4 1368, 0.39364, 0.3631 0,3. 23648, 

1  3. 70718,0.30330,0. 69185,0.68532,0. 71979,0.68837, 

2  3.65925,3.56762,3. 53349,0.57175,2.55269,0.52363, 

3  3.53853,3.46877,0.30734,0.83930,0.83961,8.57962, 

4  3.56135,1.85134/ 

DATA  B3/ 3. 37317,0. 86506,0.85649,0.85396,0. 85703, 

1  1 .  20436,0.85532,0.98931,0.87818,1.31431, 1.0:389, 

2  1.  00623,0.87335,0.90888,0.91929,3.91359,0.91265, 

3  0.90861,0.91132,3.96145,3.96383,0.94894,0.95976, 

4  0.97377,0.36292/ 

DATA  B4/0. 96440,0.96154,0.95995,0.961  10,0.  96562, 

1  1.00030,0.98034,1. 17560, 1.00645,1.09682,1.07209, 

2  1.34178, 1.30139,0.99831 ,0.99455,0.99168,0.98767, 

3  0.99542,0.98406,1.01558,1.10890,1.30529,1.01101 , 

4  1.22902,0.97741/ 

3FC(R,U,A0,A1,A2,A3,A4,B0,B1,B2,B3,34,DO)  ■  1  +  < 

1  A0*U*R#DEXP (-DBIE < J#R/B0> ) +A 1 ♦ ( U*R/31 ) #*2#DEXP ( -DBLE l U*R/B  1 )  )  + 

2  A2  «  (U#R/B2)##2  *  DEXP  (-D9LE  (U*R/32) )  + 

3  A3  *  (U»R/B3>#*2  *  DEXP  ( -DBLE (U*R/33> )  + 

4  A4»  (U#R/84)*#2#DEXP (-DBLE  (U«R/84) ) )/ (DO* DEXP (-DBLE (U*R)  ) ) 

C - USE  LOG-LOG  INTERPOLATION 

LQQINT  <X,X1,X2,Y1,V2)  *  EXP  (  (ALOG  ( Y2I-ALQG <  Y 1 ) ) * 

1  (ALOG(X)  -  ALOG4X1 )  >  /  (ALGG(X2)— ALOQ(Xl) )  +  ALOG  (YD) 

R  ■  RI4FP  /  MUCONC  ( EKE V ) 

C  WRITE!#, 1000)  ' BCONC  2 ' , EKEV 

:*###**  DO  NOT  LOOK  UP  CONSTANTS  IF  YOU  ARE  NOT  CHANGING  ENERGY 
IF (EKEV. EQ.E0LD2) THEN 
J  »  JOLD2 
GOTO  212 
END  I F 

E0LD2  *  EKEV 
C 

J  =  2 

C - us*  !ow**t  two  data  points  for  energies  b*low  table 
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IF  -V.LT.E(l))  GO  TO  212 
233  IF  t  .■  -  E  ( J  ) )  212,211,213 

213  J  =  -  +  1 

I F ( J . L T . 2 5 »  GOTO  223 

C - l; '3  highest  two  data  points  tor  energies  above  table 

GOTO  212 

711  BC0NC=3FC(R,1UC!J) , A3 ( J ) . A  1 (J) , A2 ( J ) , A3  <  J )  ,A4 (J) ,93(J)  ,  B 1  (J)  , 

1  a  2 ( J ) , S3 ( J ) , B  4  C  J ) ,DD(J> ) 

URCONC  =  R  *  MUC(J) 

J0LD2  =  J 
RETURN 

212  CONTINUE 

X2  =  BFC(R,MUC(J) , A  3 ( J ) , A 1 (J) , A  2 ( J ) ,  A3  ( J ) , A  4 ( J  >  ,  B3  ( J )  ,51  (J)  ,  B  2  ( J  >  , 

1  S3 ( J ) • 34 ( J )  ,  DO  ( J )  ) 

K  =  J  -  1 

X  1  =  3  F  2  ( R  ,  M  U  C  ( K )  ,  A  3  ( K )  ,  A 1  ( *  >  .  A  2  ( K )  ,  A  3  ( K  >  ,  A  4  <  K )  ,  9  3  ( K )  ,81  (K)  ,  B  2  ( K )  , 
1  3  3  ( :< )  ,94(K)  ,  DO  (K)  ) 

BC3NC  *  LOGIN!  (EKEV,£(J-1) ,  E  ( J  ) ,X1 ,  X2) 

L'RCGNC  =  R  *  LOG  I  NT  (EKEV,E(J-1)  ,E(J)  ,MUC  (J-l  >  ,  MUC  ( J )  ) 

J0LD2  *  J 
RETURN 
END 
C 

REAL  FUNCTION  HUCONC <EKEV) 

C - this  fucntion  calculates  mu  over  rho  for  concrete  in  cm2/gr 

REAL  £(25)  ,  fl'JORHO (2*3 )  ,  lOGINT 

C - russ  attenuation  data  from  Nuc  Scl  &  F.ng  *6  pr  267,  1975 

C 

DATA  E/ 15. ,23. ,33. ,43. ,53. ,  63. ,83. ,  130. , 133. ,230.  ,  330.  , 

1  433. ,533. ,633. ,333. , 1333. ,1533. ,  2333. ,3333. ,4333. ,5333.  , 

2  6333. ,9333.  ,13333. ,  15333.  / 

DATA  KU0RH0/9. 0130,3, 4453, 1. 1193,0.5599,3.3638,3.2734,3.2304, 

1  3. 1734,3.1399,3.1253,3.1373,3.3933,3.3373,3.3337,0.0739, 

2  3. 36  37, 3. 3519, 3. 344  9, 3. 3365, 3. 33 1  9, 0.3293, 3. 3273, 3. 0245, 

3  3.3231,3.3215/ 
r 

C - USE  LOG-LOG  INTERPOLATION 

L0GINT(X,X1,X2,Y1,Y2)  *  EXP  (  ( ALOG  ( Y2 ! -ALOG  ( Y  1 )  )  «■ 

1  (AL03(X)-AL03(X1))/(AL03(X2)-AL0G<X1)>  +  ALOG  <  Y 1 ) ) 

C 

C 

J  =  2 

C-- - use  lowest  two  data  points  fo r  energies  below  table 

IF(EKEV.LT.Ed) )  GO  TO  212 
233  IF ( EKEV  -  E<J) )  212,211,213 

210  J  ■  J  ♦  l 

IF ( J. LT. 26)  GOTO  230 

C - use  highest  two  data  points  for  energies  above  table 

GOTO  212 

211  MUCONC  *  MUCRHO(J) 

RETURN 


212  CONTINUE 

MUCONC  *  LOGINT  (£KEV,E(J-1) ,E(J) - MUOPHO ( J- 1  ) , MUORHO ( J )  ) 

RETURN 

END 
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FALLOUT  PROGRAM  DESCRIPTION  AND  LISTING 
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Tha  Fallout  program  calculates  the  DRFs  -for  a  -fallout 
source  whan  tha  detector  is  cantered  insida  a  structure. 

Tha  prcgr am  modal  is  discussed  in  Section  2.7.  Required 
input  includas  tha  equivalent  height  and  radius  of  the 
structure,  the  mass  thickness  of  tha  roof  and  walls  and  the 
three  ground  roughness  factors  w,  d/w  and  ¥  . 

When  the  Fallout  program  is  run,  it  will  ask  for  input 
data  as  shown  in  Figure  27.  Tha  first  input  required  is  a 
title  and  data  for  the  output  listing.  The  program  then 
asks  for  tha  equivalent  radius  and  height  of  the  structure, 
fallowed  by  the  mass  thicknesses  of  the  wall  and  roof.  The 
program  then  prompts  the  user  with  a  few  suggested  values 
for  the  three  ground  roughness  factors  (see  Section  2.6.2). 
Other  values  for  the  ground  roughness  factors  may  be  used. 
The  program  then  calculates  the  resultant  DRF's  for  the 
1-Hr.  bomb  fallout  spectra,  the  PWR-2  (RSS)  fallout  spectra, 
the  TMI  fallout  spectra  and  the  SL-i  fallout  spectra 
followed  by  the  DRFs  for  energies  from  15  KaV  to  15  MeV.  A 
typical  output  listing  in  shown  in  Table  14. 
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B  >4  »1 1  out 

INPUT  TITlE/DATE  OF  OUTPUT 
S*«ll  Nood  Houtt •  16  Stpt  1994 

INPUT  ESUIVAlENT  RADIUS  OF  STRUCTURE  IN  METERS  J. 
INPUT  EQUIVALENT  HEI8HT  OF  STRUCTURE  IN  METERS  5. 
NALL  MASS  THICKNESS  IN  gr/c»»»2  IB. 

ROOF  MASS  THICKNESS  IN  gr/cA*»2  IB. 

Thr»i  ground  Ficton  iuit  0*  input 
N  -  Th#  ch«r4ctori»f ic  trough  mdth 
0/N  -  Tn»  rttio  o<  41*t  »ur<*c»  to  trough 
*na  RSI  -  Th»  chArpctoriltlC  trough  *ngl«. 


Th»  Follomng 

*r»  luggtitid 

vllud: 

SURFACE 

N 

D/N 

Stootn  p 1 4n t 

B. 

1 . 

PivfP  Ar V4 

B. 3941 

1 . 

Lt«n 

2.3693 

B.  959 

Sr»v»lltd  Artt 

7.4678 

B.897 

Ploutd  Fi»ld 

25.9207 

B.939 

P51  •  45.  dtgr«ii 


INPUT  PS  I  iPigrtt •)  45. 
INPUT  N  2.1690 
INPUT  D/N  .959 


FIGURE  27  l¥Ei£«l_£gmey£g!:_£RT_Pi32i£¥_gyring  Iniiiatisn 

E*llsyt_PrSflC*mA 


Tho  program  also  calculates  the  unprotected  exposure  in 
R  pjr  Curie  of  fallout  par  square  meter.  The  exposures  from 
the  roof  source  and  ground  source  are  calculated  and  summed 
to  give  the  total  exposure.  The  fraction  of  the  total 
exposure  from  the  roof  source  is  given  as  the  roof 
contribution  (ROOF  CDNT. ) .  This  is  followed  by  the  DRFs  and 
the  ratio  of  the  DRFs  to  the  DRF  for  the  1— Hr.  bomb  fallout 
spectra.  The  last  column  of  data  is  the  exposure  one  would 
receive  from  an  interior  fallout  source  of  1  Curie  per 
square  meter.  The  program  assumes  one  gamma  per 
disintegration  for  the  1  Hr.  bomb  fallout,  the  PWR-2  fallout 
(R3S-Fal 1  out )  and  the  individual  gamma  energy  calculations,. 

The  program  is  written  in  Microsoft  FORTRAN  3.13.  This 
is  a  version  of  FORTRAN  77.  There  are  no  non-standard 
FORTRAN  commands  in  the  program.  The  program  was  run  on  a 
Il22i35_l0Sicy!3?20i§_EC2f222.i.°03l_£2!DByt3r  (PC)  which  uses  an 
§083  processor  and  an  §037  co-processor .  The  program  was 
run  under  the  MS-DOS  2.11  operating  system.  How*va:~ ,  it 
should  run  on  any  computer  capable  of  using  FORTRAN  77.  The 
program  takes  approximately  45  minutes  running  time. 

Caution,  the  program  is  written  is  FORTRAN,  the  input 
reguires  decimal  points.  The  program  may  not  converge  to  an 
answer  for  low  energies  (15  and  20  KeV)  for  wall  thickness 
greater  than  about  20  gr/cm3.  This  is  because  the  exposures 
are  less  than  B.43E-37,  which  is  the  limit  of  the  single 
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precision  data  type  used  -for  real  variables. 

The  MAIN  program  directs  the  input  and  output  of  the 
program  data.  Subroutine  SIMPS  carries  out  the  integrations 
by  using  Simpson's  rule.  The  number  o-f  increments  is 
automatically  increased  until  the  desired  accuracy  is 
obtained.  Subroutine  PROTECT  speeds  the  integration  o-f  the 
disk  source  by  dividing  the  disk  into  rings  which  converge 
to  integrals  -faster.  Function  DISKIN  -finds  the  kernel  for 
the  hypothetical  infinite  smooth  disk  source.  Function 
GROUND  calculates  the  kernel  for  the  fallout  using  the 
methods  developed  in  Section  2.7.  Function  DISKRF 
calculated  the  kernel  for  the  roof  source. 

•unction  AMU  finds  the  mass  attenuation  coefficient  for 
air.  Function  BUF  finds  the  buildup  factor  for  air. 

Function  MUCONC  finds  the  mass  attenuation  coefficient  for 
the  ground,  roof  and  walls.  Function  BCONC  finds  the 
buildup  factor  for  the  ground,  roof  and  the  walls.  Function 
RGENTG  finds  the  factor  for  converting  flux  to  exposure. 

The  listing  of  the  Fallout  program  follows. 


137 


PROGRAM  HOUSE 
CHARACTER  CT I TLE*9 1 
INTEGER  IN,  OUT 

REAL  E(23>  ,  MFP1,  LININT,  INSIDE 
REAL  EGROUP (23,4) ,  FRACT (23,4) 

EXTERNAL  CISKIN,  GROUND,  DISKRF 
COMMON  /ENERGY/  EKEV<W,  PSI,  DOVERW 
COMMON  ./SIZE/  RTHICK,WTHICK,8HGHT,BRAD 
COMMON  / INOUT/  IN,  OUT 
COMMON  / OLD J /  EOLD,  JOLD 
COMMON  /0LDJ2/  E0LD2 , J0LD2 

CA,'AAAA  NQTE  RADII  AR£  IN  m£tERS,  ENERGIES  IN  KeV,  DENSITIES  IN  gr/cm*2 
C 

DATA  E/ 1 5. ,23. ,30. ,40. ,53. ,60. ,83.  ,130.  ,  i 53.  ,230.  ,300.  , 

1  433. ,503. ,633. ,662. ,333. , 1009. , 1173. , 1332. , 1500. , 

2  2333. ,3233. ,4033. ,5033. ,6003. ,8033. , 1 3300 .,  1 5000 .  / 

DATA  EGROUP  /  25 . , 75. , ! 50 . , 253 . , 350 . , 530 . , 700. , 

1  903.  ,1  165.  ,  1523. , 1833. ,2253. ,2753. ,3503.  ,  4530.  ,5*3.  , 

2  50.  ,303.  ,753. , 1503. ,2500. . 15*0.  , 

3  723. ,637. ,364. ,234. ,83.  ,15*3.  , 

4  81.  ,234. ,364. ,637. ,662. ,  723.  ,  14*3./ 

DATA  FRACT  / . 327 1  , .  0 1 37 , . 0737 , . 04  76 , .  0929 , .  1 373 , . 1  7 1  7 , . 1 627  , 

1  .0839,. 0957,. 0299,.  0397,. 01 48,. 3042,. 0001  ,5*0.  , 

2  .  1696,.  243  1,. 27  63,.  2 133,. 0933,  15*0. , 

3  .  316, . 363, . 82, . 354 , .  326 , 1  5*3 .  , 

4  .0245,. 0529,. 7736,. 0642,. 0401,. 0151, 14*3./ 

C*******  LINEAR  INTERPOLATION  FUNCTION 

LININT(X,X1,X2,V1,Y2)  3  ( Y2-Y1 ) *  ( X-X 1 ) / ( X2-X  1 )  ♦  Y1 

Qaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa^aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 

CAATTEMPT  TO  CONTROL  EPSON  FX-30  PRINTER  FROM  FORTRAN  TO  COMPRESS  PRINT 
OPEN  ( 1 ,  FILE5  '  LPT  1  '  .STATUS3' NEW'  .ACCESS5* '  DIRECT 
1  , FORM3 'UNFORMATTED  '  ,  RECL=  1 ) 

WRITE  (1)27,15 
CLOSE  (1) 

OPENX 1  ,FILE  =  'LPT1  ') 


'AAAAAAAAAAAAAAA/ 


NAAAAAAAAA/ 


VAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 


PI  3  3. 141592654 
OUT  3  1 
EOLD  3  0.0 
E0lD2  3  0.0 
IN  3  0 
NSPEC  3  0 


cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C X7.XXX  Routines  For  inputinq  end  printing  given  data. 

c  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

WRITE!*, 293 ) 

290  FORHATdX,  '  INPUT  TITLE/BATE  OF  OUTPUT') 

READ ( * , 29 1 ) CT I TLE 
WRITE (OUT , 807 ) 

WRITE (OUT ,292) CTITLE 
WRITE (OUT ,808) 


i 


139 


J 


t 

A 


j 


1 


« 


t 


WRITE ( OUT , 304 ) 

333  WRITE!*, 331) 

331  FORMATUX,  '  INPUT  EQUIVALENT  RADIUS  OF  STRUCTURE  IN  METERS  '  ,  \  > 

332  READ  (*,333)  BRAD 

333  FORMAT (F13. 6) 

WRITE  t*,3ll ) 

311  FORMATUX, 'INPUT  EQUIVALENT  HEIGHT  OF  STRUCTURE  IN  METERS  '  ,  \ ) 
READ  (*,303)  3HGHT 
WRITE  <* , 321 ) 

321  FORMAT  (IX,  'WALL  MASS  THICKNESS  IN  gr/c,***2  '  ,  \ ) 

READ  (*,303)  WTHICK 
WR I TE ( * , 33 1 ) 

331  FORMAT  (IX,  'ROOF  MASS  THICKNESS  IN  gr/cm**2  ',\> 

READ  (*.303)  RTHICK 


WRITE  ( * , 600 ) 

WRITE!*, 610) 

READ (* ,303)  PSID 
psi  *  PSiD/iaa. *pi 
WRITE (*,611) 

READ ( * , 303 )  W 
WRITE!*, 612) 

READ ( * , 333)  DOVERW 

600  FORMAT ( 1 X ,' Three  ground  Factor*  must  be  input',/, 

1  ‘  W  -  The  characteristic  trough  width',/, 

2  '  D/W  -  The  ratio  oF  Flat  surFace  to  trough',/, 

3  '  and  PSI  -  The  characteristic  trough  angle.',/,/, 

4  '  The  Following  are  suggested  values:’,/, 

5  '  SURFACE  W 

6  '  Smooth  plane  •  0. 

7  '  Paved  Area  0.3941 

8  '  Lawn  2.3698 

9  '  Gravelled  Area  7,4678 

A  '  Plowed  Field  25.8237 

B  '  PSI  a  45.  degrees  ',/,/) 

610  FORMATUX, 'INPUT  PSI  (degrees) 

611  FORMATUX,  'INPUT  W  ' ,  \ ) 

612  FORMATUX,  'INPUT  D/W  '  ,\> 

WRITE!  CUT , 304 ) BRAD 
WRITE (OUT, 314) BHGHT 
WRITE (OUT ,324) WTHICK 
WRITE (0UTf334) RTHICK 
WRITE! OUT , 621 ) 

WRITE (OUT, 622)  PSID 
WRITE! OUT ,623)  W 
WRITE (OUT , 624 )  DOVERW 
WRITE (OUT ,803) 

C****************************************************************** 


D/W',/,/, 

1. ',/, 

1.  Vi 

0.959' ,/, 

0. 897  , / , 
0.339' ,/,/, 

i') 


C*************** INPUT  TITLE  FOR  LISTING*************************** 

807  FORMAT(1H0,25X, 'HOUSE  WITH  SEMI-INFINITE  DISK  SOURCE  IN  AIR', 
1  '  WITH  GROUND  FACTORS.  24  JULY  1984') 


I 


s 
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l 


] 


5 


1 


j 


9 


I 


1 


C*»H*T*l**HH**H*H*HH»n*H*»HH*H**HH*«Ht*HH#Hn*IH 

3  32  F0RilAT(lX,25X,F7. 1 ,  5X ,  7  (69. 4, 2X)  .4X.G9.4) 

333  FORMAT ( 1H3, 25 X ,  'E(KeV) '  , 3 X , 

1  '  UNPROTECTED' ,3X, 'ROOF' ,6X, 

2  ' GROUND ',  3X ,'  TOTAL  EXP.  ’  ,2X , 'ROOF  CQNT. ' , 3 X ,  ' DR F  ' , 3 X , 

3  'RATIO  TO  1.12  Hr.  ',2 X,  'INSIDE') 

233  FORMAT (1H3 , 25X ,  'Integration  to  12  MFPs  in  air.') 

334  FORMAT  < 1 X  f  2  5  X ,  'Equivalent  radius  of  structure  =  ', 

1  F9 . 3 ,  '  meters') 

314  FORMAT  (1X,25X,  'Equivalent  height  of  structure  - 
1  F9. 3  ,  '  meters  '  ) 

324  FORMAT (1X,25X,  'Wall  mass  thickness  =  ', 

1  F9.3,  '  g r / c m 2  '  > 

334  FORMAT (1 X ,  25X , 'Rcof  mass  thickness  = 

1  F  9 . 3 ,  '  g  r  /  c  m  2  ' ) 

834  FORMAT ( 1 X , 25X ,' Detector  is  1  meter  above  smooth  plane') 

621  FORMAT ( i X , 25X ,' The  Sround  Factors  are:') 

622  FORMAT  < 1 X , 3  3  X ,  'The  trough  angle  PSI  =  *  tF12. 4,  '  degrees') 

623  FORMAT (1 X  ,  32X The  charactori stic  trough  width  W  ■  '  ,  F 1 2 . 4 ) 

624  FORMAT (1X,30X,  'Th®  ratio  of  flat  to  trough  D/W  =  '  , F 1 2 . 4 ) 

535  FORMAT  <1X,25X,'1.12  Hr,  '  , IX , 7 (G9. 4 , 2X ) , 4X , G9. 4) 

513  FORMAT ( 1 X , 2 5 X ,  ' RSS-Fal 1  out ' , 1 X , 7 ( G9 .  4 , 2X ) , 4 X , 09 . 4 ) 

506  FORMAT ( 1 X , 25X ,  '  TM I -Fal 1  out ' , 1 X , 7 ( G9 .  4 , 2X ) ,4X,G9.4) 

507  FORMAT  X 1 X ,  25 X ,  '  SL1 -Fal 1  out ' , 1 X , 7 (G9. 4 , 2X ) ,4X,G9.4> 

291  FORMAT  ( A9 1 ) 

292  FORMAT ( 1 H0 , 25X , A9 1 ) 

0333083 0090330033003303303320333338333333309339330330333383930338008338 
033333  ROUTINE  TO  FIND  DRFS  FOR  SPECTRA 
DO  500  NSPEC  ■  1,4, +  1 
SUMGK  *  0. 

3UMGE  =  0. 

SUMUK  =  0. 

3UMUE  =  0. 

SUMRK  =  0. 

SUMRE  «  0. 

TEX  IN  »  0. 

DO  531  I  EKE V  »  1,23 

EKE V  ■  E6R0UP(IEKEV, NSPEC) 

FRCT  =  FRACTIIEKEV,  NSPEC) 

IF (FRCT. EQ. 3)  GOTO  531 
WRITE!*, 806)  EKEV 

C . LET  YOU  KNOW  PROGRAM  IS  RUNNING 

EXPOSE  *  RQENTG (EKEV) 

CALL  PROTECT ( 0. , DISKIN , UNPKRN) 

CALL  PROTECT  X BRAD , GROUND, PKRN) 

CALL  SIMPS (0. , BRAD, l.E-5, PROOF, DISKRF) 

CALL  S IMPS ( 0  , BRAD, l.E-5, INSIDE, DISKIN) 

TEXIN  ■  TEXIN  +  INSIDE  *  EXPOSE  *  FRCT 

SUMUK  «  SUMUK  +  UNPKRN  *  FRCT 

SUMUE  ■  SUMUE  ♦  UNPKRN  *  EXPOSE  »FRCT 


I 
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SUJ1RE  *  SLIMRE  +  PROOF  *  EXPOSE  *  FRCT 
501  CONTINUE 


5 


] 


4 


* 


i 


) 


i 


I 


TQTKRN  «  SUMGK  +  SUMRK 


TOTEXP  «  SUMGE  ♦  SUMRE 
DRF  «  TCTEXP/SUMUE 
I F  ( NSPEC . EQ .  1)  DRFBMB  ■  DRF 
RATIO  »  DRF/DRF3M8 
RCONT  »  SUMRE/TOTEXP 

IF (NSPEC. EQ. 1) WRITE (OUT ,505) SUMUE, SUMRE,  SUMGE, TOTEXP, 
1  RCONT, DRF, RAT  10 , TEX  I N 

IF (NSPEC. £0.2) WRITE (OUT, 51 0) SUMUE, SUMRE, SUMGE, TOTEXP, 
1  RCONT, DRF, RATIO, TEXIN 

IF (NSPEC. EQ. 3) WRITE (OUT, 506) SUMUE, SUMRE,  SUMGE , TOTEXP , 
1  RCONT, DRF, RATIO, TEXIN 

IF (NSPEC. EQ. 4) WRITE (OUT,  507) SUMUE, SUMRE, SUMGE, TOTEXP, 
1  RCONT. DRF, RATIO, TEXIN 
500  CONTINUE 


c  xnxnnmnmnmumrmnnmmmmnmmmmmnnn 

CmU  ROUTINE  FOR  GENERATING  ENERGY/DRF  TABLE 

cnnmnmmmuumxnmnmmmnmnnmmnnnmm 

DO  801  IEKEV  a  28, 1,-1 
EKEV  =  E  ( IEKEV) 


806 


WRITE!*, 306)  EKEV 

FORMAT (IX,  '  E  ( K  a  V )  *  ',510. 5, \) 

■LET  YOU  KNOW  PROGRAM  IS  RUNNING 
EXPOSE  »  RQENTG (EKEV) 

CALL  PROTECT (0. , DISK  IN , UNPKR) 

PROTECT  (BRAD,SRGUND,PRK5 
SIMPS (0. ,8RAD,l.E-5, PROOF, DISKRF) 
SIMPS (0. , 3RAD , 1. £-5, INSIDE, DISK  IN) 
=  INSIDE  *  EXPOSE 
«  UNPKR  *  EXPOSE 
PRK  *  EXPOSE 


CALL 
CALL 
CALL 
EXIN 
EXPUNP 
EXPPR« 

EXPRF  ■  PROOF  *  EXPOSE 
TOTEXP  »  EXPPR  ♦  EXPRF 

IF (TOTEXP.LT. l.E-33)  THEN 
RCONT  ■  .  8888888E33 

ELSE 


RCONT  ■  EXPRF/TOTEXP 

END  I F 

IF (EXPUNP. LT. 1 . E-30)  THEN 
DRF  -  . 0888888E33 

ELSE 


f 


DRF  «  TOTEXP/EXPUNP 

ENDIF 

RATIO  «  DRF/DRFBMB 

WRITE (OUT, 802) EKEV, EXPUNP, EXPRF, EXPPR, TOTEXP, RCONT, DRF, RAT  10, EX  IN 
801  CONTINUE 


1 
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233  CONTINUE 

C  J )  5  ]  3  3 3  3 3 :» 3 3 3 3  ? 3 3 3 3 3 3 i 2 3 3 f  3 3 9 3 3 § 3 3  3  j  3 3 -3 § 3 3 3 3 3 >3 2 § 3 3 3 3 9 g g 3 3 1 gggagggggg 3 3 
STOP 
END 

C  “  3  U  b 

SUBROUTINE  PROTECT  ( R 1  ,  D I SK  ,  TOTKRN! 

PEAL  TOLER (13) ,DMFP(13) , KERN , MFP 
C  EXTERNAL  DISK , GROUND , DI SKIN 

CCMM0N/ENER3Y/EKEV ,  W ,  P S I , DOVERW 

COMMON  /SIZE/  RTHIC'<(WTHICK,SHGHT,SRAD 

DATA  TQLER/l.E-4, 1 .E-4, 1 . E-4 , 1 . E-4 , 1 . E-4 , 1 . 2 -4 , 

1  l.E— 4, 1,£  —  4.1.E-3, l.E-J/ 

DATA  DMFP/. 333, .31, .05,, l,. 2, .5,1.  ,2.  ,5.  ,12./ 

C . RHO  =  DENSITY  OF  AIR  AT  5TP  3.331293  gm/ca**3 

RHO  =  3. 331293 
A.TJEKEV  =  AMU  ( EKEV ) 

MFP  *  AMUEKEV  ♦  RHO  *123. 

0 . DISTANCES  ARE  IN  METERS ,  M'J  IN  INVERSE  CENTIMETERS,  XJ  IS  IN  MFP3 

TOTKRN  =  3. 

HI  =3 

DO  u  3  5  J  a  1,13 
XJ  =  DMFP(J) 

TGL  *  TOLER  ( J ) 

C***»*  IF  R 1  IS  GREATER  THAN  XJ  DO  NOT  INTEGRATE  BACKWARDS 
XJMi  «  R1  *  MFP 
IF (XJ.LE.XJMi;  GOTO  S35 
N 1  =  N 1  *  1 
D2  =  XJ/MFP 

IF(Nl.EQ.l)  THEN 
XXJM1  *  XJMI 
ELSE 

XJMI  =  DMFP(J-l) 

END  IF 

Dl  =  XJMI /MFP 

CALL  SIMPS  (Dl ,D2,TOL, KERN, DISK) 

DS  *  XDOSE  *  KERN 
TOTKRN  «  TOTKRN  +  KERN 
90S  CONTINUE 
RETURN 
END 
C 

FUNCTION  GROUND  (R) 

C**»*#***  POINT  KERNAL  FOR  A  DISK  SOURCE 

C*««*t**t  This  function  calculates  the  kernal  for  a  cylindrical  wall 
C***t*#**  with  a  ground  roughness  shielding  factor 
REAL  HYP,  MUCONC ,  MFPGRN,  MUC.INDIR 
REAL  MFPWAL  ,  I  COS,  MU A IR, MFPT ,MFPAIR , MFPBU 
COMMON  /ENERGY/  EKEV , W,  PSI,  DOVERW 
COMMON  /SIZE/  RTHICK , WTHICK , BHGHT , BRAD 

C . STANDARD  REFERENCE  POINT  FOR  A  DOSE  IS  1  METER  FROM  GROUND 

PI  ■  3.14 1592654 
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HEIGHT  »  1. 

HYP  «  SORT  (HEIGHT**2. ♦  R*R ) 

I  COS  ■  HYP/R 

MUAIR  =  AMU  < ZKEV )  *  3.331293  *  133. 

MFPA I R  *  MUAIR  ♦  ICOS  #  (R-3RA0-.2) 

C - Assume  mall  is  .2  meters  thick 

RHODRT  «  1.3 

SPEC.  GRAY.  OF  DIRT  IS  APPROX  1.3  SOURCE!  MECH.  ENG.  HANDBOOK 
MUC  «  MUCONC ( EKEV ) 

MFPWAL  »  ICOS  *  MUC  #  MTHICK 
MFPBU  ■  MUAIR  #  SRAD 
MFPT  «  MFPAIR  +  MFPWAL  ♦  MFPBU 
IF (MFPT. GT. 12. )  THEN 
GROUND  =  0. 

RETURN 

END  IF 

C+++++  Find  BETA  from  R  and  Height 
IF (R. GT. 0. )  THEN 

SETA  *  ATAN (HE IGHT /R ) 

ELSE 

BETA  •  PI/2. 

END  I F 
C 

TAVE  «  0. 

A  a  W 

DIRECT  »  1. 

INDIA  «  3. 

IF (W. LE. 2. )  GOTO  20 

C  +  +  +  +  +  if  BETA  is  greater  than  Psi  ,  there  is  of  rough  ground  effect 
IF (BETA. GT. PSI )  GOTO  20 
C+++++  find  A  from  M,  PSI  and  BETA 

C+++++  A  is  that  part  of  the  trough  that  is  not  in  shaddow 
A  -  W  *  SIN(RETA)  *  COS  ( PS  I ) /SIN (BETA  +  PSI) 

C+++*»  Find  the  average  thickness  of  earth  inairect  part  passes  throu 
TAVE  ■  M/4.  *  SIN(PSI) / (COS  (PSI)*5IN(BETA) ) 

130  CONTINUE 

D  *  DOVERW  *  W 

DIRECT  ■  (A  +  D) / (W  +  D) 

INDIR  «  (W  -  A)/  (W  +  D) 

MFPGRN  ■  MUC  *  RHODRT  #  TAVE 
IF (MFPGRN. GE.  12. )  THEN 
TERAIN  =  0. 

GOTO  20 

END  I F 

C--“ - Find  Build-Up-Factor  for  fraction  that  encounters  the  ground 

BUF 1  «  BCONC(EKEV, MFPGRN) 

BUF2  ■  BUF ( EKEV , MFPAIR  +  MFPGRN) 

BUF3  «  BUF (EKEV f MFPGRN ) 

BUF4  *  BCONCtEKEV, MFPGRN  +  MFPA1 R+MFPWAL ) 

BUF5  ■  BCONC ( EKEV , MFPGRN  ♦  MFPAIR) 

BUF6  »  BUF(EK£V,MFP3RN+MFPT) 
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BUF7  =  BUF(EKEV,MFPGRN+MFPAIRrMFPWAL) 

3UT 1  «  SUFI  *  SUF2/BUF3  #  SUF4/8UF5  *  9UF6/3UF7 
TEPAIN  *  INDIP  *  EXP  i - ( MFP6RN+MFPT ) )  *  BUT  1 
23  CONTINUE 

C -  Find  Build-Up  Factor  Tor  Traction  that  missed  the  around 

BL'Fl  =  BUF  ( EKEV, MFPAIR) 

3  Ur  2  =  3CCNC(EKEV,MF?UAL+MFPAIR) 

DUF3  =  3C0NC  (EKEV , MFPA I R) 

3 U F 4  »  BUF (EKEV, MFPAIR+MFPUAL+MFPBU) 

BUF5  =  BUF ( EKEV, MFPA IR+MFPWAL) 

3 U T 2  =  BUF 1  *  BUF2/3UF3  *  BUF4/3UF5 
D I S K 1 =  DIRECT  *  BUT2  *  EXP(-MFPT) 

GROUND  *  .5  *  R/HYP**2.  *  (DISKI  +  TERAIN) 

i'ETURN 
END 
C 

FUNCTION  DISKRFIR) 

C»»******  POINT  KERNAL  FOR  A  DISK  SOURCE 
C********  This  Function  calculates  kernel  Tor  rooT 
REAL  HYP, ZSINf MUCCMC, MFP 
COMMON  /ENERGY/  EKEV,  W,  PSI ,  DOVERW 
COMMON  /SIZE/  RTHIC:<,NTHICK,BHGHT,3RAD 
IF (BH3HT.LT. 1 . )  THEN 

PAUSE  'PROGRAM  ABORT  --  BUILDING  TOO  SHORT’ 

END  I F 

A  =  BHGHT  -1.  +  .2 

C . STANDARD  REFERENCE  POINT  FOR  A  DOSE  IS  1  METER  FROM  GROUND 

HYP  =  SGRT ( A*A  +  R  *  R ) 

ISIN  =  HYP/A 

AIRMFP  a  AMU  (EKEV)  *  HYP  *  .301293  *  130. 

CONMFP  =  MUCONC  (EKEV)  *  RTHICK  *  ISIN 
MFP  =  AIRMFP  +  CONMFP 
IFIMFP.3T.33. )  THEN 
DISKRF  =  3. 

RETURN 

END  I F 

3U  »  BCONCfEKEV, CONMFP)* BUF(EKEV,MFP)/BUF(EKEV, CONMFP) 

DISKRF  »  .3  ♦  BU  *  R/HYP**2.  *  EXP(-MFP) 

RETURN 

END 

C- 

FUNCTION  DISKIN (R) 

C****«*»*  POINT  KERNAL  FOR  A  DISK  SOURCE 

C****«*»*  This  Function  calculates  kernal  Tor  in-finite  smooth  plane 
REAL  MFP,  HYP 

COMMON  /ENERGY/  EKEV,  W,  PSI,  DOVERW 

C . STANDARD  REFERENCE  POINT  FOR  A  DOSE  IS  I  METER  FROM  GROUND 

HYP  »  SORT ( 1 .  +  R*R) 

MFP  =  AMU (EKEV)  *  .001293  *  100.  *  HYP 
BU  »  BUF ( EKEV ,  MFP) 

SA  ■  1. 


144 

DISKIN  =  .5  *  BU*SA*£XP(-MFP)*R/HYP**2. 

RETURN 

END 

C- 

SU3RQUT I NE  SIMPS (LOWER, UPPER,  TOL,  SUM ,  F) 

C  NUMERICAL  INTEGRATION  BY  SIMPSON'S  RULE. 

C 

INTEGER  IN,  OUT,  PIECES,  I,  P2 
REAL  X,  DELTA,  LOWER,  UPPER,  SUM,  TQL 
REAL  ENDSUM,  ODDSUM,  SUM  1 ,  EVSUM 
COMMON  /  I NOUT /  IN,  OUT 
PIECES  *  2 

DELTA  =  (UPPER  -  LOWER)  /  PIECES 
ODDSUM  »  F  (LOWER  +  DELTA) 

EVSUM  =  3.3 

ENDSUM  »  F  (LOWER)  +  F(UPPER) 

SUM  =  (ENDSUM  4  *  ODDSUM)  *  DEwTA  /  3.3 
5  PIECES  =  PIECES  *  2 

P2  =  PIECES  /  2 
SUM  1  *  SUM 

DELTA  a  (UPPER  -  LOWER)  /  PIECES 
EVSUM  *  EVSUM  +  ODDSUM 
ODDSUM  a  0.3 

DO  10  I  »  i,  P2 

X  a  LOWER  +  DELTA  *  (2  *  I  -  1) 

ODDSUM  a  ODDSUM  +  F ( X ) 

10  CONTINUE 

SUM  »  (ENDSUM  ♦  4.0  #ODDSUM  ♦  2.0  *  EVSUM)  ♦  DELTA  /  3.0 
C**»*  WRITE  (♦,130)  LOWER, UPPER, SUM, PIECES 
103  FORMAT (IX,  *  LOWER ' , G 1 4 . 7 ,  ’  UPPER ' , G14. 7, *  SUM ',01 4. 7,  'PIECES' , 13) 

IF ( ABS  ( SUM  -  SUM  1 !  .ST.  ABS  (TOL  ♦  SUM))  GOTO  5 
RETURN 
END 
C 
C 

C- 

FUNCT I  ON  AMU(EKEV) 

C  MASS  ATTENUATION  IN  cm*2/gr 

REAL  £(26)  ,  MUQRH0<26) ,  LOGINT 

C . AIR  MASS  ATTENUATION  DATA  FROM  RADIOLOGICAL  HEALTH  HANDBOOK 

C . ,  JAN  1973,  PG139 

C 

DATA  £/  10. ,  15. ,20. ,30.  ,40.  ,30. ,60. ,80. ,100. ,  150. ,200. ,300.  , 

1  400. ,300. ,600. ,800. ,  1000. , 1500. ,2000.  ,3000.  ,4000.  ,5000. , 

2  6300. ,8000. ,10030. ,15000. / 

DA’ A  MUORHO/  4.99,1.55,0.752,0.349,0.2^3,0.209,0. 188,0. 167,0. 154, 

1  0.  1  36,0.123,0.  107,0.0954,0.0870,0.0805,0.0707,8.0636, 

2  0.3513,0.0445,0.0353,3.0308,0.0275,0.0252,0.0223, 

3  0.0204,0.0131/ 

C 
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C - USE  LOG-LOG  INTERPOLATION 

LOG  I  NT ( X , X 1 , X2 , YI , Y2)  =  EXP  ( ( ALDG ( Y2) - ALOS ( V I ) ) * 

1  (ALCG(X)-ALOGai)  ) / (ALOG (X2) -ALGG  (XI )  )  +  ALOS  ( Y1 )  * 
C 


C - use  lowest  two  data  points  for  s n e r g i e s  below  table 

IF t  E  K  E  V .  L  T .  £  ( 1 ) )  GO  TO  212 
203  IF (EKEV  -  E  !  J!  )  212,211,213 
213  J  =  J  +  1 

IFtJ.LT.26)  GOTO  233 

C - use  highest  two  data  points  for  energies  aoove  table 

GOTO  212 

211  A ,1U  =  M'JORHO  ( J ) 

RETURN 

212  CONTINUE 

All)  =  L03INT  (EKEV ,  E  ( J  —  i )  ,  E  !  J :■  ,  MUORHO  <  J- 1 )  .MUORHO (J)  ) 

RETURN 

END 

C- 

FUNCTION  R0ENT3 (EKEV ! 

r>  v  ■/  •/  •/  •/  v  v  v  »/  '/  4/  v  v  v  •/  v  •/  •/  •/  •/  ♦/  v  *'  •/  v  V  v  v  v  v  */  v  y  v  v  v  y  •/  V  y  v  v  V  »/  »/  v  V  v  •  • «/  f/vv  ■/  f  v  «  v  v  y  y  •/  v  •/ '/  y  V  V  V  •/ 

REAL  E  (26) ,  MUGRHO ( 26) ,  LCGINT 

C - Air  mass  absorption  data  from  NSRD5-N3S  29,  196?  pg  23  U  21 

DATA  E/  13.  ,  15.  ,20. ,33. ,43. ,53. ,60. ,83. , 133. , 153. ,203. ,303. , 

1  433. ,533. ,633. ,933. , 1333. , 1533. ,2333. ,3333, ,4333. ,5333. , 

2  6333. ,9333.  ,  13233. , 1  3333. / 

DATA  1UCRHG/  4. 61, 1.23, 3. 511, 3. 148, 3. 3669, .3436, .0335, .024 3, 

1  .3234,  .3253, .0263,. 32S8, .0295,  .0296,  .3295,  .3299,  .3273, 

2  .3254,.  3234,. 3235,.  01 36,. 31  74,.  3 164,.  3 1  52,.  3 1  45,.  3 1  32/ 


L - USE  LOG-LOG  INTERPOLATION 

L0GINT(X,X1,X2,Y1,Y2)  *  EXP  ( (ALOG (Y2 ) -ALOG t Y 1 )) ♦ 

1  ( ALOS  f  X  > -ALOS  (XI) ?/ (ALOG ( X2) -ALOG ( X 1 ! )  +  A L 0 G  ( Y 1 ) ) 

C 
C 

J  =  2 


0 - use  lowest  two  data  points  for  energies  below  table 

IF(EKEV.LT.Ed)  )  GO  TO  212 
233  IF (EKEV  -  E  ( J)  )  212,211,213 
213  J  *  J  +  1 

IFtJ.LT. 26)  GOTO  233 

C - use  highest  two  data  points  for  energies  above  table 

GOTO  212 

211  R0ENT3  =•  MUORHO  ( J !  *  EKEV  *  244.3365  / 1003. 

RETURN 

212  CONTINUE 

RQENTG  »  LOG  I  NT  (EKEV,E(J-1) ,E(J) ,MU0RHQ(J-1) , MUORHO  ( J ) )  * 
1  EKEV  #  244.3365  /  1 020 . 

RETURN 


END 
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C- 

FUNCT I  ON  BUF (EKEV.MUR) 

SEAL  £  ( 25 )  ,  A 1  ( 25 )  ,  A2C2),  ALPHA1  (25)  ,  A  L  P  H  A  2  (25)  ,  ALPHA3  ( 25 ) 
SEAL  MUR ,  LOG  I  NT,  SEMINT 
COMMON  /OLD J /  EOLD,  JOLD 

DATA  E/ 15. ,23. ,33. ,43.  ,53. ,63. ,33. ,103.  ,153.  ,203.  ,333.  , 

1  433. .533. ,633. ,833. , 1333. ,  1533. ,2333.  ,3303.  ,4333.  ,5333.  , 

2  6333. ,3333. ,13333. ,15333./ 

DATA  A1  /  1.253E1,4.?63,1.039E1, 1. 183E2,5.  136E2, 

1  1.641 £3, 1.477E3,1.533E3,1.242E3,1.236E3, 

2  1 . 2  5 1 E  3 , 1 . 182E3,1.232E3.4.316E3,1.102E3, 

3  1.128E3,2.941E2,4.159E2,1.162E2,1.923E1, 

4  1. 25 1 E1  ,  1.34  7E 1,1. 31  IE  1 ,3.339,6.661/ 

DATA  A2'-3.984E-l,-6.3?5E*l,-6.?24t-l,-1.462E2, -6. 199E2, 

1  -2.  712E3.-2. 543E3.-2. 736E3 , -2. 317E3 , -2 .  149E3, 

2  -1.75aE3,-2.3!?E3,"!.664E3,-4.533E3,-1.333£3, 

3  ~1. 17  4E3, -2. 1 79E2 , -2.63702,-6. 743E 1,-1.  41  7E1, 

4  -6.271 ,-3.? 19, -3. 346,-2.363,-1.496/ 

DATA  ALPHA  1 /  — 2.539E  — 2,-1. 35SE-T , -3. 1 74E-2 , -2 , 352E-2 , -4 . 23 1 E-2 , 

1  -4. 333E-2,-7.333£-2,-3. 193E-2, -0.  536E-2  , -7.  792E-2 , 

2  -5. 541E-2.-3. 95'3E-2,-2. 343E-2, -1.751  E-2, -1.  141E-2, 

3  -1 .33lE-2,-4.964E-2,-3.734E-2,-2.395E-2,-2.575£-2, 

4  -4.  129E-2 , -5. 1 76E-2 , -4. 734E-2, -4. 339E-2 , -5 . 832E-2/ 

DATA  ALPHA2/3. 955E-1 ,5. 535E-1 ,9. 333E-1 ,6. 613E-3, -1 . 352E-2, 

1  -2. 539E-2 , -4. 1 39E-2  , -4 . 753E-2 , -5. 325E-2 , -4 . 77  3E-2 , 

2  -3. 475E-2, -1 . 6730-2,-1. 351 E-2,-1. 426E-2, -2. 14 IE-3, 

3  -5. 425E-3,-3.933E-2  ,-2.  974E-2 , -7. 339E-3 , 2 .  3S4S-2. 

4  5.  333E-2, 7.  3943-2,3.  32 1 E-2, 9. i46E-2, 1.31 3E-1/ 

DATA  ALPHA3/  -2 . 659E-2 , 3 . 524E-2 , 6.  434E-2 , 6 .  739E-2 , 6 . 537E-2 , 

1  6.54oE-3, 1. 45PE-3,-6.323£-3,-9.997E-3,-9.  999F-3, 

2  1.52SE-2,1.292E-2,2.636E-2,3.632E-2,4.243E-2, 

3  7. 5l3E-2,-6.392E-2,-4. 662E-2 , -3.  262E-2 , -4 .  132E-2, 

4  -S. 231 E-2, -6. 235E-2, -5. 432E-2, -5.524E-2, -  5.  774E-2/ 

C*##*#  BUILD  UP  FACTOR  DATA  FROM  NUC  30!  U  END  73  P374  1931 

3F  !  A  ,  3 , 0  ,  D  ,  E ,  UR)  =  A»EXP(-C*UR)  •*-B*EXP(-D*UR)  +  <  1-A-B)  *EXP(-E*UR> 

C - USE  LC3-L03  INTERPOLATION 

L D 3 : T  (  X  ,  X 1 ,  X 2 ,  V 1 ,  Y 2 )  =  EXP  ( ( AL03  ( Y2) -ALOG  ( Y 1 )  )  e 
1  i  A  L  3  3  <  X  >  ►,  L  0  G  <  X 1 )  )  /  (AL0GIX2)  -  ALC3  ( X 1 ) )  ♦  AI.OG(Yl)) 

SEMINT  (X .  XI ,  X2,Y1  ,Y2)  =>  EXP  (  (ALOG  <  Y2> -ALOG  ( V 1 ))  * 

1  ( X  -  X 1 ) /  (X2-X1 )  ♦  ALC3 ( Y1 ) ) 

DO  NOT  LOOK  UP  CONSTANTS  IF  YOU  ARE  NOT  CHANGING  ENERGY 
IF  (EKE1/,  £3.  EOLD!  THEN 
J  =  JOLD 
GOTO  212 
ENDIF 

EOLD  *  EXEV 
C 

J  -  2 

C - use  lowest  two  data  points  for  energies  below  table 

IF(EKEV.LT.Ed)  )  GO  TO  212 
233  IF  <  E  K  £  V  -  E(J>)  212,211,213 
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210  J  =  J  +  1 

I F ( J . LT . 25 !  SOTO  203 

C - use  highest  two  data  points  for  energies  above  table 

GOTO  212 

211  IF (MUR. GT. 43. !  GOTO  213 

BUF  =  BF ( A 1 (J) , A 2 ( J ) , ALPHA  1 (J) , ALPHA2  ( J ’ ,ALPHA3(J)  , M U R ) 

C  WRITE (1,930) EKEV.SUF, J,A1 ( J ) , A  2 ( J ) , ALPHA  1 (J) , A  L  P  H  A  2  ( J )  ,  ALPHAS  ( J ) 
JOLD  =  J 
RETURN 

212  CONTINUE 

IF (MUR. GT. 43. )  GOTO  213 

X2  =  BF  ( A 1  CJ) ,A2(J) , ALPHA  1 ( J ) , ALPHA2 ( J ) , ALPHA3 ( J ) , MUR ) 

X 1 =BF  ( A 1  ( J - 1 ) , A 2  <  J - 1 ) , ALPHA  1 (J-l) ,ALPHA2(J-1! , ALPHA3 ( J- 1 ) , MUR ) 

BUF  =  LOG  li'JT  (EKEV ,  E  ( J- 1 )  ,£(J)  ,  X 1  ,X2) 

JOLD  *  J 
RETURN 

213  CONTINUE 

C . IF  MFPS  ARE  GT  43  INTERPRET  FROM  END  OF  RANGE 

X2=  BF  v  A 1 t J) ,A2(J) . ALPHA  1 (J) , ALPHA2 ( J ) , ALPHAS (J) ,40. ) 

X 1  =  3 F ( A 1 (J-l) , A 2 ( J - 1 >  , ALPHA  1 (J-l) ,ALPHA2(J-1) , ALPHAS  ( J- 1 )  ,  40 .  ) 
BUF2  =  LQGINT  <EKEV,E ( J-l ) ,E(J) ,X1,X2) 

X2=  BF ( A 1 (J) ,A2(J) , ALPHA  1 (J) , ALPHA2 ( J ) , ALPHA3 ( J ) , 35. ) 

X1=BF(A1 (J-l) ,A2(J-1) , ALPHA  1 (J-l) , AI.PHA2 ( J - 1 ) , ALPHAS (J-l) ,35. ) 

BUF 1  *  LOGIN T(EKEV,E(J-1) ,E(J) ,X1,X2) 

BUF  =  SEM I  NT (MUR , 35. ,40. , BUF 1 , BUF2) 

JOLD  =  J 
RETURN 
END 
C- 

FUNCTION  3CQNC (EKE^, RMFP ) 

C . THIS  FUNCTION  CALCULATES  THE  BUILD  UP  FACTOR  IN  CONCRETE  OR 

C . CONCRETE  LIKE  MATERIALS  FROM  THE  NSS  SOLUTION  USING  12  PARAMETERS 

C . BY  E I SENHAUER  AND  SIMMONS  REF  NUC  SCI  AND  ENG  56,263-270,1975 

C 

C . R  IS  IN  MFPS 

C 

REAL  E  (25)  ,  DO ( 25 ) , MUC ( 25) , LOG  I  NT , MUCONC 
REAL  A  8  (25) , A 1 (25) ,A2(25) , A  3 (25) , A  4 (25) 

REAL  B  0  (25)  ,B1 (25) ,B2(25) , B  3 (25) ,94 (25) 

R E A L * 8  BFC 

COMMON  /0LDJ2/  E 0 L D 2 ,  J0LD2 

DATA  E/  15. ,20. ,30. ,40. ,50. ,60. ,80. , 100. , 150. ,200. ,300.  , 

1  400. ,533. ,600. ,800.  , 1300. , 1500. ,2000. ,3000. ,4000. ,5000. , 

2  6000. ,8000. ,10003. , 15000. / 

DATA  MUC/8. 01 00, 3. 4450, 1. 1190,0.5580,0.3608,0.2734,0.2004,0. 1704, 

1  0.1399,0.  1250,0.1373,0.0958,0.0873,0.0807,0.0709, 

2  0.0637,0.0519,0.0448,0.0365,0.0319,0.0290,0.0270,0.0245, 

3  0.0231,0.0215/ 

DATA  DO/0.1589,0.1487,0. 1327,0.1200,0.1128, 

1  0.1  1  18,0.  1215,0.1376,0.1791,0.2149,0.2690, 

2  0.3088,0.3409,0.3674,0.4035,0.4404,0.4961, 
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3  3.5283,3.5679,0.5938,3.6138,0.6241,0.6425,3.6563,0.6676/ 
DATA  A3/1. 1 135E-2, 2.  4283E-2, 0.36645,0.  1  1734,3.  15913, 

1  0.13644,0.22372,0.24631,0.27760,3.29820,0.31837, 

2  3.  32929,0.33559,3.33965,0.34413,0.34678,0.  351  10, 

3  3.35445,0.35691,0.35425,0.34788,0.33878,3.31686, 

4  0.29369,0.24038/ 

DATA  A 1  /-0. 03 1  1 583,-0. 0024613,-0.0063622,-0. 0092903, 

1  -0. 01 5356, +0. 025044,-0.010144,0.040674,0.021783, 

2  0.029893,0.012225,0.0095543,0.0044976,0.0073127, 

3  0.0091422,0.0098983,0.010016,0.0099221 ,0.0078632, 

4  0.010371,0.040628,0.330809,-0.013397,-0.016610,-0.0029308/ 
DATA  A2/5.483E-04, 1 . 2409E-03, 3. 9089E-03 , 9 . 6343E-03 , 

1  1 . 5 1 69E-02 , 1  .  1833E-02 , 3. 7832E-02 , 3. 2359E-02 , 2 . 9878E-32 , 

2  4 . 8446E-02  , 

3  5. 5031E-02 , 5. 0474E-02, 2.2319E-02,2.  7012E-02, 2.  2857E-02 , 

4  2. 0876E-02, 1 . 9/S5E-02  , 1 . 9206E-02,2.  1801E-02 , 3. 9089E-02 , 

5  -9 . 3943E-03 , -8.  8533E-03 , 2. 8838E-02 , 2. 726 1 E-02 , - 1 .  9 1 0 1 E-02/ 
DATA  A3/2. 1 503E-04 , 5. 0984E-04 , 1 .  7224E-03 , 4 .  0328E-03 , 

1  7.  5754E-03, 1 . 3153E-03 , 3. 0337E-02, 2. 6207E-02, 3. 1915E-02, 

2  8.2357E-02, 1.0486E-01,  1.  1 029E-01 , 3. 4039E-02 , 3 .  1744E-02, 

3  4.4123E-02,5.0435E-02,5.2275E-02,4.9103E-02,3.9008E-02, 

4  2.9853E-02,2.888BE-02,2. 5395E-02 , 1 . 8330E-02 , 1 . 4079E-02 , 

5  2.2B17D  02/ 

DATA  A4  23828-05,1. 0829E-04, 3. 9905E-04,1.3308E-03, 

1  2.  2527.-03,0.00000000, 1.2457E-02, 3.  151 4E-06, 5. 8041 E-02, 

2  -  1.460BE-04, -8. 64 1  1E-04, -6. 3551 E-08, 1.0539E-01 ,9.43828-02, 

3  7. 3809E-02, 5. 9498E-02, 4. 0714E-02.3. 2172E-02 , 2. 2187E-02 , 

4  2. 0070E-03, 6. 7348E-05, 3. 5701E-03, 1 . 5406E-03, 3. 6397E-07 , 

5  9 . 8532E-03 / 

DATA  B0/0. 57471, 0.56384, 0.54853, 0.53850, 0.53787, 

1  3.54758,0.57172,0.60062,0.65913,0.69464,0.73934, 

2  3.  76455,0.78166,0.  79363,0.83851,3.81671,3.  82277, 

3  0.81861,0.8046  7,0. 79201 ,0.78176,0.77269,0.75988, 

4  0.75116,0.73974/ 

DATA  B2/0. 72706, 0.71518, 0.69814, 0.68724, 0.68395, 

1  0.90637,0.69136,0.86298,0.77287,0.84634,0.82898, 

2  0.  82998,0. 76467,0. 79645,0.83369,0. 80482,0. 79882, 

3  0.  79546,0.81708,0.85749,0. 72880,0.  63341,0.35214, 

4  0.86119,0.55119/ 

DATA  B 1/0. 42393, 0.41 368, 0.39364, 0.360! 0,0. 20648, 

1  0.70718,0.80330,0.69185,0.68532,0.71979,0.68807, 

2  0.65925,0.56762,0. 58849.0.57175,0.55269,0.52360, 

3  0.53858,0.46877,0.30734,0.83930,0.83961,0.57962, 

4  0.56105,0.85134/ 

DATA  83/0. 87317,0.86506,0.85649,0.85396,0.85703, 

1  1.00436,0.85532,0.98931 ,0.87818,1.01431,1.01389, 

2  1.30620,0. 87305,0.90888,0.91929,0.91359,0.91265, 

3  0.90861 ,0. 91102,0.96145,0.96388,0.94894,0.95976, 

4  0.97377,0.86292/ 

DA™  B4/0. 96440, 0.96 154, 0.95995, 0.96 1  1 0,0.  96562, 

1  1.00000,0.98034,1. 17560,1.00645,1.09682, 1.07209, 
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2  1.34173,1.33139,0.99331,3.99455,0.99165,0.93767, 

3  3.93542,3.93436,1.31553,1.10393,1.33529,1.01131, 

4  1.22932,3.97741/ 

BFC(R,U,A3,A1,A2,A3,A4,33,31,32,B3,84,D0)  =  1  +  ( 

1  A0*U*R*D£XP(-DBLE(U*R/80) ) +A1 * (U*R/ 8 1 ) **2»DEXP < -DBLE ( U*R/ B 1  > )  + 

2  A2  »  (U»R/32)**2  *  DEXP (-D3LE ( U* R / B 2 ) )  + 

3  A3  ♦  (U*R/33) **2  ♦  DEXP ( -DBLE (U*R/33 ) )  + 

4  A4*(U*R/84)#*2*DEXP(-DBLE(U*R/84)  )  )  / ( DO*DE XP ( -DBLE ( U*R )  )  ) 

C - USE  LOG-LOG  INTERPOLATION 

LOG  I  NT ( X , X 1 , X  2 , Y 1 , Y2 )  =  EXP  (  (ALOG ( Y2 ) -ALOG ( V 1 )) ♦ 

1  (ALOG(X)-ALOG(Xl) ) / (ALOG ( X2) -ALOG (XI ) )  +  ALOG ( V 1 ) ) 

R  =  RHFP  /  MUCONC  ( EKE  V ) 

C  WRITE!*, 1333)  'BCONC  2  ' , EKE V 

C******  DO  NOT  LOOK  UP  CONSTANTS  IF  YOU  ARE  NOT  CHANGING  ENERGY 
IF (EKEV. EQ. E0L02) THEN 
J  =  J0LD2 
GOTO  212 
ENDIF 

E0LD2  *  EKEV 
C 

J  =  2 

C - use  lowest  two  data  points  for  energies  below  table 

IFIEKEV.LT. Ell!)  GO  TO  212 
230  IF (EKEV  -  E  ( J ;  !  212,211,210 

210  J  =  J  +  1 

IF ( J. LT, 25)  GOTO  200 

C - use  highest  two  data  points  tor  energies  above  table 

GOTO  212 

211  BCONC»BFC(R,MUC(J) ,A0(J) , A 1  ( J ) , A2 ( J ) , A3 ( J ) , A4 ( J ) , B0  ( J ) , B 1  ( J )  , 

1  B2(J) , S3 ( J ) , 84  <  J ) , DO  ( J ) ) 

URCONC  =  R  *  MUC1J) 

JOLD2  =  J 
RETURN 

212  CONTINUE 

X2=BFC ( R , MUC ( J ) ,A0(J) ,A1 (J) ,A2(J) ,A3(J) ,A4(J) ,80(J> ,81 (J) ,82 (J) , 

1  B3 ( J ) , B4 ( J ) , DO  ( J ! ) 

K  =  J  -  1 

X 1  =  BFC ( R , MUC (K ) ,A0(K) , A 1 (K) ,A2(K) ,A3(K) ,A4(K) ,B0(K) ,81 (K) ,B2(K>  , 
1  B  3  <  K ) ,B4(K) , DO (K) ) 

BCONC  =  LOGINT (£KEV,E(J-1) ,E(J) ,X1 ,X2) 

URCONC  =  R  *  LOGINT(EKEV,E(J-l) ,E(J) ,MUC(J-1) , MUC  < J  > ) 

•J  0  L  D  2  *  J 
RETURN 
END 
C- 


REAL  FUNCTION  MUCONC ( EKEV ) 

C - this  fucntion  calculates  mu  over  rho  For  concrete  in  cm2/gr 

REAL  E ( 25 ) ,  MUORHO (25 )  ,  LOGINT 

C - mass  attenuation  data  from  Nuc  Sci  It  Eng  56  pr  267,  1975 

C 

DATA  E/  15. ,23. ,30. ,40. ,50.  ,60., 80. ,100. ,150. ,200. ,300.  , 
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1  400. ,500. ,600. ,300. , 1000.  ,1500. ,2000. ,3000. ,4000. ,5000. , 

2  6300. ,3030. ,1  0330. ,15330.  / 

DATA  MUORH0/3. 0103, 3. 4453, 1.1180, 3. 5588, 0.3609, 0.2734, 0,2304, 

1  0.1704,0. 1399,0.1250,0.1373,0.0958,0.0873,0.0307,0.0709, 

2  0. 0637,0.0519,3.0448,0.0365,3.0319,0.0290,0. 0270,0. 0245, 

3  0.0231,0.0215/ 

C 

C - USE  LOG-LOG  INTERPOLATION 

LOG  I  NT ( X , X 1 , X2 , Y 1 , Y2)  *  EXP ( ( ALOG ( Y2 ) -ALOG ( V 1 ) ) * 

1  (ALOG(X)-ALOG(Xl) ) /  (ALOG  (X2) -ALOG  <X1 ) )  ♦ALQG(Yl)) 

C 

J  =  2 

C - use  lowest  two  data  points  for  energies  below  table 

IF(EKEV.LT.Ed) )  GO  TO  212 
200  IF ( E K E V  -  E  ( J  )  »  212,211,213 

210  J  =  J  +  1 

IF  < J. LT. 26)  GOTO  200 

C - use  highest  two  data  points  for  energies  above  table 

GOTO  212 

211  MUCONC  =  MUQRHO ( J ) 

RETURN 

212  CONTINUE 

MUCONC  =  LOGINT(EKEV,E(J-l) ,E(J) , MUORHO ( J- 1 ) , MUORHO ( J ) ) 

RETURN 

END 
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APPENDIX  C 

CLOUD  PROGRAM  DESCRIPTION  AND  LISTING 
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The  Cloud  program  calculates  the  DRFs  for  a  cloud 
source  when  the  detector  is  centered  at  ground  level  inside 
a  structure.  The  program  modal  is  discussed  in  Section  2.3. 
Requirad  input  includes  tha  equivalent  radius  of  tha 
structure,  and  tha  mass  thicknaas  of  tha  walls. 

When  tha  Cloud  program  is  run,  it  will  ask  for  input 
data  as  shown  in  Figure  29.  The  first  input  required  is  a 
title  and  date  for  the  output  listing.  The  program  then 
requires  the  equivalent  radius  to  tha  inside  of  the  assumed 
herni-spherical  structure,  followed  by  the  mass  thickness  of 
tha  wall.  Tha  program  than  asks  the  "Do  you  wish  to  specify 
a  maximum  cloud  height?"  fallowed  by  "Do  you  wish  to  specify 
a  maximum  cloud  radius?".  Normally  the  answer  to  both 
questions  would  be  no  (n).  A  yes  (y)  answer  to  either 
question  would  allow  the  study  of  the  effects  of  limited 
cloud  size  on  the  exposures  or  DRFs  (see  Section  2.7).  The 
program  then  calculates  the  resultant  DRF'm  for  the  PWR-2 
(RSS)  cloud  spectra,  the  1-Hr  bomb  fallout  spectra,  and  the 
TMI 
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Bleloud 

I  INPUT  TITLE/DATE  of  output 

'  Sa« 1 1  Hood  Houtt.  16  S»pt  198* 

I  INPUT  RADIUS  TO  INSIDE  OF  hEHI-5PhERE  IN  HETERS  5. 

INPUT  NASS  THICHNESS  OF  NALL  IN  gr/c»»»2  1*. 

DO  YOU  NISH  TO  SPECIFY  A  HAXIHUH  CLOUD  HEISHTT,  Y  or  N 

Y 

Input  HAXIHUH  CLOUD  HEI0MT  in  MtKi  1808. 

DO  YOU  NISH  TO  SPECIFY  A  HAXIHUH  CLOUD  RADIUS’’,  Y  or  N 

Y 

I  Input  HAXIHUH  CLOUD  RADIUS  in  Ntttrt  1500. 

V _ 


« 


« 


FIGURE  2B  IyBiS*l_£222ytf r_£RT_gi3gl A^_During_Initi*t i on_of 
£ISy^_EC2flE*!Dx 
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i  cloud  spactra  followed  by  the  CRFs  -for  energies  -from  15  KeV 

A 

to  15  MeV.  A  typical  output  listing  in  shown  in  Table  15. 
Tha  program  calculates  tha  unprotected  axposura  in  R 
1  par  Curia  o f  airborne  radioactive  material  par  cubic  meter. 

The  program  then  -find*  the  protected  exposure  from  tha 
exterior  cloud  source,  followed  by  tha  exposure  from  an 

]  interior  cloud  containing  1  Curie  per  cubic  meter  of 

j 

radioactive  matarial.  The  DRFs  are  calculated  from  t’^e 
unprotected  and  protected  exposures.  Tha  protected 
exposures  are  then  compared  to  the  protected  exposure  from  a 

a 

PWR-2  (RSS)  cloud  spactra.  Tha  last  column  of  data  is  tha 

exclusion  factor  discussed  in  Section  2.9.2.  The  program 

i  assumes  one  gamma  par  disintegration  for  tha  1  Hr.  bomb 

J 

fallout,  the  PWR-2  cloud  (RSS-Cloud)  and  the  individual 
gamma  energy  calculations,. 

j  Tha  program  is  written  in  Micrgsoft  FORTRAN  3.13.  This 

is  a  version  of  FORTRAN  77.  There  are  no  non-standard 
FORTRAN  commands  in  the  program.  The  program  was  run  on  a 
^  !22<5!f_l05iCy03Di5_PCof255iOD3l_£92!2yTilC  (DC)  which  uses  an 

3303  processor  and  an  3087  co-processor .  Tha  program  was 
run  under  the  HS-DC3  2.11  operating  system.  However,  it 
I  should  run  on  any  computer  capable  of  using  FORTRAN  77.  The 

program  takes  approximately  25  minutes  running  time. 

Caution,  the  program  is  written  is  FORTRAN,  the  input 
§  C29yiC22  decimal  points.  The  program  may  not  converge  to  an 
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answer  for  low  energies  (15  and  20  KaV)  for  wall  thickness 
greater  than  about  23  gr/cm3.  This  '.s  because  the  exposures 
ara  lass  than  S.43E-37,  which  is  tha  limit  of  the  single 
precision  data  typa  used  for  real  variables. 

The  MAIN  program  diracts  the  input  and  output  of  the 
program  data.  Subroutine  SIMPS  carries  out  the  integrations 
by  using  Simpson's  rule.  Tha  number  of  increments  i3 
automatically  increased  until  tha  desired  accuracy  is 
obtained.  Subroutine  PROTECT  speeds  the  integration  of  the 
hami -sphar i cal  source  by  dividing  the  source  into  shells 
which  converge  to  integrals  faster.  Function  DHEMI  finds 
the  kernel  for  the  unprotected  cloud  source.  Function  DSPHER 
calculates  the  kernel  for  the  protected  cloud  source. 

Function  AMU  finds  the  mass  attenuation  coefficient  for 
air.  Function  BUF  finds  the  buildup  factor  for  air. 

Function  MUCONC  finds  the  mass  attenuation  coefficient  for 
the  walls.  Function  3C0NC  finds  the  buildup  facto*-  for  the 
walls.  Function  ROENTG  finds  the  factor  for  converting  flux 
to  exposure. 

The  listing  of  the  Cloud  program  follows. 


3 


CHARACTER  CTITLE*33,Y0RN*1 
I N  T  Z  5  E R  IN,  OUT 
REAL  E (23) 

REAL  HMAX,RMAX,E3R0U?<:3,4)  , FRACT  (22,4) 

EXTERNAL  F,  D H E H I , DSPHER 
COMMON  /ENER3Y/EKEV 

COMMON  /SIZE/  WALLDS,W THICK, RMAX,HMAX,R1 
COMMON  /EXPOS/  TOTEX? ,  TOTUNP  ,  EXPOSE 
COMMON  /INCUT/  IN,  OUT 
COMMON  / 0 L D J /  EOLD,  JOLD 
COflMCN  /OLD J 2 /  E0LD2 ,  JOLD2 

DATA  E/  15. ,  23. ,33. ,43. ,53.  ,60.  ,93.  ,  133.  ,  153.  ,233.  ,333.  , 

1  4  23. ,533.  ,632. ,662. ,333. ,  1223. ,  1  173.  ,  1332.  , 1323.  , 

2  2323. ,3323. ,4333. ,3233. ,6333.  ,  3333.  ,1  3333.  ,  1  5333.  / 

DATA  ESROUP  / 53 ., 333 .,  753. , 1 533 .  ,  2533 .,  1 5*3 .  , 

1  25. ,  75. ,153. ,253. ,353. ,533. ,733.  ,933.  ,1  1  65.  , 

1  52  3. , 1333. ,2253. ,  2753. ,3533. ,4533.  ,5*0.  , 

26  33. ,2433. ,2193. ,  1863. ,  1553. ,899.  ,953.  , 

6  13. ,527. ,363. ,253. ,  233. , 191  .  ,  166.  ,81 .  ,23.  ,4*3.  , 

31 ., 164. ,234. ,364. ,637. ,662. ,723.  , 13*3.  / 

DATA  FRACT  /. 3693 ,. 3359 ,.  1  476 ,. 1 34 1 ,.  0732  ,  1 5*3 .  , 

1  .3271  ,.  31 37,.  3737,. 3476,.  0929,.  1373,.  1717,.  1627, 

2  .3889,  .0957, .3299, .  0397, .01 43,.  0042,  .0331 ,5*3.  , 

3  . 333239, .  0331  9, , 32164, . 33191  ,.  03127, 

4  .33113,. 33239, .00732,. 01 63, .233454, 

5  .  2223,. 33252,. 033 1  9,. 33637,.  2764,.  333637, 4*3.  , 

6  .  3696,. 33332,. 030254,. 03331 9,.  003363,. 03334  2,.  3333 16, 13*0.  / 


C»**************INPUT  TITLE  FOR  EACH  LISTING*********************** 

337  F0RMAT(1H3,25X, 'CLOUD  SOURCE  WITH  WALL  WRITTEN  22  JUNE  1934') 

C  +  +  +  +  +  CHAN3E  EPSON  F X -83  PRINTER  TO  COMPRESS  MODE 

OPEN  (1 , F I LE=  '  LPT  1 ' , STATUS* ' NEW ' , ACCESS* ’ DIRECT ' 

1  .FORM* 'UNFORMATTED' ,RECL»l> 

WRITE  (1)27,15 
CLOSE  (1) 

OPEN ( 1 , FILE3 ' LPT  1 ' ) 

C 

OUT  =  1 
EOLD  =  3.3 
EOLD2  »  3.3 
IN  «  0 
WALLDS  =  .2 

C - Assume  walls  and  roof  are  .2  meters  thick. 

C  7.  X  X  7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. X 1.  X  7.  X  X  7. 7. 7. 7. 7. 7. 7. 7. 7. 7.  X  X  7.  X  X  X  X  7.  X  X  7.  X  7.  X  7.  X  X  X  7.  X  X  X  X  X  7.  X 
WRITE  (*,  293) 

293  FQRMATUX,  'INPUT  TJTLE/DATE  OF  OUTPUT') 

READ(*,291)CTI TLE 

291  FORMAT  (A83) 

292  FORMAT ( 1 H0 , 25X , A33) 
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1 


n 


*4 


1 


3 


1 


1 


I 


C  X  7. 7.  %  7. X 7. 7. 7. 7. 7. 7. 7. 7.  X  X  7.  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X 


303  WR I TE ( * , 30 1 ) 

331  FORMAT! IX, ' INPUT  RADIUS  TO  INSIDE  OF  HEMI-SPHERE  IN  METERS  '  ,  \ ) 

332  READ  (*,303)  R1 
303  FORMAT (F13. 6) 


WRITE  (*,304) 

304  FORMAT  < 1 X ,  'INPUT  MASS  THICHNESS  OF  WALL  IN  gr/cm*+2  ',\> 

READ  (*,303)  WTHICK 

293  WR I TE  ( * ,  294 ) 

294  FORMAT  (IX, 

1  'DO  YOU  WISH  TO  SPECIFY  A  MAXIMUM  CLOUD  HEIGHT?,  Y  or  N  * ) 

RE  AD ( * , 295 ) YCRN 

295  FORMAT  ( A 1 ) 

I F ( YORN . EQ . 'Y'.OR. YORN.EQ. 'y'>  THEN 
WRITE!*, 296) 

296  FQRMATUX,  ’Input  MAXIMUM  CLOUD  HEIGHT  in  meters  '  ,  \ ) 

READ (* , 333) HMAX 

ELSEIF (YORN. EQ.  'N  ' . OR. YORN, EQ. 'n ’ )  THEN 


297 

298 


299 


HMAX  =  1 E+33 
ELSE 

GOTO  293 
END  I F 
WRITE!*, 298) 

FORMAT (IX, 

1  'DO  YOU  WISH  TO  SPECIFY  A  MAXIMUM  CLOUD  RADIUS?,  Y  or  N ' ) 

RE  AD  <  * , 295) YORN 

IF (YORN. EQ.  '  Y  ' . OR. YORN, EQ. ' y ' )  THEN 
WRITE!*, 299) 

FQRMATUX,  'Input  MAXIMUM  CLOUD  RADIUS  in  meters  '  ,  \  > 

READ  ( * , 303 ) RMAX 

ELSEIF (YORN. EQ.  ' N  ' . OR. YORN. EQ. *n ' >  THEN 
RMAX  =  1E+30 
ELSE 

GOTO  297 
END  I F 

WR  I  TE ( OUT , 807) 

WRI TE (OUT , 292) C  TITLE 
WR I TE ( OUT  ,  808 ) 

WRITE (OUT , 809)  R1 
WRITE  (OUT, 813)  WTHICK 
WRITE (OUT, 81 1 )  HMAX 
WRITE (OUT, 812)  RMAX 
WRI TE  (OUT , 303 ) 

WRITE (OUT , 804 ) 

C8808*83988008§88809088000833908888a38998339838888808a880808@3000@@99@0 
C88983  ROUTINE  TO  FIND  DRFS  FOR  SPECTRA 
DO  503  NSPEC  *  1 , 3 , + 1 
DOSUNP  *  0. 

DOSCLD  *  3. 

DOSINT  ■  0. 

DO  501  I  EKE V  *  1,20 
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EKE V  *  E3RCUP  ( IEi'EV,NSPEC) 
r  ROT  *  FRACT  ( IEKEVt.M5PEC) 

IF (FRCT. EQ. 2)  GOTO  531 
WRITE (*,336)  EKEV 

C - LET  YOU  KNOW  PROGRAM  IS  RUNNING 

CALL  PROTECT (EKEV) 

CALL  GIMPS (3. , R 1 , 1 . E-5 , DINT , DHEMI ) 

DOSUNP  =  DOSUNP  +  FRCT  *  TOTUNP 
D03CLD  =  DOSCLD  +  FRCT  *  TCTEXP 
DOS  I N T  *  DCS INT  +  DINT  *  EXPOSE  *  FRCT 
531  CONTINUE 

DRF  =  S0SCLD/D3SUNP 
IF(NSPEC.EQ. 1)  DRF5M3  =  DRF 
RATIO  =  DRF  /  D  Ft  F  3  M  3 
EXCLUDE  =  DOSCLD/ (DOSCLD  +  DCSINT) 

IF (NS? EC. EQ. 2) WRITE (OUT ,535) DOSUNP,  DOSCLD,  DC  SI  NT, 

1  DRF , RATIO, EXCLUDE 

IF (NSPEC.EQ. 1) WRITE (OUT, 5131 DOSUNP, DOSCLD, DOS  I  NT, 

1  DRF, RATIO, EXCLUDE 

IF ( NSPEC. Efl. 3) WRITS (OUT, 536) DOSUNP,  DOSCLD, DOSINT, 

1  DRF, RATIO, EXCLUDE 

IF (NSPEC.EQ. 4) WRITE (OUT, 537) DOSUNP,  DOSCLD,  DCSINT, 

1  DRF, RATIO, EXCLUDE 
533  CONTINUE 

C33339933899939§§99399393333333,3993393333339939399393§83§3399330933330 
C 3 3 3 3 9  routing  For  ggngrtating  energy  table  93333993333399899335389333 
DO  331  I EKEV  »  23, 1,-1 
EKEV  =  EilEKEV) 

WRITE  (*,936)  EKEV 

336  FORMAT  (IX, 'E(KsV)  =  '  ,  G 1 3 . 5 , \ ) 

C . LET  YOU  KNOW  PROGRAM  IS  RUNNING 

DOSUNP  =  3. 

n  n  -  r  i  -  -» 

w  U  j  L  b  J  u  » 

DOSINT  =  3. 

CALL  PROTECT (EKEV) 


CALL  SIMPS (2. ,R1 , 1. E-5, DINT, DHEMI) 
DOSUNP  =  DOSUNP  +  TOTUNP 
DOSCLD  =  DOSCLD  +  TOTEXP 
DOSINT  =  DOSINT  +  DINT  *  EXPOSE 
DRF  =  DCSCLD/CCSUNP 
RATIO  *  DRF/DRF3MS 


EXCLUDE  =  DOSCLD/ (DOSCLD  ♦  DOSINT) 

ERRLIM  =  1 . £-33 
IF (DOSUNP. GT. ERRLIM)  THEN 
DRF  =  DOSCLD/DOSUNP 
ELSE 


DRF  *  3. 

END  IF 

QNEMEX  »  1.  -  EXCLUD 

WRITE (OUT, 832)  EKEV , DOSUNP , DOSCLD , DOS  I  NT 
1  ,  DRF  ,RA  ."10,  EXCLUDE 


1  SB 


801  CONTINUE 

333  CONTINUE 

332  FORMAT  ( 1 X , 25X , F7. 1  ,4X,6(GrJ.4,4X)  ) 

333  FORMAT ( 1 H0 , 25X , 'E(KeV) ' ,5X, 

1  'R/nr/Ci*m3',4X,'R/hr/Ci*m3',4X,'R/hr/Ci*m3’,4x, 

2  '  DRF  ' , 9  X , ' R  /  R  !  R  S  S !  '  ,6X, 

3  'Exclusion') 

334  FORMAT ( IX , 25 X , 12X, 'Unprotected  '  ,3X, 'Protected ' ,5x , 

1  'from  i  nt 5X  ,  1 4X  ,  1 4X Factor  ' ) 

928  FORMAT  * 1H0,25X Integration  to  12  MFPs  in  air.') 

809  FORMAT  <1X,25X,  'Equivalent  radius  of  structure  *  '  ,  F  9 . 3 ,  '  meters') 
813  FORMAT ( 1 X , 25X ,' Mass  thickness  ot  walls  =  '  , F  9 . 3 ,  '  gr/cm**2 ' ) 

811  FORMAT ( 1 X , 25X ,' Max i mum  cloud  height  *  '  ,  F 1 2 . 3  ,  '  meters') 

312  FORMAT  ( 1 X , 25X ,' Max i mum  cloud  radius  =  '  , F 1 2 . 3 ,  '  meters') 

535  FCRMAT(1X,25X,'1.12  HR.  '  ,  2X ,  6  ( 8 1 3 . 4 , 4X ) ) 

513  FORMAT ( 1 X , 2 5 X , 'R5S-CL0UD'  ,2X, 6(313. 4,4X1 ) 

536  FORMAT!  1 X , 25X ,  'TMi-CLOUi)'  ,2X, 6(313. 4, 4X1  ) 

537  FORMAT ( IX, 25 X, 'SL1 -CLOUD'  ,2X,6(313.4,4X)) 

STOP 

END 

C 

SUBROUTINE  PROTECT (EKSV) 

REAL  TCIER(13) ,DMFP(10! ,KERN,MFP 
EXTERN  1L  DHEM I , D3PHER 
COMMON  /SIZE/  WALLD3,WTHICK,RMAX,HMAX,Rl 
COMMON  /EXPOS/  TQTEXF  ,  TOTUNP , EXPOSE 
DATA  T0LER/1.E-4, l.E-4,  l.E-4, l.E-4,l.E-4, 1 .E-4, 

1  l.E-4,l.E-4,l.E-3,l.E-2/ 

DATA  DMFP/,003,.01  ,,05,. 1  , . 2, . 5, 1. ,3. ,5. ,12./ 

C . RHO  3  DENSITY  OF  AIR  AT  3TP  0.001293  gm/cm**3 

RHO  -  0.301293 
AMUtKEV  *  AMU(EKEV) 

MFP  =  AMUEKEV  *  RHO  *100. 

EXPOSE  =  RCENTG ( EKEV ) 

C . DISTANCES  ARE  IN  METERS,  MU  IN  INVERSE  CENTIMETERS,  XJ  IS  IN  MFPS 

TOTEXP  *  3. 

TOTUNP  *  3. 

N 1  =0 
N  2  3  3 

DO  305  J  =  1,13 
XJ  *  DMFP(J) 

TOL  *  TOLER  (  J  ) 

C***a*  IF  R 1  IS  GREATER  THAN  XJ  DO  NOT  INTEGRATE  BACKWARDS 
X J M 1  «  ( R 1  ♦  WALLDS!  #  MFP 
IF(XJ.LE.XJMl)  GOTO  805 
N 1  ■  N 1  ♦  1 
D2  «  X J /MFP 

I F ( D2. GT. RMAX )  THEN 
02  »  RMAX 
N2  «  N2  ♦  1 


END  IF 


IF(Nl.cQ.l)  THEN 
X  X  J  M  1  =  X  J  M  i 
X2JM1  =  a. 


X  J  M 1  =  DMFP(J-l) 

X 2 J M 1  *  DMFP(J-l) 

END  IF 

D 1  =  XJ.il/MFP 
D21  =  X2JM1/MFP 
D22  =  D2 

CALL  SIMPS ( DI ,D2,T0L, KERN, DSPHER) 
CALL  SI1PS(D21 ,D22)T0LIUN£HLDtDHEMI) 
TCTEXP  =  TOTE XP  +  EXPOSE  ♦  KERN 
TQTUNP  =  TOTUNP  +  EXPCSE  ♦  UNSHLC 
335  CONTINUE 
RETURN 
END 


SUBROUTINE  SIMPS (LONER,  UPPER,  TOL,  SUM,  F! 


C  NUMERICAL  INTEGRATION  BY  SIMPSON'S  RULE. 
INTEGER  IN,  OUT,  PIECES,  I,  P2 
REAL  X,  DELTA,  LOWER,  UPPER,  SUM,  TOL 
REAL  ENDSUM,  ODDSUM,  SUMI,  EVSUM 
COMMON  /INCUT/  IN,  OUT 
C 


PIECES  «  2 

DELTA  *  (UPPER  -  LOWER)  /  PIECES 
ODDSUM  »  F  (LOWER  +  DELTA) 

EVSUM  *  3.0 

ENDSUM  »  F (LOWER)  +  F  (UPPER) 

CUM  =  (ENDSUM  ♦  4  *  ODDSUM)  *  DELTA  /  3.0 
I  PIECES  =  PIECES  ♦  2 

P2  *  PIECES  /  2 
SUMI  =  SUM 

DELTA  =  (UPPER  -  LOWER)  /  PIECES 
EVSUM  »  EVSUM  +  ODDSUM 
ODDSUM  =  0.0 
DO  10  I  *  1,  P2 

X  »  LOWER  +  DELTA  *(2*1-1) 

ODDSUM  *  uDDSUM  +  F(X) 

10  continue 

SUM  «  (ENDSUM  +  4.0  *ODDSUM  +  2.0  *  EVSUM) 

*  »  DELTA  /  3.0 

IF ( AB3 (SUM  -  SUMI)  .GT.  ABS (TOL  *  SUM))  GOTO  5 

RETURN 

END 


C 

FUNCTION  DHEMI  (R) 

POINT  KERNAL  FOR  A  DISK  SOURCE 
REAL  MFP 
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COMMON  /ENERGY/  EKEV 
COMMON  /SIZE/'  HALLOS, T,RMAX,HMAX,R1 
MFP  =  AMU(EKEV)  *  3.331293  *  130.  ♦  R 
3U  =  9L'F  ( EK'EV ,  MFP) 

SA  =  1 . 

D  H  EMI  =  . 5*3U*3A*EXF (-MF?) 

IF(S.ST.HMAX)  DHEMI  -  OHEMI  *  HMAX/R 
RETURN 
END 
C 

FUNCTION  AMU(EKEV) 

C  MASS  ATTENUATION  IN  c.**2/gr 

REAL  E  (26)  ,  MUORHO (26) ,  LOGIN'’ 

C . a  I R  MASS  ATTENUATION  DATA  FROM  RADIOLOGICAL  HEALTH  HANDBOOK 

C . ,  JAN  1973,  P3139 

C 

DATA  E/12. ,15. ,  ?  ? . ,33. ,43.  ,13. ,63., 33.  ,133.  ,153. ,203. ,393. , 

1  43?  ,  533. ,  600. ,Sb3. ,  1333,  ,  1533. ,2330.  ,  3330. ,  4330. ,5303. , 

2  6333. ,3333. ,13233. , 15223. / 

DATA  MUORHO/  4.99  1. 55, 3. 732, 3. 349, 3. 249, 3. 233, 3. 1SB, 3.  167,3.  154, 

1  3.1  36,3.  123,3. 1  0 ■ ,  3 . 3954 , 3 . 0870 , 3.  3333 , 0 .  3737 , 0 . 0636 , 

2  3.351  3,3.  3445,3. 3353,3.  3333,0.3275,3.3252,3.  3223, 

3  3.3234,3.3131/ 

C 

C - USE  LOG-LOG  INTERPOLATION 

C**t**LININT(X,Xl,X2,Yl,Y2)  =  ( Y2-Y1 ) * ( X-X 1 )  /  (  X2- X 1 )  +Y1 
LQ3INT ( X , X 1 , X 2 , Y 1 , Y 2 )  =  EXP  (  ( AL03 ! Y2) -ALOG ( Y 1 ) ! * 

1  'ALOG(X) -ALCG(Xl) ) / (AL0G(X2) -ALOG (XI ) )  ♦  ALOG ( Y 1 ) ) 

C 

C 

J  *2 

C - use  lowest  two  data  points  /or  energies  below  table 

IF(EKEV.LT.Ed)  )  GO  TO  212 
230  I F ( £ < E V  -  Z  ( J )  !  212,211,213 
213  J  *  J  +  1 

I F ( J . L T , 2 6 )  GOTO  233 

Z - use  highest  two  data  points  /or  energies  above  table 

GOTO  212 

211  AMU  »  MUORHO ( J ) 

RETURN 

212  CONTINUE 

AMU  «  LOG  INT  ( E^EV ,  E  ( J  - 1 )  ,E(J)  ,  KUO:>HO  ( J-l )  ,  MUORHO  (J)  ) 

RETURN 

END 

C 

FUNCTION  DSPHER(R) 

C*»««»M#  POINT  kernal  for  a  sherical  source 
real  MU a  I R, MFP 1 ,MFP2,HFP3,MUCDNC,MFP 
INTEGER  IN,  OUT 
COMMON  /ENERGY/  EKEV 
COMMON  /SIZE/  HhLLDS,T,RMAX,HMAX,R1 


!  *.  ■ 


COMMON  /INCUT/  IN,  CUT 

C . ASSUME  RCOF  IS  .2  METERS  THICK 

MUAIR  =  AM'JiEKEV)  ♦  Z .  33 1 293 

MFP1  =  (R  -R I  -  xIAULDS!  *  MU A I R  »  133. 

MFP2  =  T  *  MUCCNC(EKEV) 

MFP3  *  R 1  ♦  MUAIR  ♦  133. 

MFP  =  f-FP  1  +  MFP2  ♦  MFP3 
B 1  *  8UF ( EKEV , MFP 1 ) 

B2  «  BCONC(EKcV,MFPl*MFP2)/3CONC(EKcV,MFPi) 

S3  =  BUF(EK£V,MFP1+MFP2*MFP3) /BUF<£KEV,MFP1tMFP2> 
BU  =  81*32*83 
SA  =  1. 

IF( (MFP2.5T. 43. ) .OR. IMF?1*MF?3. 3T. 43. ) )  THEN 
DSPHER  =  3. 

RETURN 

ELSE 

DSPHER  =  .5*BL'*SA«£TP(-MFP) 

END  IF 

IF (R. GT. UMAX) DSPHER  *  DSPHER  *  HMAX/R 

RETURN 

END 


cmx 


FUNCTION  R0T.MT6  (EKEV) 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


:  x  •/.  %  x  x  %  i  x  •/.  xxxx  i  %  xxxx  x  %  i  •/.  x  •/.  x  x 


C  This  function  converts  BUF  adjusted  flux  to  R*m**2/hr/Ci 

zxxxxxxxxx  %  x  r.  x  x  x  x  %  x  x  xxxxxxxxxxxxxxxxxx  x  x  x  x  x  x  x  x  xxxxxxxxxxxxxxxxxx  x  x  x  %  x  x  •/. 

REAL  E  ( 26 ) ,  MU0RH0'26) ,  LC3INT 

C - Air  mass  absorption  data  from  N5RD5-NBS  29,  1969  pg  20  &  21 

DATA  E/  10. , 15. ,23. ,30. ,40. ,33. ,60. ,80. , 100. , 150. ,200.  ,330.  , 

1  430.  ,330.  ,633.  ,333.  ,  1033.  ,  1533.  ,2303.  .,3330.  ,  4030.  ,5000.  , 

2  6333. ,8330. ,13033. ,15203./ 

DATA  MUORHO/  4 ,  6 1  ,  1 .  23 , 3 .  5 1  1 ,  0.  1  48 , 0. 0669 , . 0406 , . 0305 , . 0243  , 

1  .0234,. 0250,. 0263,. 3238,. 3295,. 0296 ,.0295,. 3239, .0278, 

2  .0254,. 3234,. 0235,. 01 86,. 0174,. 01 64,. 0152, .01 45,. 0132/ 


C - USE  LOG-LOG  INTERPOLATION 

LOG  I  NT  ( X , X 1  ,X2,Y1,Y2)  -  EXP i (ALOG (Y2) -ALOG  ( Y1 )  )  * 

1  (ALOG(X)-ALOG(Xl) )/ (ALOG (X2> -ALOG (XI ) )  +  ALOG  ( Y 1 ) ) 
C 


J=2 

C - use  lowest  two  data  points  for  energies  below  table 

IF (EKEV. LT. E ( 1 ) )  GO  TO  212 
230  IF (EKEV  -  E ( J) )  212,211,210 
210  J  =  J  +  1 

IF ( J. L7, 26 )  GOTO  200 


C - u 5 e  highest  two  data  points  for  energies  above  table 

GOTO  212 

211  R0F.NTG  =  MUORHO  (J !  *  EKEV  *  244.3365  / 1000. 

RETURN 

212  CONTINUE 

ROENTG  =  L0GINT(EKEV,E(J-1) , E ( J > , MUORHO ( J- 1 J , MUORHO ( J ) )  * 


l  EKEV  *  244.3365  / 1000. 
RETURN 
END 
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FUNCTION  BUF  (EKEV, MUR) 

REAL  E  (25)  ,  AK25),  A2  ( 25  >  ,  ALRHA1  (25)  ,  ALPHA2  <  25 )  ,  ALPHA3  ( 25 ) 
REAL  MUR ,  LOG  I  NT ,  SEMINT 
COMMON  / 0 L D J /  EOLD,  JOLD 

DATA  E/ 15. ,20. ,30. ,40. ,50. ,63. ,30. ,130.  ,  150.  ,200.  ,300.  , 

1  400. ,500. ,600. ,800.  ,  1000. , 1500. ,2000. ,3000. ,4000. ,5030. , 

2  6300. ,8033. ,13333. , 15303. / 

DATA  A1  /  1  .  258E1  ,4.960, 1.039E1, 1. 183E2.5.  106E2, 

1  1.641E3,1.477E3,1.503E3,1.242E3,1.206E3, 

2  1.251E3, 1. 1S2E3, 1.232E3,4.316E3,1.  102E3, 

3  1. 1 28E3 , 2 . 94 1 £2 , 4 . 1 59E2 , 1 . 1 62E2 , 1 . 928E 1 , 

4  1. 251E1 , 1.047E1 , 1.31 1E1 ,8.889,6. 661/ 

DATA  A2/-3. 994E-1  , -6 . 395E- 1 , -6. 924E- 1 , - 1 . 462E2 , -6 . 189E2, 

1  -2.7l2E3,-2.543E3,-2.736E3,-2.317E3,-2. 149E3, 

2  -1 . 756E3, -2. 319E3 ,  —  1  .  664E3, -4. 588E3, -1 . 308E3 , 

3  -1. 174E3 , -2. 179E2,-2.687E2,-6.740E1 .-1.417E1 , 

4  -6.371,-3.919,-3.346,-2.360,-1.496/ 

DATA  ALPHAl/-2.539E-2,-1.058E-3,-3.  1 74E-2 , -2 . 852E-2 , -4. 23 1 E-2 , 

1  -4.388E-2,-7.333E-2,-8.  1 90E-2 , -8.  536E-2 , -7 . 78PE-2 , 

2  -5. 5ClE-2,-3. 850E-2, -2. 843E-2,-l . 751E-2, -1 ■ 1 41E-2 , 

3  -1.331E-2,-4.964E-2,-3.784£-2,-2.395E-2,-2.575E-2, 

4  -4. 129E-2, -5. 1 76E-2 , -4 . 734E-2 , -4 . B39E-2 , -5 . 032E-2/ 

DATA  ALPHA2/3.955E-1 , 5 . 505E- 1 , 9 . 030E- 1 , 6 . 6 1 3E-3 , - 1 . 952E-2 , 

1  -2. 589E-2 , -4. 1 09E-2  , -4 . 753E-2 , -5 . 025E-2 , -4 . 773E-2 , 

2  -3. 475E-2 , - 1 .6788-2,-1.35 1 E-2, -1.406E-2, -2. 141E-3, 

3  -5.925E-3,-3.988E-2,-2.974E-2,-7.339E-3,2.054E-2, 

4  5.030E-2,7.3?4E-2,3.32!E-2,9. 146E-2, 1.018E-1/ 

DATA  ALPHA3/-2.659E-2,3.524E-2,6.484E-2,6. 739E-2 , 6 . 537E-2 , 

1  6. 546E-3, 1.458E-3, -6. 3238-3,-9. 9978-3,-9. 998E-3, 

2  1 . 528E-2, 1 . 292E-2,2. 686E-2.3. 682E-2,4. 240E-2, 

3  7 . 813E-2, -6. 392E-2,  -  4. 662E-2,-3. 262E-2, -4. 182E-2, 

4  -5. 231E-2  ,-6. 205E-2,  -5. 452E-2,-5. 524E-2 ,-5. 774E-2/ 

C*#***  BUILD  UP  FACTOR  DATA  FROM  NUC  SCI  &  END  78  PG74  1981 

BF ( A , B , C , D , E , UR) =  A*EXP(-C*UR)+B*EXP(-D*UR)+(l-A-B)*EXP(-E*UR) 

C - USE  LOG-LOG  INTERPOLATION 

LOG  I  NT ( X , X 1 , X2 , Y 1  ,  Y2 )  =  EXP ( ( ALOG ( Y2 ) -ALOG ( Y 1 ) ) * 

1  (ALOG  (X)  -  ALOG  (Xi )  >  /  ( ALOG ( X 2 )  -  ALOG (XI))  +  ALOG  (YD) 

SEMINT (X,X1,X2,Y1,Y2)  =  EXP ( (ALOG < Y2) -ALOG ( Y 1 )) * 

1  ( X-X 1 ) / ( X2-X  1  >  +  ALOG(Yl)) 

DO  NOT  LOOK  UP  CONSTANTS  IF  YOU  ARE  NOT  CHANGING  ENERGY 
IF (EKEV. EQ. EOLD) THEN 
J  =  JOLD 
GOTO  212 
END  I F 

EOLD  =  EKEV 


C 


J  =  2 
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C - use  lowest  two  data  points  for  energies  below  table 

I F  ( E  K  E  V .  L  T .  E  ( 1 ) )  50  TO  212 
233  I F (EKEV  -  E  ( J ) )  212,211,213 
213  J  »  J  +  1 

IF  <  J . LT . 25 )  GOTO  233 

C - use  highest  two  data  points  for  energies  above  table 

GOTO  212 

211  I F ( MUR . GT , 40.  )  SOTO  213 

BIJF  =  SF(A1  (J)  ,A2(J)  ,  ALPHA  1  (J)  ,ALPHA2'J>  ,  ALPHA3  ( J  >  ,  MUR ) 

0  WRITE(1,900)EKEV,3UF,J,A1 ( J ), A2 (J ), ALPHA  1 (J ) , ALPHA2  <  J  >  , ALPHA!  (J) 


SOLD  =  J 
RETURN 

212  CONTINUE 

IF (MUR, GT. 43. )  30T0  213 

X2=  BF ( A 1 ( J > , A2 ( J ) , ALPHA  1 (J) , ALPHA2  <  J ) , ALPH A3 ( J ) , MUR ) 
X1=BF(A1  (J-l )  ,A2(J-1) , ALPHA  1 ( J - 1 ) , ALPH A2 ( J - 1 ) , ALPHAS (J-l ) , MUR ) 
3UF  =  LOGINT  (EKEV , E (J-l ) ,E(J) ,  X 1 ,  X2 ) 

JOLD  =  J 
RETURN 

213  CONTINUE 

C . IF  MFPS  ARE  GT  43  INTERPRET  FROM  END  OF  RANGE 

X2=  BF <  A 1 (J) ,A2(J) , ALPHA  1 ( J ) ,ALPHA2(J) , ALPH A3 ( J ) , 40 . ) 

X1  =  BF (A1 (J-l ) ,  A2 ( J- 1 ) , ALPHA  1 (J-l) , ALPHA2 ( J- 1 ) , ALPHA3 ( J - 1 ) ,43. ) 
5UF2  =  LOGINT  (EK£V,E ( J-l ) ,E(J) ,X1 , X2) 

X2=  BF ( A 1 ( J ) ,A2(J) , ALPHA  1 ( J ) ,ALPHA2(J) ,ALPHA3(J) ,35.  ) 

X 1 =  B F ( A 1 (J-l) , A2 ( J  —  1 ) , ALPHA  1 (J-l) , ALPH A2 (J-l) , ALPHA3 ( J- 1 )  ,35.  ) 
3UF1  =  LOGINT  (EKEV, E ( J - 1 ) ,E(J) , X 1 , X  2 ) 

BUF  =  SEMINT  (MUR , 35.  ,40.  ,  BUF 1 ,  BUF 2 ) 

JOLD  =  J 
RETURN 
END 
C 
C 

FUNCTION  BCONC (EKEV , RMFP ) 


C . THIS  FUNCTION  CALCULATES  THE  BUILD  UP  FACTOR  IN  CONCRETE  OR 

C . CONCRETE  LIKE  MATERIALS  FROM  THE  NBS  SOLUTION  USING  12  PARAMETERS 

C . BY  E I SENHAUER  AND  SIMMONS  REF  NUC  SCI  AND  ENG  56,263-270,1975 

C 

C . R  IS  IN  MFPS 


C 

REAL  E  (25)  ,  DO', 25)  , MUC  (25)  , LOG INT  , MUCONC 
REAL  A 0 (25) , A 1 1 2 5 ) ,A2(25) , A 3 (25) , A 4 ( 2 5 ) 

REAL  B0(25) ,B1  (25) ,B2(25) ,33(25) , B 4 (25) 

R E A L * 8  BFC 

COMMON  /0LDJ2/  E0LD2,  J0LD2 

DATA  E/  15.  ,20. ,30. ,40. ,50. ,60. ,80. , 100.  , 150. ,200. ,300. , 

1  400. ,500. ,600. ,803. ,  1000. , 1500. ,2000. ,3000. ,4000. ,5000.  , 

2  6000. .8000.  ,10003. ,  15000. / 

DATA  MUC/9. 01 00, 3. 4450, 1. 1180,0.5588,0.3608,0.2734,0.2004,0. 1704, 

1  0.1399,0.1250,0. 1 073, 0.0958,0. 0873, 0.03U7, 0.0709, 

2  0.3637,0.0519,0.0448,0.0365,0.0319,0.0290,0.0270,0.0245, 
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3  0.0231,0.0215/ 

DATA  DO/0. 1539,0. 1497,0.  1327,0.  1233,0.  1  128, 

10.1 118,0.1215,0.1376,0.1791,0.2149,3.2693, 

2  0.3339,0.3409,3.  3674,0. 4085,0.4434,3.496  1  , 

3  0.5280,0.5679,0.5938,0.6158,0.6241,9.6425,0.6560,0.6676/ 
DATA  A3/  1.  1  1 358-2, 2. 4293E-2, 0.06645,0.  1  1734,0.  15910, 

1  0.  186  44,0.22372,0.  24651 ,0.27760,9.29820,0.31837, 

2  0.32929,0.33559,0.33965,0.34413,0.34678,0.35110, 

3  0.35445,0.35691 ,3.35425,0.34788,0.33878,0.31686, 

4  0.29369,0.24038/ 

DATA  A 1  /-0. 00 1  1 583, -0.0024  6 18, -0.0060622, -0.0092903, 

1  -0.0 15356, +0.025044, -0.01 01 44, 0.043674, 0.02 1  783, 

2  3.329893,0.012225,0. 3395543,0.3044976,0.0373127, 

3  3.3391422,3.0393983,0.013316,3.3399221,0.3078632, 

4  0.013371,0.040628,3.330839,-0.013397,-0.016610,-0.0029308/ 
DATA  A2/5.483E-04,1.2439E-03, 3. 9389E-03, 8. 63438-03, 

1  1.5169E-02, 1. l833£-02,3. 7832E-02 , 3. 2359E-02 , 2 .  9378E-02 , 

2  4.8446E-02, 

3  5.503l£-02,5.0474E-02,2.2018E-02,2.7012E-02,2.2857E-02, 

4  2. 0876E-02, 1.9785E-02, 1.9206E-02, 2.  1 3018-02,3. 9089E-02, 

5  -9.3940E-03,-8.8533E-03,2.883SE-02,2.726lE-02,-l.  9101E-02/ 
DATA  A3/2.  1503E-04.5. 0984E-04 , 1 . 7224E-03 , 4 . 0328E-03 , 

1  7.5754£-03,1.3153E-03,3.0337E-92,2.6207E-02,3.  1915E-02, 

2  B.2357E-02, 1.0486E-01 , 1.  1 029E-01 ,3. 4039E-02 , 3 . 1  744E-02, 

3  4.  41  23E-02,5.0435E-02,S.2275E-02,4.9l03E-02,3.9008E-02, 

4  2.9853E-02,2.8888E-02,2.5095E-02,1.8330E-02, 1.4079E-02, 

5  2. 2817E-02/ 

DATA  A4/4.2382E-05,1.0829E-04,3.  9905E-04,  1.0308E-03, 

1  2. 2527E-03.0. 00000000,  1.2457E-02, 3.  151  4E-06.5. 804  1 E-02, 

2  -1.  46088-04,-8.641  IE-04,-6. 3551E-08, 1 . 0539E-01  ,  9. 4382E-02, 

3  7.330?E-02,5.9498E-32,4.0714E-02,3.2172E-02,2.2187E-02, 

4  2. 0073E-03 , 6. 734GE-05 , 3. 5701E-03, 1 .  5406E-03, 3.  6397E-07 , 

5  9 . 8532E-03/ 

DATA  80/0.57471,0.56384,0.54853,0.53850,0.53787, 

1  0.54758,0.57172,0.60062,0.65913,0.69464,0.73934, 

2  0.76455,0.78166,0.79363,0.83851,0.81671 ,0.82277, 

3  0.81861,0.80467.0.79201,0.78176,0.77269,0.75938, 

4  0.75116,0.73974/ 

DATA  B2/0. 72706, 0.71 51  8, 0.6981 4, 0.68724, 0.68395, 

1  0.90637,0.69136,0.86298,0.77287,3.84634,0.82398, 

2  0.82998,0.76467,0.79645,0.80369,0.80482,0.79382, 

3  0.79546,0.81708,0.85749,0.72880,0.63341,0.85214, 

4  0.86119,0.55119/ 

DATA  B 1/0. 42393, 0.4 1368, 0.39364, 0.3601 0,0. 28648, 

1  0.70718,0.80330,0. 69185,0.68532,0.71979,0.68807, 

2  0.65925,0.56762,0.58849,0.57175,0.55269,0.52360, 

3  0.50858,0.46877,0.30734,0.83930,0.83961  ,0.57962, 

4  0.56105,0.85134/ 

DATA  B3/0.  8731 7, 0.86506, 0.85649, 0.85396, 0.85703, 

1  1.00436,0.85532,0.98931,0.87818,1.01431 ,1.01389, 

2  1.00620,0.87305,0.90888,0.91929,0.91859,0.91265, 
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3  3.90341,3.91132,3.96145,0.96338,3.94394,3.95976, 

4  3.97377,3.36292/ 

DATA  34/3.96443,3.96154,3.95995,3.96113,3.96562, 

1  1.30333,0.93334,1.17560,1.03645,1.09632,1.07239, 

2  1.34173,1.33139,3.99331,0.99455,3.99163,3.98767, 

3  3.93542,3.93436,1.01538,1.13393,1.03529,1.31131, 

4  1.22932,0.97741/ 

3rC!R,U,A0,Al,A2,A3,A4,B0,81,82,B3,34,D0)  -  1  +  ( 

1  A3*U*R*DEXP(-D3LE(U*R/80) ) +A1 * <U*R/B1 ) *#2*DEXP (-DBLE (U«R/B 1 ) )+ 

2  A2  «  (U*R/B2) **2  *  DEXP ( -D8LE (U*R/B2) )  + 

3  A3  *  ( U * R / B 3 ) **2  *  DEXP ( - D B LE < U * R / B 3 > )  + 

4  A4*(U*R/B4)**2*DEXP(-DBLE(U#R/B4)  ) )  /  ( DO  *  D  E  X  P  (-DBLE  (U*R)  )  ) 

r - USE  L02-LC3  INTERPOLATION 

103 1  NT ( X , X l , X2 , Y 1 , Y2)  *  EXP  ( !  ALOG ( Y2) -ALOG ( Y 1 ) ) * 

1  <ALQG<X!-AL03(Xi))/(AL03(X2)-AL03<Xl))  +  ALOB(YD) 

R  =  RMFP  /  MUCONC ( EKEV) 

C  WRITE!*, 1330)  'BCONC  2  ' , EKEV 

C  ***** *  DO  NOT  LOOK  UP  CONSTANTS  IF  YOU  ARE  NOT  CHANGING  ENERGY 
IF (EKEV. EC. E0LD2) THEN 
J  =  JULD2 
GOTO  212 
END  IF 

E0LD2  «  EKEV 
C 

Js2 

C - use  lowest  two  data  points  for  energies  below  table 

IF (EKEV. LT. E  ( 1 ) )  GO  TO  212 
233  IF (EKEV  -  E < J) )  212,211,210 
213  J  -  J  +  1 

IF ( J. LT. 25)  GOTO  233 

C - use  highest  two  data  points  for  energies  above  table 

GOTO  212 

211  BCONCaBFC(R,f1UC(J),.40(J)  ,  A 1  ( J ) ,  A2  ( J )  ,  A3  ( J )  ,  A4  ( J )  ,  B0  ( J )  ,  B!  ( J )  , 

1  32 (J ) , S3 ( J ) , B4  ( J ) , DO  ( J )  ) 

URCONC  ■  R  *  flUC(J) 

J0LD2  «  J 
RETURN 

212  CONTINUE 

X2  =  BFC(R,MUC(J) , A3 ( J ) ,A1 (J) , A  2 ( J ) , A3 ( J ) , A  4 ( J ) ,B0(J) ,B1 (J) ,B2(J) , 

1  33 ( J ) , B4  ( J )  ,  DO  ( J ) ) 

K  ■  J  -  1 

X  1  «BFC  (R ,  flUC  (K)  ,A0(K)  ,A1  (K)  ,A2(K)  ,A3(K)  ,  A  4  ( K )  ,B0(K)  ,B1  (K)  ,B2(K)  , 
1  B3 ( K) , B4 ( K)  ,  DO  !K ) ) 

BCONC  «  LQGINT!EKEV,E(J-1),E(J),X1,X2) 

URCONC  ■  R  *  LOGINT(EKEV,E(J-l) ,E(J) , MUC ( J  —  1 ) , MUC ( J ) ) 

J0LD2  »  J 
RETURN 
END 
C 

REAL  FUNCTION  MUCONC  (EKEV) 

C - this  fucntion  calculates  mu  over  rho  for  concrete  in  cm2/gr 
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REAL  E ( 25 ) ,  HU0RHO<25> ,  LOGINT 

C - mass  attenuation  data  from  Nuc  Sci  &  Eng  56  pr  267,  1975 

C 

DATA  E/15. ,20. ,33. ,43. ,53. ,63.  80. , 1 33.  ,  1 53.  , 230 .  ,  333.  , 

1  433. ,503. ,630. ,833. ,1303. , 1 533. , 2333 ., 3030 . ,4330. ,5330.  , 

2  6333. ,8300. ,12330.  , 15330.  / 

DATA  MUQRHO/8. 01 00, 3. 4450, 1. 1180,3.5588,3.3638,0.2734,0.2334, 

1  3.1704,3.1399,3.1250,0. 1373,3.0953,0.0873,3.3807,0.0709, 

2  0.3637, 0.05 19, 0.0443, 0.^365, 0.031 9, 0.0290, 0.0270, 0.0245, 

3  0.0231,0.0215/ 

C 

C - USE  LOG-LOG  INTERPOLATION 

LOG  I  NT ( X , X 1 , X2 , Y1 , Y2)  =  EXP ( ( ALOG (Y2) -ALOG ( Y  l ) ) * 

1  (ALOG ( X ) -ALOG (XI) )/ (ALOG ( X2) -ALOG ( X 1 ) )  +  ALOG ( Y 1 ) ) 

C 

C 

J32 

C - ui»  lowest  two  data  points  for  energies  below  table 

IF(EKEV.LT.Ed) )  GO  TO  212 
230  IF (EKEV  -  £ ( J ) )  212,211,210 
210  J  »  J  +  1 

IF ( J. LT. 26)  GOTO  200 

C-- - use  highest  two  data  points  for  energies  above  table 


GOTO  212 

211 

MUCONC  ■ 
RETURN 

MUORHO(J) 

212 

CONTINUE 

MUCONC  « 

RETURN 

END 

LOGINT (EKEV, E(J-l) ,E(J> , MUQRHO ( J- 1 ) , MUORHO  ( J ) ) 

